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HAUPTAUFSÄTZE 


Das Problem derRingwellen auf einer Flüssigkeitsoberfläche. 


Von Hans Widenbauer in München. 


Übersicht. Es werden «die Oberflächenwellen unter gewissen Einschränkungen der 
Variablen r und # berechnet und diskutiert, welche entstehen, wenn ein kreiszylindrischer 
Stein in Tiefwasser geworfen wird. Ausgegangen wird von den üblichen Gleichungen der 
Hydrodynamik und mittels des Fourierschen Integraltheorems eine einmalige Störung (dis- 
kontinuierliche Anfangsform) zur Zeit t==0 dargestellt ($ D). Die Integrationen ($ 2) werden 
unter den eingangs erwähnten Einschränkungen von r und f nach einer exakteren Methode 
(„Sattelpunktmethode*) bewerkstelligt, als es bisher bei ähnlichen Wellenproblemen geschah. 
In$ 3A wird die Oberflächenform einer auf einen Punkt konzentrierten Störung, in $ >3B 
(lie Oberflächenform bei endlichem Radius der Störung diskutiert. 


$ 1. Ableitung der Gleichungen der Ringwellen. 

Das Problem, um welches es sich handelt, ist folgendes: 

Eine unendlich ausgedehnte Flüssigkeit von einer bestimmten Tiefe, welche anfänglich 
in Ruhe ist, wird durch eine Ursache irgendwelcher Art an einer Stelle in Bewegung ge- 
setzt. Es treten Wellen auf, die konzentrisch nach außen laufen. Diese Erscheinung ist uns 
aus dem täglichen Leben ganz vertraut; sie entsteht, wenn wir etwa einen Stein in eine 
Wasserfläche werfen. Gefragt ist nun nach der Oberfläche an einer bestimmten Stelle, bzw. 
zu einer bestimmten Zeit. 

Ohne sich allzu weit von der Wirklichkeit zu entfernen, kann man in der mathematischen 
Behandlung des Problems weitgehende Vereinfachungen treffen. Eine Reibung der Flüssig- 
keit kommt nieht in Betracht; als äußere Kraft ist nur die Schwere vorhanden: die Flüssig- 
keit ist wirbellos, weil anfänglich in Ruhe; ferner werden wir sie als inkompressibel an- 
nehmen. Es existiert daher ein Geschwindigkeitspotential 9, welches der Gleichung 


genügt. 
Zur mathematischen Behandlung des Problems legen wir ein Zylinderkoordinatensystem 
(r,a, 9) zugrunde. Die Störungsstelle denken wir uns im Ursprung; die ungestörte Wasser- 
oberfläche werde dureh die {r, a}-Ebene gebildet; die positive y-Achse gehe senkrecht nach oben. 
Die zu lösende Differentialgleichung des Problems ist 
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Die Grenzbedingungen sind folgende: 
Ist die Wassertiefe, so muß in der Tiefe y=  h die Normalgeschwindigkeit 
0 
sein. 


Die Gleichung der Wasseroberfläche werde im folgenden mit = 1) bezeichnet. Dann 
haben wir als weitere Grenzbedingung, daß der Druck an der Oberfläche y =, verschwinden 
muß. Dies liefert aus den Eulerschen Gleichungen unter Beschränkung auf kleine Eır- 
hebungen , die bekannte Bedingung 


Of 


"erner muß die Geschwindigkeit, mit der die Flüssigkeitsteilchen in normaler Richtung 

fortschreiten, gleich der Normalkomponente der Geschwindigkeit der Oberfläche sein: 

04 


Of 


0 


Als Lösung des Problems erhält man bei Kreissymmetrie und periodischer Zeit- 
abhängigekeit die stehende Zylinderwelle') 


n=cosot-J, (kr) 


sinot Cosk(y-+h) 
Cos kh 


(1). 


J, (k r) 
Dabei ist die Frequenz o© mit der Wellenzahl % und der Wassertiefe h durch dieselbe 
„Dispersionsformel* 

gk-Tgikh) 
verknüpft wie bei den ebenen Wasserwellen (2-dimensional). 


Nun stellen wir die Gleichung der Ringwellen auf, welche bei einer einmaligen 
Störung der Wasseroberfläche entstehen. Wir denken uns zur Zeit ?=0 einen kreiszylin- 
drisch begrenzten Stein in das Wasser geworfen. Dieser 


EEE r Anfangszustand wird mathematisch der Forderung ent- 
RG | 7777777 Sprechen, daß zur Zeit #=-0 der Wasserspiegel in einem 
C Per u 7 Kreis vom Radius d um das Stück © gesenkt sei (Abb. 1): 
rzd 
(2). 
Abb. 1. | 


Um nun auf Grund dieser Anfangsbedingung die Gestalt der Oberfläche zu einer be- 
liebigen Zeit zu erhalten, benützen wir die Fouriersche Integraldarstellung °) 


F(r) \k-J,(kr)dk\aF (a) J, (ka) da. 


In unserem Falle ist 
(Cd 
\aFla)J,(ka)da U\aJ,(ka)da= 1; (k.d)°). 


Damit erhalten wir aus (1) 
„= —(d- \ 


Cosk(y+ h) 


gk- Tg(kh) 


'; Vgl. bei unendlieher Wassertiefe Lamb: Hydrodynamies, IV. Ed., p. 423. 
"), Riemann-Weber: Part. Diffel., VII. Aufl., Bd. 1 (1930), S. 420. 
»), Riemann-Weber: Part. Diffel., VIII Aufl., Bd. 1 (1930), S. 415. 
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Auf entsprechende Weise kann man verfahren, wenn an Stelle der Anfangserhebung 
für den Druck der Zustand (2) zur Zeit t=V0 gefordert würde. Im Gegensatz zur früheren 
Lösung (3) würde sich die Oberfläche aus sin » #-Gliedern, das Geschwindigkeitspotential aus 
eos »#-Gliedern aufbauen. Wir wollen darauf nicht weiter eingehen und uns auf den etwas 
anschaulicheren Fall einer Anfangsform der Oberfläche beschränken. 


$ 2. Integration der Ringwellen bei Tiefwasser. 
Im folgenden soll der Fall betrachtet werden, daß die Wassertiefe h groß gegen die 
Wellenlänge i( ug: ist („Tiefwasser“). Die Dispersionsformel geht dann über in 
o’—ygk. 

Die interessanteste Frage ist die nach der Gestalt der Oberfläche. Wir wollen daher 
nur die Integration der Oberfläche »; durchführen, obwohl es mit den Methoden, welche wir 
bei der Auswertung des Integrals 7 anwenden werden, nicht schwierig wäre, auch das Po- 
tential zu integrieren. 

Es handelt sich um folgendes Integral, wobei schon &’=gk gesetzt wurde: 


n=— (Cd 


/ur Auswertung schreiben wir zunächst nur die Besselsche Funktion Oter Ordnung als 
Integral 


J,(kr) = 


7T 
\ eikrecosa da 


Da %k später von O bis & integriert werden muß, führen wir aus Konvergenzgründen 
den Integrationsweg in der «-Ebene so, daß Ne ficosa} <O, was auf dem Weg von 1 über 2 
nach 3 in Abb. 2 der Fall ist. 
Führen wir nun eine neue Integrationsvariable » = ykd ein 
und setzen für cosygkt Ne Lei \ so erhalten wir 


o und o sind dabei reine Zahlen mit der Bedeutung 
r /3 d- 
= = . . . . . . . . +) 
| r? (9) 


Zur Auswertung des Integrals benützen wir eine in der 
theoretischen Physik mehrfach gebrauchte Näherungsmethode, die 
„Sattelpunktmethode*“ °). | Abb. 2. 


A. Integration nach n. Zunächst müssen wir das analytische Verhalten des Integranden 
in der komplexen »-Ebene betrachten. 


Wir wollen das Gebiet suchen, wo 
sei i(n’cosa+2on)= fin), 
nJ, (n’) = F(n), 


daher f(n) = s\ase “+2oe 


#4, Riemann-Weber: Part. Diffgl., VIII. Aufl., Bd. I (1930), S. 411. 
5) Vel. Riemann-Weber: Part. Diffgl., VIIL Aufl., Bd. 1 (1030), S. 444. 
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Die Grenzkurve, welche die Gebiete Re trennt, ist gegeben durch Re = 
oder in Polarkoordinaten 
2osiny 
asın (2r—+y)' 
Das ist demnach eine Kurve, welche 2 Asymptoten besitzt, die mit der reellen Achse die Winkel 


und einschließen. Nun ist aber f(n) nur vom 2. Grade. Die Kurve = 0 
stellt daher eine Hyperbel, und zwar, da «die Asyınptoten aufeinander senkrecht stehen, eine 


eleichseitige, dar (I in Abb. 3 und 4). Der Mittelpunkt derselben ist 


| 
+ 
A 


Geht man zum Komplexen über 
—N,, 
so lautet er 


N, : 
cos y 


Die Hyperbel geht für alle Werte von O0<y<7 dureh den Nullpunkt und berührt dort die 
reelle Achse, wie sich unmittelbar ergibt, wenn man die Gleichung der Hyperbel in recht- 
winkligen Koordinaten =x-+iy) hinschreibt. Sie artet in die Geraden #= und y=0 
aus, wenn y=-0(a=0), und in und y--0, wenn y=r(a=n). Ferner ersieht man, 
daß im Äußeren der Hyperbel, also in dem Gebiet, in dem der Mittelpunkt n, liegt, Re tf(n)y < 0 
(schraffiert in den Abb. 3 und 4). Für 0>x verschwindet dort der Integrand exponentiell. 
Nehmen wir im folgenden an, daß o >> 1, d. h. betrachten wir das Wellenbild in einer 
Entfernung, die groß gegen den Radius der Störung ist, so werden nur diejenigen Stellen des 
Integrationsweges einen wesentlichen Beitrag zum Integral nach » in (4a) liefern, wo Ne tf(n)} 
seinen größten Wert annimmt. Dies findet bekanntlich in den Punkten („Sattelpunkt*) statt 
dn 
0 und so abzuändern, daß er, falls dies überhaupt möglich ist, durch einen solchen Sattel- 
punkt hindurchführt. In unserem Fall existiert ein einziger Sattelpunkt. Er ist 


für die =0. Wir versuchen daher den Integrationsweg unter Beibehaltung der Grenzen 


COS A 
also der Mittelpunkt der Hyperbel 1. 
Um nun eine möglichst günstige Näherung zu erhalten, bestimmen wir den Integrations- 
weg durch den Sattelpunkt so, daß auf ihm Ne tfGr), am schnellsten zu- bzw. abnimmt. Diese 
Kurven erhalten wir (vgl. Anm. 5) aus der Beziehung 


= const = 


Mit Hilfe der Entwieklung der Funktion fin) um den Sattelpunkt n, 


fin) i| +(n cos al 


COS A 
eht sie über in 


il2are(n—n.) +7 + 


Das sind 2 aufeinander senkrecht stehende Gerade, welene mit der reellen Achse die Winkel 


.) 
y 


u und 75 einschließen. Von den beiden Geraden kommt natürlich nur diejenige als 
Inteerationsweg in Betracht, welche vollständig im schraffierten Gebiet liegt (IL ın Abb. 3 u. #), 

Nun müssen wir zusehen, ob wir unseren Weg, vom Nullpunkt ausgehend, über den 
Sattelpunkt nach + x legen können. Wir bestimmen zunächst die Kurve durch den Null- 


ınkt. auf der Nefftn)‘! am sehnellsten Diese ereibt sich aus 
- 


= 


was wieder eine gleiehseitige Hyperbel ist (III in Abb. 3 und 4), deren Asymptoten mit der 
JT y 


5 bilden. Der Mittelpunkt ist, ebenso wie bei 


> und f 
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der Hyperbel I, der Sattelpunkt »,. (Das muß natürlich so sein, denn das System der Hy- 
perbeln, dem I angehört, bildet die Niveaulinien, das System, dem III angehört, die Fallinien, 
während der Mittelpunkt », der beiden Hyperbelsysteme der dazugehörige Paß (Sattelpunkt) 


ist.) III berührt für alle Werte y = die imaginäre Achse und artet für y= |Rela) 


in die Geraden und y=-+ aus. 
Die Lage des Sattelpunktes »,, sowie die Lage der Gebiete und Wege ist abhängig von 
cos a=a-el?, Wir unterscheiden daher folgende 2 Fälle: 


Jt 
1. Fall. 


Durchläuft « in der «Ebene den Integrationsweg I nach 2 (Abb. 2), so wandert der 


Sattelpunkt », = in der komplexen n-Ebene von der negativ-reellen Achse im 2. Qua- 


COS A 
dranten nach der positiv-imaginären Achse. Mit ihm drehen sieh die Hyperbeln I und II 
und die Geraden II und I1*. Abb. 3 zeigt die Konfiguration für einen festen Wert von «a. 
Wir werden nun den Integrationsweg (in Abb. 3 stark ausgezeichnet) folgendermaßen führen: 


| / 
®) 
/ 


mag. Achse 


77 
/ 
/ 
Abb. 8, Abb. 4. 


Wir gehen vom Nullpunkt aus auf dem Teil der Hyperbel III, welcher vollständig im 
schraftierten Gebiet liegt; dann gehen wir etwa auf einem Kreis von genügend großem Radius 
zur positiv-reellen Achse und auf dieser schließlich ins Unendliche. Einen wesentlichen Beitrag 
zum Integral werden nur diejenigen Stellen des Weges liefern, welche auf dem Rande des 
schraffierten Gebietes liegen, denn dort findet für o>» kein exponentielles Abnehmen des 
Integranden statt. In unserem Falle ist das, und zwar für alle Werte von a, nur der Punkt 
» 0. Um den hiervon herrührenden Beitrag zu berechnen, führen wir als neue Integrations- 
variable 

n—=iy 


ein; somit ergibt sich, wenn wir für J,(n?) nur das 1. Glied der Entwicklung um den Null- 
punkt‘) setzen, 


2 
” | 
Schätzt man nun das Integral nach y% der Größe nach ab, so findet man unmittelbar, wenn 
1 
man die Integration nach « sogleich vorwegnimmt, daß |] wie —, gegen Null geht. Wir 


werden sehen, daß wir diesen Beitrag gegen den nun zu bereehnenden vernachlässigen können. 


2. Fall. 


Durchläuft « den Integrationsweg 2 nach 3, so bewegt sich der Sattelpunkt n, von der 
positiv-imaginären Achse im 1. Quadranten nach der positiv-reellen Achse, und mit ihm drehen 
sich I, II und Ill. Den Integrationsweg legen wir folgendermaßen (stark ausgezeichnet in 
Abb. 4): 


6) E. Jahnke nm. Fr. Emde: Funktionentafeln, S. M. 
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Wir gehen vom Nullpunkt wieder auf dem Teil der Hyperbel III, der vollständig im 
schraffierten Gebiet liegt, kehren dann um und gehen auf der Asymptote II durch den Sattel- 
punkt, dann z. B. auf einem Kreis von genügend großem Radius zur positiv-reellen Achse und 
auf dieser schließlich ins Unendliche. Den Beitrag vom Nullpunkt können wir in ähnlicher 
Weise wie in 1. vernachlässigen. Der Hauptbeitrag längs des ganzen Integrationsweges wird 
vom Sattelpunkt und dessen unmittelbarer Umgebung geliefert. Wir setzen für f(n) die Ent- 
wicklung (6) im Sattelpunkt in unser Integral und führen die neue Integrationsvariable 


1/, 4 s,e 4 
".s.e = 

cos a 


ein. Den Faktor n .J,(n”) in (4) können wir bei der Integration als langsam veränderlich an- 
sehen und in 1. Näherung durch den Wert am Sattelpunkt », ersetzen. Inwieweit das be- 
rechtigt ist, wird hernach untersucht werden. Somit erhalten wir aus (4): 


») I 
— + (n— no)? cos « 
„= eos« | 
— Me \ da:n.J (n cos « 4|. 
q or 0 1 
> 
£ cosa ” 
3 2 
. 
2 
2U 0 cos a 


B. Integration nach a. Aus der Gestalt des nun zu berechnenden Integrals erkennt 
man sofort, daß auch die Integration nach a längs des Weges 2 nach 3 mit Hilfe der Sattel- 
punktmethode bewerkstelligt werden kann, sofern man 00°>1 voraussetzt. (Über die Be- 
deutung dieser Voraussetzung vgl. unter $ 3.) Die Verhältnisse gestalten sich hier einfacher 
als bei der Integration nach », so daß wir uns kürzer fassen können. 

Wir setzen 


== . 

cos a 
Das Gebiet, wo NRetg(a)} <O, ist das schraffierte in Abb. 2. Als Sattelpunkt, welcher sich 

dgla) 
aus np — 0 bestimmt, kommt nur 
in Frage. Die Entwicklung von 4 (a) im Sattelpunkt unter Vernachlässigung höherer Potenzen 
lautet 
(a — 
gla)=+ill+ 


Aus der Beziehung = ergibt sich die Richtung, mit der wir in den 
Sattelpunkt «== hineingehen müssen. Sie schließt mit der positiv-reellen Achse den Winkel 


-j ein. Zur Ausführung der Integration setzen wir 


a—ı= —b-e 


und erhalten für das Integral in (7a) 


2 

Daher nach (Ta) 

„= 


Gl. (8) gilt zunächst für o>1 und 00°>1. Nun müssen wir noch überlegen, unter 
welchen Einschränkungen die Annahme gerechtfertigt ist, daß n .J, (n?) bei der Integration 


0? 
J, 2 
A 


nach » und 5, bei der Integration nach a in den betreffenden Integrationsintervallen 


cos a 
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| 


\angsam veränderlich sind. Wir haben uns außerdem die Frage vorzulegen: 


Wie ändert 


sich 7, wenn diese zuletzt genannten Annahmen zum Teil fallengelassen werden? Zu diesem 


Zwecke berechnen wir die nächste Näherung. 


A. Integration nach n. Wir entwickeln n J, (n’) = F(n) um den Sattelpunkt », in 


eine Reihe und gehen bis zu quadratischen Gliedern. 


F(n)=F(n)+(n n)F(n,)+ F" (n,). 


Dabei bedeuten 
F(n) =mJ, 


F' (n,) 2 N. (n,°) J, (n,?) 


| 
0 
cos a’ 
n—n,=4a Um 
COS a 
Das Integral (7) geht dann über in 
cos « s..e 4 s...e 
7 Ne d a | d (n,) I (n,)+ > cosa pr (n,). 
cos a cos a 


Die Ausführung der Integration ergibt 


06 cos a) eos 
n= Ida —; e 
„T lo Io 
cos a 
0° 0? 
3 
| COS [24 [74 A 00? 
Ne da E COS 


(9). 


B. Integration nach a. Das 1. Integral in Gl. (9) berechnen wir in derselben Weise 


wie unter A. Wir setzen also: 


cos a 
G(a) = Gla)+la— (a) + — (a): 
J, (0°) 
i 7 
a—n=—be 
Die Integration liefert 
Rei (0°) e Ne \ 2J, (0°) — e 
2 
v2, (0°) cos 0 0° 2.J,(0 sin 0 0°, 


308 
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Im 2. Summanden der Gl. (9) ersetzen wir den Klammerausdruck bei der Integration 
nach « dureh seinen Wert an der Stelle «= 7, denn wir wollen in unserer Näherung nur bis 


1 
zu Gliedern mit dem Faktor ,„ gehen. Wir erhalten 


0 
v-e 0° 
» 2 | 
| 0° (0°) 20°J, (0*) sın 00°, 
Somit ist in 2. Näherung 
(0°) cos o 0° 
(10). 
+20°.J, (0°) > J,(6°)| sin 0 0? 


Nimmt man nun an, daß auch 0 >> 0°, so kann man im allgemeinen den 2. Summanden 
der Gl. (10) gegen den 1. streichen und kommt auf Gl. (S) zurück. Diese Vernachlässigung 


ist Jedoch nieht mehr gerechtfertigt, wenn .J, (0°) von der Größenordnung — und kleiner wird 


(d.h. wenn 0° gegen eine Nullstelie von ./, (0°) rückt); denn in der Nähe einer Nullstelle von 


J, nimmt maximale Werte an | 


Zum Schlusse dieses Paragraphen geben wir Gl. (5) an, wenn man für o und o ihre 
Werte nach .Gl. (5) einsetzt, 


d gd-F 
= {m COS 
(11). 
d Ir 4 4 r? 


$ 3. Diskussion. 

A. Zunächst wollen wir den speziellen Fall, der fast allein in der Literatur diskutiert 
wird, betrachten, daß der Radius d der Störung gegen Null geht. Aus Gl. (11) ersehen wir, 
daß wenn d>V0O auch 70. Um überhaupt eine merkliche Oberflächenerhebung zu erhalten, 
müssen wir fordern, daß das Volumen V der anfänglichen Störung während des Grenzübergangs 
konstant bleibt: 


Da lim 
so geht, wenn d>VÖ, (11) über in 


(12) gilt (vel. Gl. (ID) unter der Bedingung T- Il. Man überzeugt sich leicht, daß man 


im Falle d>V0 die in (10) enthaltene zweite Näherung unberücksichtigt lassen kann. Gl]. (12) 
wurde sehon von H. Lamb’) nach der Kelvinschen „Methode der stationären Phase“ 
abgeleitet, die als physikalisch einleuchtender Ersatz der Sattelpunktmethode angesehen 
werden kann. 

Über die vorbereitenden Arbeiten von Cauchy‘°) und Poisson®) vgl. H. Burkhardt '"). 


Das entsprechende 2-dimensionale Problem, wo die Oberflächenerhebung für große pr 


lautet, wurde von H. Lamb vollständig für alle Werte von f und = in dem eben genannten 
Werke gelöst. 


MH. Lamb: Hydrodynamies, IV. Ed., p. #22. 

sS) Mein. des Sav. EKtran., t. T (IST). 

») Mem. de Acad. Roy. d. Se., t. I (1816). 

10) Jahresber. d. deu'sch. Mathematiker-Vereinigung, X, ?, 1 (1908), S. 429 IT. 
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Wir kehren nun zur Formel (12) zurück und betrachten zunächst die Erhebung der 
Oberfläche an einem festen Orte r im Laufe der Zeit. Abb. 5 stellt die Kurve „= n(f) dar, 


Die Erhebungen erfolgen mit wachsender Zeit immer schneller, während die Amplituden be- 
ständig zunehmen. Aus der Wahl der Maßstäbe der Achsen in Abb. 5 ergibt sich, daß für 
verschiedene Orte die gegenseitigen Abstände zweier entsprechender Phasen sich wie yr, 
die entsprechenden Erhebungen dagegen sich umgekehrt wie r? verhalten. 


und zwar etwa von denjenigen Werten von f an, von denen ab unsere Formel gilt | » 1) 
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Abb. 6 stellt das Wellenprofil in einem gegebenen Augenblick dar, hauptsächlich für 
den Bereich der Wellen, für den (12) gilt. Für verschiedene Zeiten verhalten sieh die hori- 
zontalen Abstände zweier entsprechender Wellenpunkte wie #, die entsprechenden Erhebungen 
umgekehrt wie #*. Daraus geht hervor, daß das Wellenprofil im Laufe der Zeit auseinander- 
gezogen wird, d. h. daß die Wellenlänge einer bestimmten Welle beim Fortschreiten sich 
immer mehr vergrößert. 

Der Umstand, daß die Amplituden der Wellen mit wachsendem f bzw. mit abnehmendem r 
ins Unendliche wachsen, ist natürlich nur eine Folge des Anfangszustandes lim d>0. Wenn 
wir hernach die Störung bei endlichem Radius d betrachten, wird diese Unendlichkeit der 
Amplituden selbstverständlich verschwinden. 

Falls der Grenzübergang d>V0 schon vor den Integrationen ausgeführt wird, läßt sich 
ohne besondere Schwierigkeit das Integral ») durch Reihenentwicklung lösen. So erhielt Lamb'') 


gt 
folgende nach Potenzen von = fortschreitende Reihe: 


>2ar 65! r 10! 

Da die Reihe nur sehr langsam konvergiert, ist sie für die Diskussion nicht zu gebrauchen. 
Jedoch auf einen Umstand, der sich aus dieser Reihe ergibt, wollen wir noch hinweisen. 
Wir sehen nämlich, daß sofort, nachdem die Störung im Nullpunkt erfolgt, die Oberfläche 
in beliebig großer Entfernnng vom Störungszentrum zu sinken beginnt. Es existiert daher 
überhaupt keine Wellenfront. (Der Anfang der Welle (entweder f klein oder r groß) kann 
aus Gl. (12) natürlich nieht erhalten werden.) Die erste Vertiefung des Wasserspiegels 
pflanzt sich mit unendlicher Geschwindigkeit fort. Diese Tatsache, auf die zuerst im 2-dimen- 
sionalen Falle Lord Rayleigh') und B. Pidduek'*) hinwiesen, ist nicht besonders über- 
'aschend, da ja die Flüssigkeit als inkompressibel angenommen wurde. 


11) ].e. Anm. 7. 

12) On the instantaneous propagation- of disturbance in a dispersive medium, exemplified by waves on water 
deep and shallow, Phil. Mag. Ser. VI, vol. 18 (1909), p. >. 

13) On the propagation of a disturbance in a fluid under gravitv, Proe. Roy. Soe. Ser. A. 83 (1010), p. 347. — The 
wave-problem of Cauchy and Poisson for finite depth and slightly ecompressible fluid, Proe. Roy. Soe., Ser. A. 86 (1912), p. 306. 
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Setzen wir ın (12) die Phase 


tr 


const, 


so ergibt sich die Phasengeschwindigkeit » 


dr gt ı/gr 
Durch Klimination von e folgt 
Sie ist demnach vom Ort und von der Zeit abhängig. Als Phasenbeschleunigung ergibt sich 
Const. 


Jeder Wert der Phase eilt demnach mit einer konstanten Beschleunigung radial nach außen. 
Die Periode r und die Wellenlänge 4 der Wellenbewegung bestimmen sich aus den 
beiden Gleichungen 


Unter Benutzung von ; = I finden wir 


(was schon Cauehy'’) und Poisson') angaben). 
Die Phasengeschwindigkeit ist in Übereinstimmung mit Gl. (13a) wieder 


Führen wir in » die Wellenlänge 4 nach Gl. (14) ein, so erhalten wir 


an’ 


| 
also dieselbe Formel wie bei den ebenen Wasserwellen. Daher wird auch die Gruppen- 


. . . . . .. 
geschwindigkeit in beiden Fällen gleich, nämlich —. 


B. Nun wollen wir die Oberfläche », diskutieren, wenn der Radius d der Störung eine 
endliche Größe hat. (Siehe Ende $ 2, Gl. (1D)) 

Zum Vergleich geben wir die Erhebung im entsprechenden 2-dimensionalen Problem an, 
das verhältnismäßig frühe von W. Burnside'*) gelöst wurde. Er fand im wesentlichen 
unter denselben Voraussetzungen wie oben 


“ür solehe Werte von r, f und d, für die “1, reduziert sich (11) auf (12). Auf 


4 
diesen Fall brauchen wir daher nieht nochmals eingehen. 

Da der cos-Faktor in (11) sich verhältnismäßig schnell ändert, können wir den Faktor 
| 
‚| 
änderlich ist. Im folgenden werden wir auf die Amplitude unser besonderes Augenmerk 
richten, während wir auf die Diskussion des cos-Faktors, welche bereits in A geschah, nicht 
weiter einzugehen brauchen. 

Zunächst betrachten wir wieder die Erhebung der Oberfläche an einer festen Stelle r 
im Laufe der Zeit. Eine so bequeme und allgemeine Darstellung der Kurve =») (f) wie 
bei A ist hier nieht möglich. Wir müssen uns begnügen, für einen bestimmten Wert von 
r, d und € die Kurve zu zeichnen. Dies ist in Abb. 7 geschehen. (Dabei müßte man eigentlich 
in der Nähe der Nullstellen von J, auf Gl. (10) zurückgehen; doch spielt die Verbesserung 


als eine Amplitude deuten, die allerdings sowohl zeitlich wie räumlich ver- 


it) Cauchy:p. 24. le. Anm. 8. 
15) Poisson:p. 163. Anm. 
16) On deep-water waves resulting from a limited original disturbanee, Proe. Lond. Math. Soc. 20 (1889), p. 22. 
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ei dem gewählten Maßstab der Abb. 7 (und auch der Abb. 5) keine bedeutende Rolle, so 
laß wir sie außer acht gelassen haben.) Die Amplitude der Wellenbewegung bleibt überall 
»ndlieh und ist selbst wieder wellenförmig. 

Die Erhebung an der festen Stelle r ist ein Maximum, wenn 


wo @, die Argumente der aufeinanderfolgenden Maxima und Minima der Besselschen 
Funktion 1. Ordnung sind (und zwar ist angenähert a, 21.54: a,= 9,975: a, = 8,535). Die 
Ordinaten der aufeinanderfolgenden Maxima der Amplituden (bei festem r) nehmen mit wachsen- 


qar 


ler Zeit t ab, und zwar für große T fr angenähert wie letzteres ergibt sich, 


wenn man in (11) für /, die bekannte asymptotische Entwicklung setzt: 


-sin 


fer COS (15). 


yyda 
| für r-245 m 

| 0,3535 m | 

7 IM, 


| 
+ 1 } 444 


| 


Abb. 7. 

07 | N N 
N 
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Abb. 8. 
Die Zeit, welche das n!® Maximum braucht, um an die Stelle » zu gelangen, ist 
Ya 
Weiter ergibt sich, daß die Stellen, für welche das Maximum im Laufe der Zeit erreicht 


wird, sich mit der konstanten Geschwindigkeit 


2 


nach außen bewegen, und zwar das 1. (und relativ größte) Maximum mit der größten, die 
nachfolgenden mit immer kleinerer Geschwindigkeit. 
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Wir betrachten andererseits die Nullstellen J,. Diesen entsprechen nach der korrigierten 
(Gl. (10) Stellen zwar nieht von vollständiger Ruhe, aber von nur geringer Bewegung. Auch 
(diese pflanzen sich mit konstanter Gesehwindigkeit fort. Würde man daher die Wasserfläche 
von oben betrachten, so sähe man außer den stark bewegten Zonen Ringe nahezu glatten 
Wassers, welche mit konstanter, aber unter sieh verschiedener Geschwindigkeit fortlaufen. 
Wenn V, die Geschwindigkeit des ie" dieser Ringe bedeutet, so gilt entsprechend Gl. (16) 


| | 


[b; sind die Argumente der aufeinanderfolgenden Nullstellen von J,.] 


Abb. 5 stellt das Wellenprofil zur Zeit ?(=Ssee) dar. Die Maxima der Amplituden 
rücken bei abnehmendem r näher zusammen, während die Ordinaten derselben nahezu konstant 
bleiben (vgl. Gl. (15)). 


Schließlich wollen wir untersuchen, welchen Einfluß der Radius d und die Tiefe C der 
Störung auf die Wellenform ausüben. 


Die Ordinaten der Oberflächenerhebung sind nach Gl. (11) direkt proportional €, was 
physikalisch verständlich ist. Eine weitere Wirkung übt die Tiefe © der anfänglichen Störung 
auf die Wellenbewegung nicht aus. 


Bemerkenswert ist, daß im Bereiche der Näherung (11) der Radius d der Störung keinen 
direkten Einfluß auf die Wellenlänge und Periode der Wellen (Gl. (14)) hat, sondern nur auf 
die Amplitude, und zwar sind für große Er; die Abstände aufeinanderfolgender Maxima 
(oder auch aufeinanderfolgender Nullstellen) der Amplitude sowohl bei festem r wie auch bei 
festem F umgekehrt proportional zu d (vgl. Gl. (15)). ‚Jedoch übt d einen indirekten Einfluß 
auf die Wellenlänge A und die Periode 7 der Wellenbewegung insofern aus, als die Amplitude, 
die ja von d abhängig ist, diejenigen Wellenlängen und Perioden, welche in der Nähe der 
Maxima der Amplitude liegen, hei vorhebt, diejenigen Wellenlängen und Perioden aber, welche 
in der Nähe der Nullstellen liegen, unterdrückt. So erhalten wir aus den Gl. (14) und (16) 
die Wellenlänge 7, und Periode 7, der Welle in der Nähe des n!"n Maximums, welches an 
einer festen Stelle r erreicht wird, 


d 


a „ 


Demnach ist ZA, direkt proportional mit dem Radius der Störung, 7, direkt proportional mit 
| d, die Phasengeschwindigkeit der Welle in der Nähe des n'!*" Maximums infolgedessen 
direkt proportional mit Yd. 


Schlußbemerkung. Obwohl die Ringwellen eine so alltägliche Erscheinung sind, hat man 
sieh in der theoretischen Physik noch wenig damit beschäftigt. Das 2-dimensionale Problem 
(Fortpflanzung in geradlinigen Wellenzügen) kann «durch die Arbeiten von Burnside') und 
Lkamb'‘) als gelöst betrachtet werden sowohl für Anfangsstörungen von endlicher Breite als 
für solehe, die auf eine Linie konzentriert sind. Dagegen hat man sich bei dem 3-dimensio- 
nalen Problem (Fortpflanzung in Ringen) meist auf Anfangsstörungen beschränkt, die auf einen 
Punkt konzentriert sind ($ 3A), so daß die bemerkenswerten Eigenschaften der Amplitude, 
wie wir sie im $ BB dargelegt haben, nicht erhalten werden konnten. Überdies geschah die 
Bereehnung der Integrale auf die mathematisch nieht ganz befriedigende „Methode der statio- 
nären Phase“). In anderen Arbeiten’) wurde angenommen, daß die Störung zur Zeit 0 
in gewisser Weise über die ganze Flüssigkeit ausgebreitet ist (kontinuierliche Anfangsform), 
was die mathematische Behandlung vereinfacht, aber der naturgemäßen Formulierung des 
Problems Abbruch tut. 


Zum Schluß möchte ich meinem verehrten Lehrer, Herrn Geheimrat Prof. A. Sommer- 
feld, für wertvolle Anregungen und bereitwillige Hilfe herzlich danken, desgleichen Herrn 
Dr. W. A. Maue für Besprechungen über die Sattelpunktmethode. 435 


17\ Burnside: |]. e Anm. 16. 
18) ,amb: |]. e. Anm. 7. 


19), Lamb: 1. e. Anm. Havelock: The propagation of groups of waves in dispersive media, Proc. Roy. 
Soe. 81 (1008), p. 308. 

>»), K. Terazawa: On deep-sea water waves eaused by a local disturbance, Proe. Rey. Soe. Lond. Ser. X 
97. On the oseillations of deepsea surface caused by a local disturbance, Science Rep. Töhoku Univ. 
p. 109. 
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Unterteilung des Tafelschritts. 


(„Untertafelung“ „Briggs-Moutonsche Aufgabe“.) 
Von F. Emde in Stuttgart, 


man nicht sämtliche Werte der Tafel nach der Definition der zu tafelnden Funktion be- 
»echnet, sondern nur einige gleichmäßig über die ganze Tafel verteilte Werte („Ausgangs- 
werte“) und darauf die dann noch fehlenden aus den schon berechneten erschließt. Man be- 
‚echnet also die Tafel, kurz gesagt, erst nur in weiten Schritten und geht dann durch Inter- 
polation zu feineren Schritten über, z. B. zu 10mal so feinen. Dafür stehen verschiedene 
Verfahren zu Gebote. 


B“ der Berechnung einer Zahlentafel erspart man sich gewöhnlich viel Rechenarbeit, wenn 


I. 

Nach Lagrange kann man durch die Punkte, die die berechneten Werte darstellen, 
oder durch einige von ihnen eine Parabel höherer Ordnung legen. Man ersetzt also die Funktion 
dureh ein Polynom. Die Vorzahlen dieses Polynoms hängen linear von den schon berechneten 
Werten (den Ausgangswerten) ab. Die wirkliche Ausführung dieses einfachen Gedankens ist 
meist sehr mühsam. Nun braucht man aber die Funktionswerte, wenigstens für die Tafel, 
nieht für beliebige Werte der unabhängigen Veränderlichen, sondern nur für wenige, über 
das grobe Intervall gleichmäßig verteilte Werte. Setzt man diese Werte ein, läßt aber die 
Ausgangswerte unbestimmt, so erhält man für jeden Teilpunkt ein für allemal eine lineare 
Funktion der Ausgangswerte mit numerischen Koeffizienten, die sich auf beliebige zu tafelnde 
Funktionen anwenden läßt. Legt man z. B. die Parabel durch 6 äquidistante Punkte y_,...9%,, 


so erhält man für den Funktionswert in der Mitte zwischen ihnen den Ausdruck (17, p. 45) 


75 125 y_,+9% , 5y_,+ 9%; 
— 1>8 128 10 64 100 


Bei neu zu bereehnenden Tafeln so zu interpolieren, empfiehlt K. Pearson (LT, p. 22). 
Man erspare es sich so, die Differenzenstaffel für die erst grob geteilte Tafel zu bilden. Zur 
Sicherung des Ergebnisses gegen Rechenfehler genüge es, erst nach der Verfeinerung des 
Tafelschritts die Differenzenstaffel zu bilden. 

Gegen diesen Rat wendet sich L. J. Comrie (L12). Wenn. man so verfahre, könnten 
beträchtliche Fehler in den Ausgangswerten unbemerkt bleiben, weil sie bei der Verfeinerung 
des Tafelschritts die neuen Nachbarwerte in Mitleidenschaft ziehen und daher im Gang der 
Differenzen nicht zum Vorschein kommen. Man müsse Produkte aus vielziffrigen Zahlen bilden. 
Bei vielstelligen Tafeln reichten dazu die Einstellwerke der Rechenmaschinen nieht aus. Die 
aufeinanderfolgenden Produkte seien von derselben Größenordnung statt abzunehmen. Das 
von Pearson empfohlene Verfahren sei sehr mühsam. Aus allen diesen Gründen sei es zu 
verwerfen. (Siehe auch L 17, p. 559.) 


Wenn man schon für die grob gestufte Tafel die Differenzenstaffel gebildet hat, so kann 
man die Vorzahlen des Ersatzpolynoms durch die Differenzen statt durch die Ausgangswerte 
selbst ausdrücken. Von jeder Ordnung kommt in dem Polynom je eine Differenz vor (bei 
Mittelwertbildungen zwei). Man gelangt so zu den bekannten Interpolationsformeln von 
Bessel und von Stirling. (Die Newtonsche kommt praktisch nur am Anfang und Ende 
der Tafel in Betracht.) Setzt man an jedem Teilpunkt den Wert der unabhängigen Veränder- 
lichen ein und läßt die Ausgangsdifferenzen unbestimmt, so erscheinen die Funktionswerte an 
den Teilpunkten als lineare Funktionen von Differenzen verschiedener Ordnung mit nume- 
rischen Koeffizienten. Diese Koeffizienten findet man in den Interpolationstafeln zusammen- 
gestellt, z. B. am Schluß der Bücher von Rice (L4) und Lindow (L 10). Siehe auch L 14, 
p. 110... 125. Ein Beispiel ist die bekannte Formel für die Interpolation in die Mitte: 


Hierbei sind die groben Differenzen in der folgenden Weise bezeichnet: 


1,0 
(05) 
+00 Y, II, 
(+05) 


1,0 Y, II, 


2 
| 
| 
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Wie ersichtlich, kann man hier wenigstens erwarten, daß im allgemeinen jedes folgende 
Produkt klein gegen das vorhergehende ist. Wenn man einer fertigen Tafel einen vereinzelten 
/Z;wischenwert entnehmen will, wird man in der angegebenen Weise verfahren '); nicht aber, 
wenn man eine neu zu bereehnende Tafel feiner stufen und die Stufen nicht zerade nur halh 
so groß machen will. Denn dafür hat man bequemere Wege. 

Wenn man von der groben Teilung der Tafel zu einer mehr als 1Omal so feinen über- 
gehen will, so interpoliert man ein oder mehrere Male in die Mitte und geht erst dann an 
die Fünf- oder Zehn-Teilung. Die Zahl der Ausgangswerte kann man natürlich auch nicht 
beliebig verringern. Comrie (L12, p.510) rät dazu, die Ausgangswerte so dieht zu wählen, 
daß die siebenten Differenzen nicht 1000 Einheiten der letzten Stelle übersteizen. Sonst 
würde es schwierig, mittels der Differenzenstaffel Rechenfehler aufzudecken. 


II. 


Denken wir uns für die feiner gestufte Tafel die Differenzenstaffel gebildet. Wenn uns 
für alle Differenzen irgendeiner Ordnung, sagen wir der n-ten, die Zahlenwerte gegeben werden, 
von den Differenzen niedrigerer Ordnung aber nur Je ein Zahlenwert, so können wir die ganze 
Differenzenstaffel und die Funktionswerte rekonstruieren. Denn wir brauchen nur zu dem 
einen gegebenen Differenzenwert (n  1)-ter Ordnung die Werte »-ter nach und nach zuzu- 
schlagen, um alle Differenzenwerte (n — 1)-ter Ordnung zu erhalten. Dann kann man ebenso 
alle Differenzen (n — 2)-ter Ordnung herstellen usw.. und schließlich die Funktionswerte an 
den Teilpunkten des groben Intervalls. Ein besonders einfacher Fall ist der, daß die Diffe- 
renzen »-ter Ordnung konstant sind, also die Differenzen (n + 1)-ter Ordnung =. 

H. Briggs und G. Mouton haben sich daher die Aufgabe gestellt, aus den Diffe- 
renzen der groben Staffel je einen Differenzenwert jeder Ordnung der fei- 
nen Staffel zu bestimmen. Diese Aufgabe ist aber erst von Lagrange allgemein ge- 
löst worden (L1)?). 

Die Anweisung zur Berechnung je einer feinen Differenz aus den groben für Zehnteilung 

des groben Schritts auf Grund der Besselschen Formel ist im folgenden tabellarisch zu- 
sammengestellt. 
Bei der fortlaufenden Summation können kleine Fehler, z. B. die Abrundungsfehler, bedeutend 
anwachsen. Um das Ergebnis gegen solehe Fehler zu sichern, muß man mit 2 bis 3 Stellen 
mehr rechnen, als im Ergebnis verlangt werden („Schutz-Zusatzstellen*). Also muß man auch 
die Ausgangswerte genauer ausrechnen, als es sonst nötig wäre. Die Zusatzstellen bean- 
spruchen natürlich die Aufnahmefähigkeit der Rechenmaschinen und können deshalb bei viel- 
stellizen Tafeln zu Umständlichkeiten führen. 

Für jedes grobe Intervall muß von Jeder Ordnung eine Differenz aufs neue bestimmt werden. 

Wenn man eine Rechenmaschine benutzt, ist die Rechnung trotzdem noch bequem, be- 
sonders, wenn zwei Personen zusammenarbeiten, von denen die eine die Maschine bedient 
und diktiert und die andere schreibt. 

Comrie zeigt (L11, p. 447), daß für die fortlaufende gleichzeitige Summierung von zwei 
Differenzen besonders die Rechenmaschine „Brunsviga Dupla* geeignet ist. Für jemand, 
der solehe Rechnungen bisher nur von Hand oder mit zwei getrennten Maschinen oder in 
zwei Reehnungsgängen mit einer Maschine ausgeführt habe, sei die Einfachheit des Rechnungs- 
ablaufs auf der Brunsviga Dupla wirklich faszinierend. Wesentlich ist dabei, daß die Maschine 
zwei Empfangswerke und eine Einrichtung hat, um Zahlen aus dem Empfangswerk in das 
Stellwerk zu werfen. 

1) Vielstellige Tafeln kann man nicht gut für lineare Interpolation einrichten, weil sie sonst zu umfangreich 
werden würden. Um dennoeh aueh bei vielstelligen Tafeln die Interpolation noch einigermaßen bequem zu machen, 
wird in England nach einem Vorsehlag von K. Pearson neben die Funktionswerte die zweite und, wenn nötig, auch 
noch die vierte und die sechste Differenz gedruckt. Einen Zwischenwert erhält man dann nach der Laplace-Everett 


schen Interpolationsformel <=1- X) 


wobei man vorteilhaft eine Interpolationstafel benutzen wird, z. B. die von E. Chappell (London, #1 Westeombe Park 
Road, 1020) oder die von A. J. Thompson (Traets for computers No. V, London 1921, Cambridge University Press). 
Siehe auch L 14, p. 126... 129 und L 17, p. 147. 

Ohne Differenzen sollten vielstellige Tafeln nieht mehr gedruckt werden. 

?) Die Ableitung der Formeln auf Grund der Besselschen Interpolationsformel gibt Rice (L4, p.83...91), auf 
Grund der Stirlingschen Radan (L3, p. 21). Die Zahlenwerte für die Zwei-, Fünf- und Zehn-Teilung des Tafel- 
sehritts auf Grund der G außsehen Formel findet man in der bekannten fünfstelligen Logarithmentafel von F.G.Gaub 


(1. 2). einige Zahlenwerte für Zehnteilung nach Bessel und Stirling bei Milne-Thomson (L 15, p. 87). Thiele (L &, 


S. 00 ...9) und Steffensen (L 9, p. 7l... 79 drücken umgekehrt die groben Differenzen durch die feinen aus. 
Banuschinger (L », Whittaker-Robinson (LS, p.53...5N) und Davis (L 14, p. 88) geben nur die Ableitung auf 


(rund der Newtonschen Formel. 
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Berechnung je einer feinen Differenz aus den groben bei Zehnteilung des Tafelschritts 
(Besselsche Interpolationsformel). 
Differenzen in oder symmetrisch zur Mitte des groben Intervalls. 


> 
10 
I <5 
1000 III < 120 
Vor a 
10000 448 1,225 < 110 
Files x | — 67 23,0 0,6 10,01 FI < 7500 
<5 
IV’ < 240 
= IV 0.50 = < 1100 
Hot 1.05 0 0 Einheiten der 
100 letzten Stelle 
+ 27, | - 
31T: +0) 
Ä 10.000 x 
| 1.22,36125 1,44 +0,05 
100 000 


Nach einer freundlichen brieflichen Mitteilung Herrn Comries kann man mit den (nach 
’atenten von Ellis gebauten) „National*- Buchungs-Maschinen (Klasse 3000) sechs Differenzen 
eleichzeitig addieren. Die Untertafelung sei durch diese Maschinen vollkommen revolutioniert 
worden. (Siehe auch L 16.) Zusatz bei der Korrektur: Mit den Multiplex-Buchungsmaschinen 
der Astra-Werke (Chemnitz) kann man wahrscheinlich noch höhere Differenzen ebenso 
addieren. Solche Maschinen kosten rund 6000 RM. Siehe Lenz, Die Rechen- und Buchungs- 
maschinen, 3. Aufl., Leipzig 1932 bei Teubner, S. 55 und 2. | 

Wir wollen das Verfahren an einem Zahlenbeispiel veranschaulichen. Die folgende 
Zahlentafel 1 gibt S-stellig äquidistante Werte der Funktion’) (L 18) 


(0,75 
ia di= 025! n 
0 
und die Differenzen bis zur 7. Ordnung. Aus diesen Zahlen soll f für = 
bestimmt werden. Aus den groben Differenzen, die in der Zeile 0,95 stehen, berechnen wir 
nach der eben gegebenen tabellarischen Anweisung die feinen Differenzen, die in der Zahlen- 
tafel 2 kursiv gedruckt sind. Die 5. und die folgenden feinen Differenzen verschwinden. 


3) Bei K. Pearson: Tables of the ineomplete Z/-funetion (London 1922, H. M. Stationery Office) findet man 
‘-stellige Zahlenwerte für unsere Funktion in der letzten (mit »y = — 4,75 übersehriebenen) Spalte der linken Seiten II8... 
126, wobei ?x?= u gesetzt ist: 

0,75) 
fürn —0,10,2....27,0 oder x—=0,47...1,92. Für noch größere x unterscheidet sieh die Funktion von I um weniger als 10 7. 
Allgemein ist 


\ in dt — 
) n n 
0 
und bei Pearson mit \ n x" und Pp=-— 
) 
Yn= — "dt=I(up)-I (\ 
l 7 
mit Zahlentafeln auf S. II8...138S mit 7 Dezimalen für 
p = — 1,0 0.95 — 0, — 0,00 


und auf S. 16% ...164 mit 5 Dezimalen für v—=0,1...6,0 und 
p= — 1,00 — 1.76 0,75 
= 4,17 4 


(L 6, 
aus. 
auf 
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Darauf erhalten wir durch Aufsummieren von rechts nach links die gesuchten 9 Werte von 
f(#). Der Vergleich mit dem gegebenen Wert von f (1.0) zeigt, daß wir f (1,0) um 6 Einheiten 
der 8. Dezimale zu klein erhalten haben. Ferner ist f(0,95) = 0,909 725 33. Unser Wert ist 
also um 43 Einheiten der 8. Dezimale zu klein. Von den S Dezimalen können wir nur 6 behalten. 


Zahlentafel 


ILS 0.817 812 69 14 883 22 15355 66 243 49 
65 388 31 1 627 23 + 703 78 433 95 
00 201 00 I14 266 17 16 510 45 165m > 089 44 + 12 677 44 20451 
IS STT 86 | 22 58 6 462 21 26 34 155 4 16056 
1.0 0,032 078 86 SI 707 45 | 16 048 24 + 2 115 78 S38 00 
32 829 62 2 5779 — 1166 
1.1 13 470.25 1 504 12 


Zahlentafel 2. 


0.00 1.883 201 60 


1 S41 167 48 

> 478 383 
2 0.504 315 14 167 8 

> 305 50 12 
3 0.890 620 64 167 95 

5 1375 11 
| 0,004 758 19 167 84 

4 969 71 - 34 
0.0909 727 90 50 23 

4 802 21 - 57 
0.014 550 11 166 93 

I 635 28 + 80 
0.0919 165 39 166 13 

4 4609 15 —- 108 
S 0.023 694 54 165 10 

t 304 05 + 126 
0.927 938 59 163 84 

140 21 

1.00 0.032 078 
IV. 


In sinnreicher Weise benutzt Comrie (L12 und L 15, p. 91) die beiden unter II und III 
eeschilderten Verfahren nebeneinander, um noch bequemer zu den einzuschaltenden Funktions- 
werten zu gelangen. Der Hauptnachteil des Briggs-Moutonschen Verfahrens III besteht darin, 
daß man zum Schutz gegen Fehlerhäufungen die Rechnung mit Zusatzstellen belasten muß. Dieser 
Nachteil wäre beseitigt, wenn es gelänge, in den gesuchten einzuschaltenden Werten die noch 
vewünschte letzte Ziffer rechts für sich vorauszubestimmen und dadureh von der fortlaufenden 
Summation unabhängig zu machen. Ohne fortlaufende Summation arbeitet das Verfahren II, 
(las jedoch als Ganzes viel mühsamer als III ist. Kann man nach II die letzte Ziffer rechts 
für sieh allein finden, ohne zugleich auch alle übrigen mitberechnen zu müssen? 

Der einzuschaltende Wert 7, bestimmt sieh aus den groben Differenzen n\. nach der 
Besselschen Interpolationsformel in der Form 


n 


a] 


also als Summe von » Gliedern. Zur letzten Ziffer von 9, wird im allgemeinen jedes dieser 
» Glieder beitragen. Der Beitrag eines Gliedes wird durch seine Einer und Zehntel hin- 
reichend genau bestimmt sein, also durch zwei Ziffern. Für diese zwei Ziffern hat Comrie 
Tafeln bis » = 4 berechnet und im „Nautical Almanac for 1951” (L13, p.536...859) ver- 
öffentlicht (auch als Sonderdruek käuflich). Sie gelten für Zehnteilung oder Fünfteilung des 
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»roben Intervalls.. Man geht in die vier Tafeln mit Dy, und v=0,1.0,2|...!09 ein. Diese 
Tafeln reichen auch dann noch aus, wenn man noch die fünfte und sechste Differenz mit- 
nimmt. Die fünfte Differenz wird nämlich durch eine kleine Änderung der dritten und die 
sechste dureh eine kleine Änderung der vierten berücksichtigt. Von der dritten Differenz ist der 
neunte Teil der fünften abzuziehen und von der vierten Differenz der fünfte Teil der sechsten. 

Hat man so die letzten Ziffern der einzuschaltenden Werte festgelegt, so bildet man 
aus ihnen in der gewöhnlichen Weise die Differenzenstaffel und erhält so die letzten Ziffern 
der Differenzen. Diese können dureh Fehlerhäufung erheblich falsch werden. Das schadet 
aber niehts, weil es nicht mehr die Funktionswerte beeinflußt. 

Alle vorangehenden Ziffern der Differenzen und der Funktionswerte bestimmt man dann 
venau nach Verfahren III, also von rechts nach links (d. h. beginnend bei der Differenz 
höchster Ordnung und ankommend beim Funktionswert), nieht wie bei der letzten Ziffer von 
links nach rechts. Nur braucht man jetzt keine Zusatzstellen mehr. Nachdem die letzte 
Ziffer für jeden Teilpunkt aller groben Intervalle festgelegt ist, kann die Summation auch 
über die Grenzen der groben Intervalle fortlaufen. Im Gegensatz zum Verfahren III braucht 
man nicht mehr für jedes grobe Intervall je eine feine Differenz von jeder Ordnung gesondert 
zu berechnen. 


V. 

Hat man die Comrieschen Tafeln nicht zur Hand, so kann man die Beiträge der groben 
Differenzen zur letzten Ziffer der einzuschaltenden Funktionswerte y, mit dem Rechen- 
schieber auf folgende Weise bestimmen. 

I. Man dividiere II, +11, durch 40 und schreibe den Quotienten in der Form 


+IL| n 


mit m =0|1|2...,; n=0|1j2 39. Mit 3 (wie überhaupt mit deutschen Buchstaben) sei 
stets eine Ziffer O0 1.2.....9 bezeichnet. (3, 3.) stehe abkürzungsweise für 3, 10' +3, 10°. Dann 
ist der dreiziffrige Beitrag Y, der zweiten Differenz zur letzten Stelle von %, 


n 
10 


2 


Die Vorzahlen «a, y entnimmt man für die Teilpunkte 0,10,2....)09 dem 
folgenden Täfelchen: 


a Y 7 
0.1 0,9 0.9 0.225 10,6 
0.2 0.8 6 1.6 0,100 + 08 
0.3 0.7 21 0,525 - 0,7 
0,4 0.6 2.4 0,600 04 
0,5 0,5 5 250,685 00 


Der errechnete Beitrag gilt als positiv, wenn // negativ ist. Bei positivem II empfiehlt es sich, 
den negativen Beitrag durch seine dekadische Ergänzung zu ersetzen (z.B. + 6,45 statt - 3,55). 


2. Man bringe I/I/,,, auf die Form 
A 
Für den Beitrag zur letzten Ziffer von 9, bekommt man 
n 3 


Unterhalb 0,5 gehe man bei negativem III, oberhalb 0,5 bei positivem /11 zur dekadischen 
Ergänzung über. 


3. Man dividiere IV, + IV, durch S0OO und schreibe den Quotienten in der Form 


n 


n=0|1|2|...17%9. 


33 95 
51 
56 
40 
22 
| | | 
= 
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Der dreiziffrige Beitrag Y, der vierten Differenz zur letzten Ziffer von 9, ist dann gegeben 
dureh 


n | 
% u v 

0.1 0.0 () 3.5 1.35 3.14 0,392 

0.3 0.7 3.5 135 71.14 0.067 

0,4 0,6 0 6.0 9,60 1,120 

0.5 0.5 5 5) 3.75 90.38 1.172 


Bei negativem /V gehe man zur dekadischen Ergänzung über. Gewöhnlich wird die Rech- 
nung ziemlich einfach werden, weil von den Ziffern 3 nur wenige nicht —0 sein werden. 
v1. 

Das Comriesche Verfahren soll an unserm früheren Zahlenbeispiel vorgeführt werden. 
Wir haben nach Sstelliger Reehnung die Zwischenwerte nur 6stellig erhalten. Wir gehen 
daher jetzt sofort von 6stelligen Werten aus (Zahlentafel 3). Zunächst sind aus den groben 


Zahlentafel>. 


| + 
0.817 513 14 884 1386 249 
65 388 1 627 —- 704 434 
0,0 0,883 201 114266 16 510 I 165 2.089 + 731 677 597 271 
IS 2 58 | 462 4206 + 27 1517 169 
1,0 0.032 079 sı 707 16 048 +3 041 2 117 — 786 
32 829 2 579 - 813 
1.1 0,964 908 15 469 1 304 


Differenzen, die in der Zahlentafel 3 in der Zeile 0,95 stehen, die Endziffern der feinen Diffe- 


renzen zu bestimmen. 
Zahlentafel 4. 


0.1 0.2 03 0.4 0.5 0.6 0,7 0,8 0,0 

Yy 10 10 10 10 10 10 10 10 10 10 
+78 TS 3 12 46 24 
>4 33 49 5%) 93 24 26 
Ill 459 35 37 32 15 0 82 68 63 72 
Il 4500 17 25 36 0 20 05 36 25 7 
Summe 19 52 SU 60 00 53 46 S7 


Von 7° 488785 brauchen nur die beiden letzten Ziffern 78 beibehalten zu werden (Zahlen- 


tafel 4, zweite Spalte). Von /1== 32555 werden Vielfache von 4000 abgestoßen, so daß 
558 übriebleibt. Ferner bekommt man 
II 162 IV= 12% 
| 
. +) 
459 1509 


/ahlentafel 4 enthält die zweiziffrigen Beiträge der vier Differenzen zu den Endziffern der ge- 
suchten Zwischenwerte, angereben in Zehnteln (d. h. hier in Einheiten der 7. Dezimale). 
Diese Zilfernpaare schreibt man entweder aus den Comrieschen Tafeln ab oder bestimmt sie 
nach den Anweisungen unsers Abschnitts V. Die abgerundeten Spaltensummen sind die End- 
zlfern der gesuchten 9 Zwischenwerte. Sie und ihre Differenzen stehen fettgedruckt in Zahlen- 
tafel 5. Ob man in einer Spalte abwärts oder aufwärts zu subtrahieren hat, erkennt man 


| 
Ir: 
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geben aus dem Gang der groben Differenzen Zahlentafel 5. 
in Zahlentafel 3. Die um Null schwan- 
kenden 3. Differenzen sind als ver- y=f(r) 7 
schwindend zu erachten. 
Nachdem die Endziffern feststehen, 7 
ist noch je eine feine Differenz aus jeder 0,30 0,853 201 ia. 167 
Spalte vollständig zu berechnen nach ’ 0.888 849 168 
einem der unter III erwähnten Verfahren. > 473 4 
Welche Werte wir dazu wählen (in > 0.804 315 167 
weleher Zeile von Zahlentafel 5), ıst an > 306 2 
sich gleichgültig. Um später den Fort- 137 
sang der Summation nicht zu unter- 4 0.904 758 168 
brechen, wählen wir am besten die 4.969 -2 
Rech- Werte, die in der Zahlentafel 5 in den 0.309 724 166 
rden. beiden Anfangszeilen 0,900 und 0,905 ’ 4 808 
6 0.014 >30 168 
stehen (kursiv gedruckt). Diese Werte 4.635 +3 
rden bestimmen wir am besten nach Stirling 7 0,919 165 165 
aus den groben Differenzen, die in 1 470 —1 
rn Zahlentafel 3 in der Anfangszeile 0,9 des 923 635 4 304 166 2 
Intervalls stehen, nämlich nach folgen- $ 0.927 939 164 u 
den Formeln: 4 140 
1,00 0.932 079 
111 
IV, 
0,00 
II on 0.825 I 0,155 
1000 
I ? 0.825 (Ill. 0,5 0,24 0,5 + 
Jiffe- 20 2 100 
und erhalten 


Jetzt kann die Zahlentafel 5 durch Aufsummieren von rechts nach links vervollständigt 
werden, ebenso wie früher beim Verfahren II. 

Das Comriesche Verfahren der Untertafelung hat so große Vorzüge, daß es vermutlich 
die ältern Verfahren verdrängen wird. 450 
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Über die Kräfte auf einen elliptischen Zylinder, der sich in 
einer idealen Flüssigkeit bewegt. 


Von A. Ratib in Istanbul. 


ei Betätigung eines Ruders (an einem Flugzeug) treten Momente auf, die zuweilen erheb- 
liche Werte annehmen können. Wir wollen dieses Problem im folgenden untersuchen, 
wobei wir es in geeieneter Weise idealisieren. 

Wir betrachten die ebene Bewegung einer inkompressiblen, idealen Flüssigkeit, die im 
Unendlichen ruht, und in der sich ein elliptischer Zylinder bewegt. Die Bewegung der Flüssig- 
keit sei wirbelfrei, doch soll eine Zirkulation um den Zylinder nicht ausgeschlossen sein. Der 
elliptische (Querschnitt des betrachteten Zylinders drehe sich um einen Punkt O seiner großen 
Achse nach irgendeinem gegebenen Gesetz; diese Drehung braucht nicht gleichförmig zu sein. 
Der Drehpunkt 0 selbst bewege sich mit einer nach Größe und Richtung konstanten Ge- 
schwindigkeit 7/,. In der Abbildung sei 0. die gegebene Richtung dieser Geschwindigkeit. 
/um Koordinatenanfang nehmen wir den Mittelpunkt 0, der kllipse, - und y-Richtung sollen 
in die große und kleine Achse fallen. © sei der Winkel zwischen der großen Achse und der 
festen Riehtung O  .,. Sind g und y Potential- und Stromfunktion für die Bewegung der 


0 


Flüssigkeit, so sind die und „-Komponenten der Geschwindigkeit eines Flüssigkeitsteilchens: 


Im Unendlichen soll die Flüssigkeit 
ruhen, dort muß also =0 sein. 
Für einen Randpunkt M der Ellipse 
erhält man als Ausdruck der Tatsache, 
daß die Normalkomponente der Relativ- 
geschwindigkeit der Flüssigkeit ver- 
schwinden muß, die Randbedingung 


v=Uy—-Vx teonst, 


Ist, 


At 
Wir führen jetzt elliptische Koordinaten Z = +7 ein, die mit « und y durch die 
Gleichungen: 
bzw. y=eSinfsiny 


verknüpft sind: dabei ist 2e die Brennpunktsentfernung der Ellipse. Die Kurven £= const 
sind Ellipsen, die Kurven »,  const konfokale Hyperbeln. Der Umriß des elliptischen Zylinders 
sei die Ellipse 

Wenn man «die klassische Lösung «des einfachen Falles, wo sich die Ellipse um ihren 
Mittelpunkt O, dreht, etwas abändert, so bekommt man das komplexe Strömungspotential für 
unseren allgemeineren Fall in der Form 


0 s) 


Hier sind und be Zins, die Halbachsen der Ellipse, ist die Konstante der 
Zirkulation um den Zylinder und 4°00, ist die Entfernung des Drehpunktes 0 von dem 
Mittelpunkte ®, der Ellipse. 

Um zu zeigen, daß der obige Ansatz die Bedingungen im Unendlichen und am Rande 
erfüllt, zerlegen wir ihn in reellen und imaginären Teil: 


B (2). 
r r -( 2 (<u 5) | so 
y' (Ubsın Var cos ı)) eosZye Aamecosye 


Man sieht ohne weiteres, daß die Bedingung e = r=0 im Unendlichen erfüllt ist. Am Rande 
der Ellipse soll nach Gl. Uy- R?/2-+ const sein. Nun ist für einen Rand- 
punkt M der Ellipse die Entfernung R von dem Drehpunkt O0 gegeben durch 
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was mit ec &, cos = acosy, bsinnnach einfacher Umformung 


ergibt. Damit wird der am Rande vorgeschriebene Wert von y 
v=Uy—Vx teonst= UÜbsiny- Vacosı 


Das stimmt mit dem aus (2) folgenden Werte (für £ 5,) 


Ubsiny- Vacos eos 2n+5 = COS 


überein, wenn die in (1) freibleibende Konstante den Wert 


0 
- 2 

4 


Irr 


erhält. Die Randbedingung an der Ellipse ist also erfüllt. 


Um die Kräfte zu berechnen, die von der Flüssigkeit auf das Hindernis ausgeübt werden, 
insbesondere das Moment, kann man entweder die klassische Methode zur Bestimmung der 
Kräfte auf einen festen Körper in einer wirbelfrei bewegten, idealen Flüssigkeit anwenden, 
oder direkt die Formel von Bernoulli benutzen. 


Erste Methode: Wir nehmen zunächst an, die Zirkulation sei Null. Ist o die Dichte der 
Flüssigkeit, 7 die kinetische Energie der Strömung pro Längeneinheit des Zylinders, so ist 


= o\ \ \ ds, 


vw 


wo n die in die Flüssigkeit weisende Normale am Rande ist. Ist: ds das Längenelement des 


dp dy 


elliptischen Umrisses, so wird also 2T=oiypdy. 


Nun kann man Potential- und Stromfunktion in der folgenden Formel schreiben: 
Uy, Vy,toy,;,. 
Man erhält die Funktionen q; und ; ohne weiteres aus den Formeln (2), indem man nach 
U, V und © ordnet. Eine einfache Zwischenreehnung führt zu dem Ergebnis 


(siehe Lamb, $ 121, p. 151). 

Das Moment, welches pro Längenheit des Zylinders von der Flüssigkeit ausgeübt wird, ıst 
d/oT ‚oT ‚oT 


OU 


Die Durchführung der Rechnung ergibt 


.) 


To 


Durch analoge Rechnung erhält man die Kraftkomponenten 


X 
:— TU, sin® +a 
do 


4 
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Zweite Methode. Formel von Bernoulli: Für bewegte Bezugsachsen und für unseren 
speziellen Fall nimmt die Bernoullische Druckgleichung die folgende Form an (siehe 
20, p. 18: 


Wir setzen dabei voraus, daß keine äußeren Kräfte wirken; F'(f) ist eine willkürliche Funktion 
der Zeit. Ferner ist 


(uw — U’+(e— V). 


Dann wird das Moment des Flüssigkeitsdruckes pro Längeneinheit des Zylinders durch das 
um die Ellipse herum zu nehmende Integral 


ds 


geliefert, wo p, und p, die Projektionen des Druckes auf die «- und y-Richtung sind. Dies 
Integral wird in unserem Spezialfalle 


d 


In (Rk)dn, 


N= 


wo mit Hilfe von Formel (4) 


Po = P (So, N) 
und mit Formel (3) 
B,=Bi},n) 


einzusetzen ist. Die Auswertung des Integrals ergibt 


N eT’sin209 12.0? | cos 
Wenn die Ellipse sich ohne Drehung parallel verschiebt und eine Zirkulation nicht auftritt, 
so ergibt sich 


Wenn die Rotation gleichförmig ist, verschwindet der mittlere Ausdruck in Formel (5). 


Man kann die Zirkulation z. B. dureh die Korderung bestimmen, daß am Ende der 
großen Achse der Ellipse die Geschwindigkeit zu null wird. Dadurch bestimmt sich 7" zu: 


# 
\ 


2/dt' 

Man sieht, daß man mit dieser oder einer gleichwertigen Hypothese eine Momentenformel 
bekommt, die sieh zusammensetzt: 1. aus Ausdrücken, die von der Rotation nieht abhängen : 
2. aus Ausdrücken, die in der Winkelgesehwindigkeit linear sind; 3. aus Ausdrücken, die in 
der Winkelbeschleunigung linear sind. 


Um den Spezialfall der ebenen Platte zu erhalten, genügt es, in der Formel 5) e=a 
zu setzen. Die Geschwindigkeit wird dann an beiden Kanten der Platte unendlich. Man 
kann über die Zirkulation /'" so verfügen, daß die Geschwindigkeit an einer der Kanten end- 
lieh bleibt. Soll sie (siehe Abb.) bei ‚I oder bei .1’ endlich bleiben, so muß 


sein. 

Es gibt einen Sonderfall, wo die Geschwindigkeit an beiden Kanten endlich bleibt. 
Dafür ist erforderlich, daß 


ige" 


ist, wo ©, der Wert von © für f=-#, ist. 449 
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Über eine spezielle Markoffsche Kette am Galtonbrett. 
Von H. Münzner. 
(Aus dem Institut für mathematische Statistik der Universität Göttingen.) 


as Galtonbrett wird in der Wahrscheinliehkeitsreehnung zur Veranschauliehung der 
Binomialverteilung verwendet. Es handelt sich dabei um ein schräg aufgestelltes Brett, 
auf dem ein System von Stiften in horizontalen Reihen derart angeordnet ist, daß ein Netz 
von gleichschenkligen Dreiecken gebildet wird. In dieses Gitter von Stiften läßt man oben 
an einer bestimmten Stelle kleine Kugeln aus Stahl, Holz oder einem anderen Material laufen, 
die unter dem Einfluß der Schwere nach unten rollen und sieh in den auf der unteren Seite 
des Stiftefeldes angebrachten Fächern sammeln. Die Kugeln müssen so groß sein, daß sie 
in jeder Stiftereihe auf einen der Stifte stoßen und sich dort entscheiden, ob sie in schräger 
Riehtung rechts oder links von diesem Stift weiterrollen. 
Ist die Wahrscheinlichkeit, nach rechts oder links zu fallen, für jede Kugel an jedem 
Stift gleiel ki ‚so werden, wenn s Reihen von Stiften vorhanden sind, in den s-+ 1 in Frage 


kommenden Auffangfächern, die wir der Reihe nach mit @=-0,1,2...s bezeichnen wollen, 
gerade | 

Ns! 1 \° 

(1) 

Kugeln zu erwarten sein. N bedeutet dabei die Gesamtzahl der benutzten Kugeln, und ga, s) 
stellt die bekannte Binomialverteilung dar, die etwa von s= 12 angefangen in recht guter 
Annäherung durch die Gauß verteilung h 


(2) 
ersetzt werden kann. 

Bei zahlreichen praktischen Versuchen am Galtonbrett, die im Institut für mathematische 
Statistik zur Durchführung von Funktionaltransformationen') ausgeführt wurden, hieß sich 
eine derartige Gaußkurve nicht genau erzielen, sondern es traten immer wesentlich Nlachere 
Kurven?) auf. Das gab zur Vermutung Anlaß, daß das Galtonsche Prinzip der gleichen 
Wahrscheinlichkeiten, die für das Rechts- oder Links-Hinunterrollen der Kugeln gelten sollten, 
in Wirklichkeit nicht erfüllt war. 

Das Studium vieler Kugelbahnen, die von H. Bitterling mit Hilfe von berußtem 
Papier festgehalten wurden, ergab dann auch, daß die Kugeln eine größere Tendenz haben, 
die einmal eingeschlagene schräge Links- oder Reehtsriehtung beizubehalten, als diese Richtung 
zu wechseln, was sich leicht aus dem Beharrungsvermögen der Kugeln erklären läßt. Je 
lacher das Galtonbrett, das drehbar befestigt war, gestellt wurde, um so stärker war diese 
Tendenz des Beibehaltens der einmal eingeschlagenen Richtung. 

Damit war offenbar das Wahrscheinlichkeitsschema, das zur regulären Bernoulli- 
verteilung führt, gestört, und es ergab sich die Notwendigkeit, der Berechnung der Verteilungs- 
funktion (a, s), die sich unten in den Fächern bildet, ein anderes Wahrscheinlichkeitsschema 
zugrunde zu legen, bei dem der genannte Tatbestand berücksichtigt wurde. 

P und © seien die Wahrscheinlichkeiten für das Links- und Rechts-Hinunterfallen am 
ersten Stift, p sei die Wahrscheinlichkeit, daß die einmal eingeschlagene Richtung bestehen 
bleibt, und 9=1-- p die Wahrscheinlichkeit dafür, daß ein Richtungswechsel eintritt. Unter 
einmal eingeschlagener Richtung soll dabei nur die unmittelbar vorausgehende Richtung ver- 
standen werden, von der allein die folgende Entscheidung abhängen soll, während wir keinen 
Einfluß von den früheren Richtungen her annehmen wollen. Auch soll das Springen der 
Kugeln über den linken oder rechten Nachbarstift hinweg vernachlässigt werden, da es bei 
den vorliegenden Versuchen verhältnismäßig selten eintrat. 

Bezeichnen wir das Links-Hinunterfallen einer Kugel mit dem Ereignis 0, das Rechts- 
Hinunterfallen mit dem Ereignis 1, so sind P und 9 die Anfangswahrscheinliehkeiten für die 
Ereignisse O0 bzw. 1, während die Wahrscheinlichkeit für 0 unter der Voraussetzung, daß 
unmittelbar vorher ebenfalls das Ereignis 0 eingetreten war, gleich „p, = p Ist und die Wahr- 
scheinlichkeit für 0, nachdem vorher das Ereignis 1 stattgefunden hat, gleich ‚p, Ist. 
Entsprechend sind „p,= und ‚p, = p definiert. 


I) F. Bernstein: Verallgemeinertes Galtonbrett zur Durehführung von Funktionaltransformationen, Ztsehr. 
f. Phys., Bd. 77, 1932. 

2) Über die gleiche Erfahrungen beriehteten K. Odenhausen: (Ztsehr. f. math. u. naturwiss. Unterricht 1951) 
und W. Seitz mit K. Hamacher-Odenhausen (Physikalische Zeitschrift 35. Jahrg. 1934). 
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Wird die Summe sämtlicher Ereignisse bei s Versuchen mit a bezeichnet, so stellt die 
Wahrsecheinlichkeit (a, s), daß bei s Versuchen gerade a eintritt, die gesuchte Häufigkeits- 
verteilung dar, die bei s Reihen in den s— 1 Auffangfächern des Galtonbrettes zu erwarten 
ist, sofern der obige Ansatz den wirklichen Vorgängen entspricht. 

Um die Wahrscheinlichkeit y(as) zu berechnen, zerlegen wir sie in die beiden 
Summanden ,(a,s) und ,(as). Dabei sind und y,(a,s) die Wahrscheinlichkeiten 
dafür, daß eine Kugel bei s Reihen von links her bzw. von rechts her in das a-te Fach gerollt 
kommt. Für diese Wahrscheinlichkeiten gelten offenbar die Differenzengleichungen 

s)| 


(3), 
v,(a,s+1) pyp,(a,s) + s)| 


aus denen dann nach der bekannten Methode sich die erzeugenden Funktionen F,(r, y) und 
y) herleiten lassen, deren Summe 


gerade die gesuchte Wahrscheinliehkeitsverteilung y'(a, s) erzeugt. 


1 
Für lautet diese erzeugende Funktion 


1 5) 
1—-zyp) i—yp— 


Fi, y) 


woraus sich durch Reihenentwieklung nach = und y die Wahrscheinlichkeit y'(«a, s) sofort 
angeben läßt. Nach einfachen Umformungen ergibt sich dann: 


7 | l | ‚| für <a<s 
| 


Tabelle I enthält für » 0,1: 03:05: 0.7: 08:09: 095, 0,99 und 0,999 bei s— 10 die 
verschiedenen nach Formel (6) berechneten Wahrscheinlichkeitsverteilungen. Dabei sind die 
Wahrscheinliehkeiten mit 10000 multipliziert. 


Tabelle 1°). 


1 2 4 5 
p 10 fe) 7 6 
0 1 63 1680 6512 
0.3 0 662 2391 3718 
0,5 10 US 410 1172 2050 2460 
0.7 202 469 1207 1490 1594 
671 875 1040 1151 1190 
1958 621 666 715 124 
0.05 311 402 409 413 416 418 
0,00 4568 36 36 9 
0000 4959 10 | 10 10 10 10 


Für kleine p haben wir steile gloekenförmige Verteilungen, die mit wachsendem » immer 
Nacher werden und schließlich in #-förmige Verteilungen mit einem Hügel in der Mitte über- 
gehen, der immer niedriger wird und im Grenzfall p=1, wo wir den 2-Kategorientyp haben, 
ganz verschwindet. Bei p 0,999 läßt sich dieser Hügel auch noch nachweisen, wenn wir 
die Berechnung der Ordinaten auf neun Dezimalstellen vornehmen. Für «= 1,2...9 stimmen 
die ersten sieben Dezimalstellen alle überein und lauten 0,0009960, während für die achte 
und neunte Dezimalstelle sich ergab: 


= 2/1 4|5|/6| 9 
12 15 16.15. 12.07.00 
Krst im Grenzfall p 1 verschwindet dieser Hügel vollkommen. 


*) Diese Tabelle stellt nur einen kleinen Auszug der umfangreichen numerischen Bereehnungen dar, die von 
NH. Patzer durehgeführt wurden. 
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Bei verschiedenen Anfangswahrscheinlichkeiten P und @ ergeben sich entsprechende 
schiefe Verteilungen. 

Bei wachsendem s tritt der ««-förmige Verteilungstyp immer mehr zurück, und wir er- 
halten schon für verhältnismäßig große p glockenförmige Verteilungen. 

Von besonderem Interesse ist der Grenzfall s>®». Wir können hier das Resultat 
sofort angeben, wenn wir beachten, daß das vorliegende Wahrscheinlichkeitsschema eine 
Markoffsche*) Kette mit der Grunddeterminante: 


darstellt. Lassen wir zunächst die Sonderfälle p=1; 0,5; und O0 außer acht, so können wir 
ohne weiteres das allgemeine Markoffsche Resultat für derartige Wahrscheinlichkeits- 
verkettungen auf unsere spezielle Wahrscheinlichkeitskette übertragen. Markoff beweist, 
daß mit wachsendem s die Wahrscheinlichkeitsverteilung unabhängig von den Anfangs- 
bedingungen gegen eine Gauß verteilung’) geht. Demnach gilt 


ur la, 
so z\a,s) 
wobei 
| sa)? 
Ss): 
my2sa 


ist und sich a und m aus dem Nenner @ der erzeugenden Funktion y) (4) 
folgender Weise ergeben: 


(9). 


Für unseren Spezialfall erhalten wir nach Formel (4): 
und daraus: 
l p 
(10) 


Den Formeln (S) und (10) entnehmen wir nun das Ergebnis, das wir bei großen s angenähert 
für y(a,s) die Gauß verteilung 
s\2 
2q(« 
la, s): / E 
‚la,: 
setzen können. 
| 
Der Fall p 


ist nunmehr auch in dieser Formel enthalten und ergibt die gewöhn- 


liche Laplacesche Formel. Ist so erhalten wir flachere Gaußkurven, Ist p< ,, 0 


entstehen steilere Gaußkurven. Bei p=0 haben wir den Ein-Kategorientyp, bei p° 1 
den Zwei-Kategorientyp, während bei p=1- e, wobei « eine auch noch so kleine feste Zahl 
bedeuten kann, schon Gaußverteilung vorhanden ist. 
Man kann leicht zeigen, daß Formel (11) mit Ausnahme der Fälle p 0 und p=1 auch 
dann gültig ist, wenn die Anfangswahrscheinlichkeiten P und 0 voneinander verschieden sind. 
Das in den Formeln (6) bzw. (11) enthaltene theoretische Resultat wurde nun mit den 
wirklichen Vorgängen am Galtonbrett in folgender Weise verglichen. 


An einem Galtonbrett. dessen Stifte in 18 Horizontalreihen in Form von gleichschenkligen Drei- 
ecken von 1 em Basis und 1 em Höhe angeordnet waren, ließen wir von oben an einer bestimmten 
Stelle 1000 Stahlkugeln von 2,381 mm Durchmesser einlaufen. Die Neigung des Brettes betrug 20 Grad, 
sie wurde so gering gewählt, weil dadurch die Kugeln mehr die Tendenz hatten, in der einmal ein- 
zeschlagenen Riehtung zu verharren und so ein größeres p erzielt werden konnte. Ferner war diese 
geringe Neigung nötig, um das von der Elastizität der Kugeln herrührende Springen über 1 oder 
2 Stifte seitwärts hinweg auf ein Minimum zu reduzieren. Ganz vermeiden ließ sich jedoch dieses 
Springen nicht, so daß auch in die äußersten Randfächer, die nach der Theorie keine Kugeln mehr 
enthalten dürfen, einige Kugeln hineinfielen. Das Ergebnis des Versuchs, d. h. die Kugelanzahl in 
jedem einzelnen Fache ist in Spalte I der Tabelle 2 dargestellt. Spalte II enthält die Erwartungswerte 


4, A.A. Markoff-H. Liebmann: Wahrscheinlichkeitsreehnung, Leipzig 1912, S. 272 ff. 
5) Von K. Mahler stammt für das Vorliegen von Gaußverteilung für den Grenzfall dieses Wahrscheinlichkeits- 
schemas ein direkter Beweis, der wegen seines Umfangs nieht veröffentlicht werden konnte. 


; 
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und mittleren Fehler für diese Anzahlen bei Zugrundelegen der gewöhnlichen Bernoulliverteilung. 
Diese Werte weichen sehr stark von den empirischen ab. In Spalte III sind die Erwartungswerte nach 
der vorstehenden Theorie und ihre mittleren Fehler angegeben. Zur Berechnung dieser Erwartungs- 
werte diente die Formel (6), wobei für die Wahrscheinliehkeit » die Größe 0,75 verwendet wurde. 
Dieser Wert 0.75 wurde empirisch durch Auszählung der Entscheidungen der Kugeln. deren Bahnen 
wir auf berußtem Papier festhielten, gewonnen. Bei einer Gesamtzahl von 294 Entscheidungen ergab sieh 


219 
>19mal ein Gleichbleiben der Richtung und 75mal ein Riehtungswechsel, woraus wir P= 59, 775 


erhielten. 

Die Übereinstimmung zwischen Spalte III und Spalte I ist sehr gut, fast immer liegen die empi- 
rischen Werte im Bereich des einfachen mittleren Fehlers um den theoretischen Wert. Spalte IV enthält 
zum Vergleich die nach der asymptotischen Formel berechneten Erwartungswerte für die einzelnen 
Ordinaten, die nicht mehr wesentlich von dem genauen Wert in Spalte III abweichen. 


Tabelle 2. 


empirische Bernoulli- Markoffsche Kette Markoffsche Kette 
Werte verteilung (genaue Werte)  (asympt. Formel) 

10 | 4 — | 2,1 
5 0,0 38 +19 | 5.4 

13 0,1 22 +30 | 10.2 

7 | > 0.6 182 +42 16.5 

6 35 31+ 1,8 31.0 +54 28,6 
-5 50 34 45.8 + 6,6 43,0 
| 68 328 + 56 60,0 

3 | 77 70.8+ 81 79.6 + 8.6 | 77.8 

> Si] 121,4 + 10,3 93,5 + 9,3 | 33.6 

| % 166.8 + 11,8 102,5 + 9.6 104.6 

110 g 18544123 106,0 + 9,7 108.5 

| 166.8 + 11.8 1028 +96 104.6 

> 88 121,4 + 10.3 93,5 + 9,3 93.6 

3 74 708 + 81 79.6 + 8,6 17,8 

4 60 32,8+ 56 63,1+ 77 60.0 
45 11.7 + 3.4 45.8 + 6,6 43,0 

6 38 3.1%+ 1,8 31,0 + 5,4 28,6 

7 25 0,6 18,2 + 4,2 16.5 

S S 0,1 9024380 | 10,2 

5 0,0 3,8 +19 | 5.1 
10 2 — | 2,7 


Die gute Übereinstimmung zwischen der theoretisch und empirisch gefundenen Häufig- 
keitsverteilung zeigt, daß die wirklichen Vorgänge am Galtonbrett durch eine derartige 
Markoffsche Kette hinreichend genau erfaßt werden. Die Berücksichtigung weiterer syste- 
matischer Ursachen für die Abweichung von der gewöhnlichen Bernoulliverteilung hätte 
nur einen Sinn für Versuche mit sehr vielen Kugeln, weil dann die systematischen Ab- 
weichungen gegenüber den Zufallsschwankungen groß würden, während im vorliegenden Fall 
bei 1000 Kugeln die Zufallsschwankungen, die durch die mittleren Fehler charakterisiert sind, 
größer sind als möglicherweise vorhandene weitere systematische Abweichungen. 

Schließlich sei noch darauf hingewiesen, daß das vorliegende Wahrscheinlichkeitsschema 
zu denjenigen Wahrscheinlichkeitsschematen gehört, die zu dem sogenannten statistischen 
Ausgleich führen. Bei Erscheinungen, denen dieses Wahrscheinlichkeitsschema für den 
Fall p<gy zugrunde liegt, ist nämlich ein häufigerer Wechsel zwischen den Alternativereig- 
nissen zu erwarten als im gewöhnlichen Bernoullifall. Damit ist eine größere Häufigkeit 
der kleineren Iterationen und das Auftreten von steileren Häufigkeitsverteilungen verbunden; 
es sind also gerade die für den statistischen Ausgleich typischen Merkmale vorhanden. Damit 
soll nieht behauptet werden, daß die Beispiele von K. Marbe sich durch dieses Schema 
erklären lassen, sondern nur wieder einmal gezeigt werden, daß Marbes Lehre vom sta- 
tistischen Ausgleich nicht die Anwendungsmöglichkeiten der Wahrscheinlichkeitsrechnung, 
sondern lediglich die die Anwendung des Bernoullischemas erschüttert. 445 
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Die Ermittlung von Potentialströmungen durch ein 
experimentelles Verfahren. 
Von W. Barth in Völklingen (Saar). 


ie rechnerische Ermittlung von Potentialströmungen stößt in vielen Fällen auf unüber- 
D windliche mathematische Schwierigkeiten. Im folgenden soll eine experimentelle Methode 
zur Bestimmung der Potentialströmung beschrieben werden, die selbst verwickelte Aufgaben 
in einfacher Weise zu lösen gestattet. 


Für eine Strömung durch geschichtete Stoffe gelten bei kleinen Reynoldsschen 


Op Op Op 


/Jahlen zwischen dem Druckgefälle , und den Komponenten der Geschwindigkeit 
02’ 


€, die Gleichungen!) 


02 
wo k eine die Durchlässigkeit charakterisierende Materialkonstante ist. Aus der Kontinuitäts- 
bedingung 


folgt dann 


dp 


++ 


02° (), 
d.h. die Strömung ist eine Potentialbewegung; die Potentialflächen sind die Flächen gleichen 
Druckes. Diese Erkenntnis kann man dazu verwenden, um durch Messung des Druckes in 
geschichteten Stoffen die Potentialströmung zu bestimmen. Die Durchführung dieses Ver- 
fahrens soll an einem Beispiel gezeigt werden. 


Anläßlich einer Untersuchung über Strömungsvorgänge im Hochofen erschien es 
wünschenswert, den Verlauf einer idealen Strömung im lHochofengestell kennenzulernen. 
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einen Hochofen. Versuchsapparät. 


W. Barth: Der Druckverlust in geschiehteten Stoffen.“ Forschung 1933, S. 83. 
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In Abb. 1 ist ein Schnitt durch einen Hochofen wiedergegeben. Den Unterteil des Ofens 
kann man als Rotationskörper auffassen, in welchen an mehreren Stellen — entsprechend der 
Zahl der Blasformen Luft eingeblasen wird. Dieser Körper kann aus einer Anzahl gleicher 
Segmente zusammengesetzt gedacht werden, wobei aus Symmetriegründen in jedem Segment 
die Strömung ähnlich verlaufen muß. Man kann sieh daher auf die Untersuchung eines 
Seementes beschränken. 

Dieses Segment wurde in Blech ausgeführt und ist in Abb. 2 dargestellt. Der Einfach- 
heit halber wurde die Konizität des Ofens vernachlässigt. Das Gefäß wurde mit kugeligen 
(Juarzkörnern von ca. I mm & (Normsand) gefüllt. Durch eine Düse, welche durch ein grobes 
Sieb abgeschlossen worden war, kann Luft abgesaugt oder eingeblasen werden. Die Luft- 
menge wurde durch eine Blende gemessen. Die Wandungen des Versuchsgefäßes waren an 
zahlreichen Stellen durehbohrt und mit Stutzen zum Anschluß eines mit Wasser gefüllten 
!-Rohres zur Messung des Druckes versehen. Zur Messung des Druckes im Innern wurde 
ein Messingröhrehen von 2 mm Außendurehmesser, welches an ein T-Glas angeschlossen 

worden war, dureh eine Meßstelle in das Schüttgut ein- 


£bene — geführt. Es kann dabei angenommen werden, daß durch 
f —5 das Einführen dieser Röhrehen «die Druckverteilung nicht 
F® or merklich gestört wird. Der Druck wurde an allen Stellen 


— der Behälterwand und an ausgezeichneten Punkten im 
Innern des Materials ermittelt. Mit Hilfe dieser Werte 
können die Flächen gleichen Druckes dureh Interpolation 
leicht bestimmt werden. Die Ergebnisse sind in Abb. 3 
für die wichtigsten Schnittebenen zusammengestellt. 

Aus dem Abstand der Druckflächen kann auf die 
Größe der Geschwindigkeit geschlossen werden. Zu diesem 
Zwecke muß der Druckabfall pro Längeneinheit für ver- 
schiedene Gesehwindigkeiten bekannt sein. Diese Werte 
wurden experimentell ermittelt und sind in Abb. 4 zu- 
sammengestellt. Die Riehtung der Geschwindigkeit kann 
mit Hilfe der Tatsache, daß die Stromlinien die Druck- 
flächen durehschneiden müssen, bestimmt werden. 
Da die Stromlinien selbst nur schwer darzustellen sind, 
wurden in Abb. 3 nur die Linien eingezeichnet, die dureh 
Verbinden der Projektionen der Geschwindigkeiten in den 
einzelnen Schnittebenen entstehen. 

Die ermittelten Flächen entsprechen nur angenähert 
den Potentialllächen. Dies rührt daher, daß für den ver- 
wendeten Sand Gl]. (1), wie aus Abb. 4 hervorgeht, in dem 
fraglichen Bereich nur angenähert gilt. Für die dureh- 
geführte Untersuchung war diese Abweichung ohne 
Abb.3. Strömungsverlauf und Flächen gleichen Bedeutung, da für den vorliegenden Zweck nicht die 

Druckes (Potentialflächen,) Potentialströmung, sondern vor allem die Strömung beı 
größeren Reynoldsschen Kennziffern von Interesse war. 
In Fällen, in denen die Potentialflächen sehr genau fest- 


= gestellt werden müssen, empfiehlt es sich, gegebenenfalls 
ein Füllmaterial zu wählen, dessen Korngröße kleiner ist 
I (vgl.')). Beim Einfüllen des Sandes ist darauf zu achten, 
* dafs die Lagerung an allen Stellen gleich dieht ist, so daß 
SB der Widerstand des Füllmaterials an allen Stellen und 
Q 
N nach allen ltichtungen als ‚gleich angenommen werden 
sw kann. Eine gleichmäßige Verdiehtung kann evtl. dureh 
& Klopfen und Drehen erreicht werden. Randwertprobleme 


zweiter Art können dadurch gelöst werden, daß man an 
bestimmten Stellen Luft eintreten bzw. absaugen läßt, wie 
| | RE es der geforderten Geschwindigkeitsverteilung entspricht. 
4 Die beschriebene Methode läßt sich für die verschieden- 
Geschwindigkeit c der Luft Prohl dürfte in fast alle 
im freigedachten Querschnitt artigsten robleme verwenden und dürfte in fast allen 
Fällen die Bestimmung der Potentialströmung ganz er- 

Abb, 4. Druckabfall pro 100 mm Sehütthöhe 

in Quarzsand von ea. I mm Korngröße. heblieh erleichtern. 


in 


Zusammenfassung. Durch Messung der Druckverteilung in geschichteten Stoffen bei der 
Durehströmung von Luft wird die Potentialströmung ermittelt. 446 
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Vorträge der Hauptversammlung in Bad Pyrmont 
der Gesellschaft für angewandte Mathematik und Mechanik. 


Von den auf der Versammlung in Bad Pyrmont vom 9. bis 14. September gehaltenen Vorträgen 

bringen wir im folgenden die von den einzelnen Vortragenden selbst zur Verfügung en 

Vortragsauszüge. Die vollständige Veröffentlichung einzelner Vorträge bleibt vorbehalten. 


Topologische Dualitätssätze 
und Reziprozitätstheoreme der Schaltungstheorie. 
Von W. Canuer in Göttingen. 

n der Elektrotechnik, besonders in der Nachriehtentechnik, hat die Konstruktion sogenannter 
| ee Schaltungen praktische Bedeutung '). Hier soll der mathematische Kern der 
physikalischen Fragestellung herausgeschält werden, wobei sich herausstellen wird, daß eine 
enge Beziehung zwischen den topologischen Dualitätssätzen und unserem Problem, aber nicht 
einfach ein Spezialfall dieser Sätze vorliegt, Mittelbar geben die folgenden Überlegungen 
auch Aufschluß darüber, für welche Fachwerke eın Uremonascher Kräfteplan möglich ist e. 
Unter Beschränkung auf den interessanteren elektrischen Fall werde folgendes Reziprozitäts- 
theorem zum Ausgangspunkt genommen: Mit jeder Funktion fi) io, Kreisfrequenz; 
für die Schaltung, Abb. 1) ist auch (z. B. Abb. 2) als komplexer 
Weechselstromwiderstand durch positive Schaltelemente ohne Elektronenröhren realisier- 
har. Das ist in dem Satz von ©. Brune über die Natur der realisierbaren l unktionen [(4) 
enthalten ’). Das genannte Reziprozitätstheorem gestattet noch Verallgemeinerungen ; weitere 


Abb. 1. Abb. 2. Abb. 2. 


Reziprozitätstheoreme der Schaltungstheorie haben jedoch keinen topologischen Inhalt. Wohl. 
bekannt ist, wie man durch Vertauschung von „Reihe“ und „parallel“ und Ersetzen der in den 
Stromzweigen liegenden Elemente durch reziproke Elemente (C = u.BRB=R'E=6G, 2. B. Abb. l 
und 2) zu einer Schaltung eine Zweig für Zweig widerstandsreziproke („duale*) erhält. Diese 
Konstruktion ist ein Spezialfall der folgenden am Beispiel der Wh eatstoneschen Brücke 
'Tetraederkonfiguration, Abb. 3) erläuterten dualen Konstruktion. Man breite das elektrische 
Netzwerk ohne Überschneidung in die Ebene (auf die Kugelfläche) aus, wähle innerhalb jeder 
Masche einen Knotenpunkt der dualen Schaltung, verbinde diese Knotenpunkte so, daß gerade 
jede Verbindung einen Zweig der Ausgangsschaltung überschneidet und ordne je zwei sich 
überschneidende Zweige (Strecken) einander zu, derart, daß das Produkt der W iderstands- 
werte zugeordneter Zweige 1 ergibt (z. B. entsprieht in Abb. 3 dem Lin Zweig 1 das 
"—=L in Zweig 1’. Daß diese Art dualer Konstruktion tatsächlich stets zu reziproken 
Wechselstromwiderständen führt, leuchtet physikalisch sofort ein, wenn man bedenkt, daß 
die Kirehhoffschen Knotenpunktsgleichungen der Ströme für die eine Schaltung ebenso 
lauten wie die Spannungsgleichungen der Stromkreise für die duale Schaltung und umge- 
kehrt. Nun gibt es aber, was von Elektrotechnikern oft übersehen wird, Schaltungen (Strecken- 
komplexe), für welche diese duale Konstruktion versagt, z. B. der Streekenkomplex Abh. 4. 
Von der weiteren Komplikation, die durch das Auftreten gegenseitiger Induktionen entsteht, 
soll hier der Einfachheit halber abgesehen werden, obwohl der erwähnte grundlegende, von 
Brune bewiesene Satz nur richtig ist, wenn gegenseitige Induktionen zugelassen werden. 
Es erhebt sich die Frage: Welches sind die allgemeinsten Schaltungen, zu denen man auch 
der Struktur nach widerstandsreziproke (duale) konstruieren kann? Diese Frage ist äquivalent 
mit der folgenden: Für welche Schaltungen stimmen die Kirehhoffschen Knotenpunkts- 

E7 Ss. z.B. H. Bitter: Telegraphen- und Fernsprechtechnik 2, 1934, Nr. 3, S. 59 bis 65. | EN 

?, Wenn man den Cremonaschen Satz hinzunimmt, daß für jedes Fachwerk (einschl, Seilpolygon), welches sich als 


senkrechte Projektion eines Eule: schen Polyeders darstellen läßt, tatsächlich ein COremonascher Kräfteplan existiert. 
3) ‚Positive Funktionen“, s. z. B. diese Zeitschrift X, 1031, S. 436 bis 437. 
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gleichungen der Zweigströme mit den passend geschriebenen Kirchhoffschen Kreis- 
gleichungen der Zweigspannungen einer dualen Schaltung überein? Oder auch, wenn die 
nichtwesentliche Einschränkung auf zusammenhängende nichtzerfallende Streckenkomplexe 
gemacht wird: Welches ist die Klasse der Streckenkomplexe, die duale in dem Sinne haben 
(Definition der Dualität), daß die Inzidenzmatrix des Streckenkomplexes mit der Matrix einer 
passenden (nieht notwendig kleinsten) Basis der Zyklen des dualen Streekenkomplexes über- 
einstimmt. Denkt man sich in diese Zyklen Polygone eingespannt, so entsteht ein Flächen- 
komplex, von dem gefordert wird, daß seine zweite Inzidenzmatrix die gespiegelte der (ersten) 
Inzidenzmatrix des dualen Streckenkomplexes ist. Durch Berechnung der Eulerschen 
Charakteristik oder daraus, daß die Basiszyklen, also alle geschlossenen Ketten beranden 
(nullhomolog sind), folgt, daß der Flächenkomplex homöomorph einer Kugel (Eulersches 
Polyeder) ist. Von anderen möglichen Definitionen der Dualität von Streckenkomplexen ’), 
ist folgende besonders bemerkenswert, da sie an die Kirchhoffsche kombinatorische 
Berechnungsweise zusammengesetzter Widerstände anknüpft und somit ebenfalls dem physi- 
kalischen Sachverhalt gerecht wird: Ein Streekenkomplex heiße dual zu einem anderen, wenn 
eine eindeutige Zuordnung zwischen den Strecken beider besteht, so daß jedem vollständigen 
Baum des ersten Komplexes ein System von unabhängigen Strecken des dualen Komplexes ent- 
spricht. Wie nach der ersten, so ist auch nach dieser Definition die Dualität eine wechselseitige. 
Vollständiger Baum eines gegebenen Streckenkomplexes (z. B. 4,5, 6 in Abb. 3) heißt dabei ein 
System von Strecken, das keinen geschlossenen Kreis aufweist, das aber bei Hinzunahme irgend- 
einer weiteren Strecke einen geschlossenen Kreis erhält. Die Komplementärmenge eines 
Systems von Strecken, die einen vollständigen Baum bilden, heißt System unabhängiger 
Strecken (z. B. 1, 2,3 in Abb. 3). Die Spannungen in den Zweigen (Strecken) eines voll- 
ständigen Baumes bestimmen, abgesehen von einer willkürlichen additiven Konstanten, sämt- 
liche Spannungen einer Schaltung, ebenso die Ströme eines Systems unabhängiger Zweige 
sämtliche Ströme. Durch vorwiegend kombinatorische Überlegungen in Verbindung mit dem 
Jordanschen Kurvensatz folgt, daß ein Strecekenkomplex, der nach der letzten Definition 


A439 


Abb. 4. Abb. >». 
einen «dualen Streckenkomplex gestattet, sich ohne Überschneidung auf die Kugel ausbreiten 
läßt und umgekehrt, in Übereinstimmung mit der ersten Definition der Dualität. Ein praktisches 
Kriterium, wann eine duale Konstruktion möglich ist, gibt ein Satz Pontrjagin und 
Kuratowski, nach dem es notwendig und hinreichend für die Ausbreitbarkeit eines Strecken- 
komplexes auf die Kugel ist, daß weder die Konfiguration Abb. 4 noch die Konfiguration 
Abb. 5 als Teilkomplex (mit Zulassung von Zwischen-Knotenpunkten) auftritt’). Für sämtliche 
Schaltungen lassen sich Zweig für Zweig widerstandsreziproke Konstruktionen angeben, wenn 
sogenannte ideale Transformatoren zugelassen werden. Gegenüber den Dualitätssätzen von 
Poincare und Alexander handelt es sich bei der hier auftretenden Dualität nicht um 
eine duale Beziehung zwischen zwei Teilen ein und derselben Mannigfaltigkeit, sondern um 
eine Dualität zwischen zwei getrennten Komplexen. 413 


Fine Verallgemeinerung des Differenzenverfahrens 
für Differentialgleichungen. 
Von Lothar Collatz in Berlin. 
ei dem zur numerischen Lösung von  Differentialgleichungsproblemen angewendeten 
Differenzenverfahren stellt man der Differentialgleichung, indem man die Differential- 
quotienten «dureh Differenzenquotienten ersetzt, ein System von Differenzengleichungen für 
Näherungswerte F der unbekannten Funktionswerte f an diskreten Stellen des Raumes, den 


°4) Kine von gemeinsamer Einbettung in eine Mannigfaltigkeit unabhängige Definition der Dualität gibt auch 
Il. Whitney in Transaetions Amer. Math. Soe. Bd. 34. Nr. ?, S. 339 bis 362. 

5) Der hier interessierende Teil des mengentheoretisch bewiesenen Satzes von Knratowski (Fundamenta 
Mathematieae Bd. 15 (1930), S. ?71 bis SB) kann unter Benutzung einer rein kombinatorischen Definition der Dualität 
rein kombinatorisch bewiesen werden. 
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Gitterpunkten“ gegenüber. Es hat sich nun bei der großen Zahl der nach diesem Verfahren 
lurchgerechneten Aufgaben gezeigt, daß das Verfahren bei großer Maschenweite h, also 
reringer Unbekanntenanzahl und relativ kurzer Rechnung, einen guten Überblick über den 
Verlauf der Lösungsfunktion gibt, daß aber das Verfahren bei Maschenverfeinerung nur langsam 
konvergiert. Braucht man genaue Näherungswerte F, so muß man sehr kleine Werte von h 
benutzen. Dabei steigt die Rechenarbeit stark an, und das Differenzenverfahren bleibt gegen: 
über anderen Lösungsmethoden nicht konkurrenzfähig. 

Die Konvergenz läßt sich beschleunigen, wenn man die Differentialquotienten nicht 
dureh Differenzenquotienten, sondern durch „finite Ausdrücke“ ersetzt, die genauer mit den 
Differentialquotienten übereinstimmen als die Differenzenquotienten. So ist z. B. für fünfmal 
stetig differenzierbare Funktionen f(x) 


df f(x >h) 2 


— +J 
dr 12h IS max in (x—?2h,x 
h 2 max in (x, x -+h) 


Die Abweichung des erstgenannten finiten Ausdrucks vom Differentialquotienten geht bei 
Maschenverfeinerung mit der vierten Potenz von h gegen 0, die des Differenzenquotienten 
(dagegen nur proportional h. Es ist daher beim Rechnen mit derartigen stärker approximierenden 
finiten Ausdrücken bei Übergang zu kleinerem h eine raschere Konvergenz zu erwarten, als 
beim Reehnen mit Differenzenquotienten, was sich in zahlreichen Fällen bestätigt. hat. 


Ist eine Funktion f(#,.%s....,%,„) durch eine lineare Differentialgleichung L(f) = 
as as.) von mter Ordnung und durch lineare Randbedingungen festgelegt, so legen wir 
der Berechnung der Funktion f ein System von (nicht notwendig regelmäßig angeordneten) 
Punkten 1, 2, ‚p im n-dimensionalen Raum zugrunde. Die Aufgabe des Differenzenver- 
fahrens ist die Bestimmung eines Systenis von Näherungswerten F\, F\,,..., F, für die Werte 
fss f» der Lösungsfunktion f an den Punkten 1, 2, ...., p. Für die F, wird ein System 
linearer Gleichungen aufgestellt, von denen man eine so findet: Man greift einen der Punkte, 

etwa den Punkt j, heraus und bildet für ihn eine Summe der Gestalt I’ C.9 fo. Jeder 
einzelne dieser p Summanden wird nach dem Taylorschen Satz entwickelt nach f und den 
partiellen Ableitungen von f an der Stelle j bis zu Gliedern rter Ordnung (r = m); dann wird 
die Summe umgeordnet und zusammengefaßt nach f und den partiellen Ableitungen von f 
an der Stelle j. Diese Entwicklung wird verglichen mit ZL(fj), und die gliedweise Gleich- 
setzung der Koeffizienten der partiellen Ableitungen von f an der Stelle j liefert lineare Be- 
stimmungsgleichungen für die noch verfügbaren Konstanten 0%. Es gilt dann 


IN fo — Restglieder (r + I)ter Ordnung 


Dann soll Z Co‘ Fo=t; eine der Gleichungen für die F; sein. Die Frage nach der Auf- 
lösbarkeit Gleichungssystems für die Cor erledigt der Satz: zu einem beliebig 
vorgegebenen linearen Differentialausdruck L(f) finite Ausdrücke I Co fo mit jedem ge- 


wünschten Approximationsgrad r. Für die meist gebrauchten Differentialausdrücke sind finite 
Ausdrücke höherer Approximation bereits von mir berechnet worden. Weitere Gleichungen 
für die F; erhält man, indem man andere Systeme von Konstanten Co‘) verwendet, oder an 
Stelle des Punktes j andere Punkte bevorzugt, oder indem man die linearen Randbedingungen 
heranzieht. 

Die Fehler F, — f; der Näherungen genügen einem Gleichungssystem mit derselben 
Koeffizientenmatrix wie bei dem Gleiehungssystem für die Näherungen F'; selbst; es folgt hier 
allgemein aus der Möglichkeit, das Differenzenverfahren überhaupt anwenden zu können, schon 
(die Möglichkeit, eine Fehlerabschätzung für das Verfahren durchzuführen, vorausgesetzt, daß 
die Lösungsfunktion genügend oft, nämlich (+) mal partiell nach allen Variablen differen- 
zierbar ist: Eindeutigkeit der Lösung braucht nicht vorausgesetzt zu werden. 

Eine ausführlichere Darstellung werde ich demnächst veröffentlichen. 4A 
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Spannungszustand durch Drillung in dünnwandigen Kastenträgern bei 
verhinderter Endwölbung. 


Von Hans Ebner in Berlin-Adlershof. 


ei der Beanspruchung von Stäben mit nieht achsensymmetrischem Querschnitt durch ein 
Drillmoment ergibt sich bekanntlich eine Querschnittsverwölbung. Wird diese Verwölbung 
an irgendeiner Stelle des Stabes behindert, dann überlagert sich dem reinen Schubspannungs- 
zustand, der bei Drillung mit freier (Querschnittswölbung entsteht, ein zusätzlicher Längs- 
und Schubspannungszustand, der in der Nähe des an der Wölbung behinderten Querschnitts 
von maßgebendem Einfluß sein kann. 
Bei dünnwandigen Hohlstäben mit dem umschlossenen (Juerschnitt F; und der Wand- 
dieke d ergibt sich der Schubspannungszustand durch ein Drillmoment M bei freier Quer- 
schnittswölbung aus der Bredtschen Formel: 


M 
2 


Der bei verhinderter Endwölbung sich diesem Schubspannungszustand überlagernde Spannungs- 
zustand wird im folgenden für den Sonderfall des einseitig eingespannten Kastenträgers mit 
rechteckigem (Querschnitt ermittelt. Der Kastenträger hat die Länge a, die Querschnitts- 
seiten b und e; die entsprechenden Wanddicken sind d, und d.. Jede Längsscheibe hat ihr 
besonderes Koordinatensystem (Nullpunkt in Mitte der @Querwandseiten am Einspannquer- 
schnitt, Abszissenachse .- in Längsriehtung. Ordinatenachse y bzw. z in Richtung der Quer- 
wandseiten b bzw. e). 

Bei reiner Drillung findet eine Verwölbung des ganzen Kastens in der Weise statt, daß 
sich die 4 Eekpunkte eines vorher ebenen (uerschnitts gegeneinander aus dieser Ebene ver- 
schieben, wobei die Querschnitte der einzelnen Längsscheiben zwischen den Eekpunkten 
eben bleiben. Diese geradlinige Verwölbung muß am eingespannten Ende durch Aufbringen eines 
Längsspannungszustandes rückgängig gemacht werden. Die Wirkung dieses Selbstspannungs- 
zustandes läßt sich durch eine Gruppe X von 4 gleichgroßen, abwechselnd entgegengesetzt 
gerichteten Längskräften in den Ecken des Kastenträgers ersetzen. Ist der Spannungszustand 
für die Kraftgruppe X = 1 ermittelt, so kann die statisch unbestimmte Größe X mit Hilfe 
des Prinzips der virtuellen Verrückungen aus der Bedingung verschwindender Verwölbung 
am Einspannquerschnitt berechnet werden. 

In den einzelnen Längsscheiben wird ein ebener Spannungszustand vorausgesetzt. Die 
Gleichgewichts- und Verträgliehkeitsbedingungen für die Scheibenelemente sind erfüllt, wenn 
die Spannungen als Ableitungen der Airyschen Spannungsfunktion F gewonnen werden, 
wobei F der Differentialgleichung I1F=0 genügen muß. Außer dem Gleichgewicht 
am Element muß Gleichgewicht der einzelnen Längsscheiben vorhanden sein. Da bei Wirkung 
der Längskraftgruppe X 1 von den Endquerwänden auf die Längswände der Randkräfte 
b)2a bzw. e/2a übertragen werden, müssen die Resultierenden der Schubspannungen im 
(uerschnitt der Längswände mit diesen Randkräften übereinstimmen. 


Außerdem sind folgende Randbedingungen zu erfüllen: 


I. Am eingespannten Rand «=0 muß der Querschnitt der einzelnen Längsscheiben 
eben bleiben. 

2. Am freien Rand «= a muß die Längspannung o, verschwinden. 

>». An den Längsrändern müssen die in @uerrichtung auftretenden Spannungen 
6, bzw. verschwinden. 

I. Die Randdehnungen &, der Längsscheiben müssen an den Längsrändern überein- 
stimmen; da dort o,= 0.0, kann dafür auch die Bedingung gleicher Rand- 
spannungen o, gesetzt werden. 

>. Die Schubbelastungen r:(d der Längsscheiben müssen am gemeinsamen Rand 
eleich sein. 

Die ersten beiden Randbedingungen werden im Ansatz der Spannungsfunktion berück- 
sichtigt, die drei letzten werden zur Bestimmung der in der angesetzten Spannungsfunktion 
noch enthaltenen Konstanten benutzt. Die in dieser Weise erhaltene Spannungsfunktion 
lautet dann z. B. für die Längsscheibe b: 


+ D, | B, Zin „ + „ Coj km cos 
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Daraus ergibt sich dann für die Längsscheibe b infolge X=1 der Spannungszustand: 


] 


Darin ist: 


1 Zin?kne, 1 


2kmd, Rat, 


mt 


Der Spannungszustand in der Längsscheibe e ergibt sich entprechend. 

Mit dem Spannungszustand X 1 kann jetzt die statisch unbestimmte Größe X aus der 
Klastizitätsbedingung bestimmt werden, daß die gegenseitige Verschiebung der Eckpunkte 
des Einspannquersehnitts infolge des Zustandes X entgegengesetzt gleich der entsprechenden 
Verschiebung infolge des Drillmomentes M bei unbehinderter Quersehnittswölbung sein muß. 
Die Eekverschiebungen infolge X=1 erhält man am zweekmäßigsten mit Hilfe des Prinzips 
der virtuellen Verrückungen, wobei die Arbeit der von den Spannungen an den entsprechenden 
Dehnungen für die Gesamtheit der Scheibenelemente aufzustellen ist. Näherungsweise kann 
die Bestimmung der Eekverschiebungen infolge X 1 unter der Annahme linearer Verteilung 
der Längsspannungen erfolgen. Bei der Berechnung eines Kastenträgers mit praktisch vor- 
kommenden Abmessungen ergab sich bei dieser Näherung die statisch Unbestimmte rd. 
Id v. H. zu groß. 

Der allgemeinere Fall des auf Drillung beanspruchten Kastenträgers mit Zwischenquer- 
wänden und beliebigen an den Querwänden angreifenden Drehmomenten wurde unter der 
Annahme linearer Verteilung der Längsspannungen in einer in der ZFM 1953, Heft 23 und 24 
veröffentlichten Arbeit untersucht. 487 


Zerreißversuche an Flachstäben mit ebener Formänderung. 
Von H. Fromm in Danzig. 


ie Erforschung der Werkstoffeigenschaften drängt dazu, auch für plastische Stoffe die 
„technische Kohäsion* zu ermitteln, d. h. die Normal- oder Schubspannungen, welche 
zu einem Reiß- bzw. Scherbruch führen. Bei Zerreißversuchen an Rundstäben kann man aus 
Kraft und engstem (Querschnitt der Einsehnürung eine mittlere Bruchspannung errechnen. 
Die Einschnürung bedingt jedoch wie eine Kerbe eine ungleiehmäßige Spannungsverteilung, 
und zwar kann man einsehen, daß im plastischen Zustand die größte Spannung in der Mitte 
und die kleinste am Rande herrscht. Dies Ergebnis findet durch einen bekannten Versuch 
von Ludwik') seine Bestätigung. Die Ermittlung der Kohäsion fordert die Kenntnis dieser 
größten Spannung. 

Für das ebene Problem ist die Lösung dieser Aufgabe weitgehend vorbereitet. Auf 
Grund der Henekyschen?) Sätze über Gleitliniennetze hat Prandtl’) eine graphische 
Konstruktion des Netzes aus zwei gegebenen Gleitlinien angegeben, und mit Hilfe der von 
Caratheodory und Schmidt?) gegebenen analytischen Behandlung hat Schapitz?’) gezeigt, 

a) P. Ludwik: Zeitschr. d. Ver. Deutsch. Ing. 71, 1927, S. 1592. 

2) H. Henky: ZAMM 3, 1933, 8. 241. 

3) Prandtl: ZAMM 3, 1925, S. 

A) C, Caratheodory u. E. Schmidt: ZAMM 3, 1923, S. 46». 

5), E. Schapitz: Ing. Arch. 4, 1933, S. 227. 
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daß das Netz durch eine freie Oberfläche völlig bestimmt ist, und hat diese Randwertaufgabe 
gelöst. Davon benutzen wir folgende Ergebnisse: 1. Zwei von den Endpunkten A und B 
eines einsinnig um den Winkel 9 gekrümmten Randstückes ausgehende Gleitlinien krümmen 
sich bis zu ihrem Schnittpunkt € entgegengesetzt um den gleichen Winkel #9 = g/2 (Abb. 1). 
2. Eine von einem Wendepunkt des Randes ausgehende Gleitlinie schneidet die Gleitlinien 
der anderen Schar in Wendepunkten. Daraus folgt für die Einschnürung wegen ihrer doppelten 
Symmetrie, daß die Gleitlinien durch den engsten Querschnitt vom konkaven Teil des Randes 
ausgehen und daß ihr Krümmungswinkel 9 dem Winkel zwischen entsprechender Randtangente 
und Stabachse gleich ist (Abb. 2). Die Gleitlinie vom Wendepunkt des Randes geht also 
dureh die Querschnittsmitte, der somit der größte Winkel #9 zukommt. Es ist dies ge- 


INAX 


Abb. 1. Abb. 2. 


wissermaßen das Gesetz, nach dem sich die Form des Randes in Anpassung an die Dicke 
im engsten (Juerschnitt ausbildet. 3. Die Hauptspannungen o, und o, in der Netzebene wachsen 
vom Rande aus mit 9 je um 2kd, wo 2k = 0, — o, die konstant vorausgesetzte Hauptspannungs- 
differenz ist. Am Rand ist nun mit 0,0 überall o, = 0, 2k, und im engsten Querschnitt 
gilt daher 


Is bestätigt sich also, daß die größte Spannung in der Mitte liegt, sie ist 
— 


Die Voraussetzung der Henekyschen Sätze: o, 0, 2k kann in verschiedener Weise 
physikalisch bedingt sein: Sie kann erfüllt sein, wenn bei ebenem Spannungszustand (o, = 0) 
die Mohrsche Fließbedingung gilt, daß die größte Hauptspannungsdifferenz gleich 2k ist. 
Dazu müssen aber o, und o, verschiedene Vorzeichen haben. In der Einschnürung ist nun 
der Rand konkav, also müßte 9 (wie in Abb. 2) positiv, also nach (1) o, und o, positiv sein: 
d. h. die ganzen Überlegungen gelten hier nieht. In der Tat treten auch bei dünnen breiten 
Flachstäben schräg durch den Stab verlaufende Einschnürungen auf, die mehr zu einer Be- 
ziehung 6, 0, 2% passen. — Die Henekysche Voraussetzung ist aber stets bei ebener 
Verzerrung erfüllt, wenn dazu die v. Misessche Fließbedingung und Fließregel gilt. Für 
diesen Fall gelten also die Gleiehungen (1) und (la). Es genügt z. B. eine Messung des 
Wendetangentenwinkels um Zu finden. Hat man die ganze Rand- 
kurve von einem Versuch aufgenommen, so läßt sich das Gleitliniennetz zeichnen und damit 
die Spannungsverteilung ermitteln. Mit der mittleren Spannung oy im engsten Querschnitt f 
ist dann die Zugkraft P=oy f \o,df=or-\(1+9)df, also nach (la) 


+ 
max — \ (1 d f 


(2). 
Die Spannungsverteilung kann man mit guter Näherung nach einer Parabel annehmen °). 
Dann ıst 


1+9. 


3 


6 
’ 


Diese Ergebnisse setzen homogenen Stoll voraus, Verfestigung macht sie nicht ungültig. 
solange sie homogen ist und die v. Misesschen Annahmen bestehen läßt. Gerade in der 
Kinsehnürung ist aber die Verfestigung nicht homogen, sondern in der Mitte erheblich stärker 
als am Rande, so daß die Henekvyschen Sätze nieht mehr gelten. Das Gleitliniennetz sieht 
anders aus und (1) liefert zu kleine Werte, wenn #9 jetzt den Tangentenwinkel des Randes 


Verl. Fußnote 7. 
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bedeutet, der einem Querschnittspunkt durch seine Gleitlinie zugeordnet ist. Da nun Ver: 
'estigung und # beide nach der Mitte des Querschnitts hin wachsen, so wird man die Ver- 
'estigung in erster Näherung dadurch berücksichtigen können, daß man in (I) 9 durch n9 >49 
ersetzt. So folet aus (2a) 
1 n 


MAX 


Um nun diese Möglichkeit zur Abschätzung der größten Spannung und damit der 
technischen Kohäsion auszunutzen, hat Herr Baranskı auf meine Anregung im Institut für 
Mechanik der Technischen Hochschule Berlin Zerreißversuche mit ebener Formänderung aus- 
oeführt. Flachstäbe von Ss mm Dicke und bis zu 100 mm Breite wurden in der Zerreißmaschine ge- 
(lehnt, wobei die Breite dadurch unveränderlich erhalten wurde, daß die längsverlaufenden Ränder 
dureh T-Flansche verstärkt und dureh dazwischen gebaute Querstützen an der gegenseitigen 
Annäherung verhindert wurden’). Daß die Formänderung mit guter Vollkommenheit eben war, 
konnte mehrfach geprüft werden: 1. Ein aufgeritztes Quadratnetz ging in ein Rechteeknetz 
über, wobei die Längsrisse ihren Abstand von 5 mm beibehielten, auch im Tal der Ein- 
schnürung, welches senkrecht zur Zugrichtung geradlinig über die ganze Breite verlief. 2. Den- 
selben Verlauf zeigt auch der Bruch, der bei Messing 63 (ebenso wie beim Rundstab) ein 
Scherbruch ist, während bei Mn-Si-Stahl „Böhler ES* in der Mittelzone ein Reißbruch ansetzt, 
dem ein Scherbruch in den Außenteilen des Querschnitts folgt (genau wie beim Rundstab). 
Der Bruch beginnt in der Mitte des Tales und läuft im Tal dem Flanschen zu; daher ist 
die für den Bruch gemessene Spannung wirklich die den ebenen Bruch erzeugende Spannung. 
3. Laufend mit einem besonderen Gerät durchgeführte Diekenmessungen über die ganze Länge 
und Breite des Stabes ließen eine gleichmäßig verteilte Wellung erkennen, wie sie auch an 
Rundstäben *) als Folge der Inhomogenität des Stoffes beobachtet wurde. Wie an Rundstäben 
konnte auch hier die allmähliche Entstehung der ebenen Einsehnürung messend verfolgt werden, 
wobei sich im wesentlichen die von Schapitz gefundenen Ergebnisse wiederholen lassen, 
I. Die unten beschriebene Art der Spannungsermittlung liefert eine weitere Probe für die 
Ehbenheit der Formänderung. 

Bei dem Böhlerstahl ergab sich nun nach dem Bruch ein Wendetangentenwinkel von 
Fa > 924 (> 149%). Ohne Beachtung der Verfestigung würde also nach (la) und (2a) 


MAX 
OÖ max 1,24 1,07 


d. h. die technische Reißkohäsion würde nur etwa 7"o über der über den (Querschnitt 
gemittelten Spannung liegen. Demgegenüber wird man also den Einfluß der Verfestigung 
nicht vernachlässigen dürfen, der für den vorliegenden Fall ungefähr richtig geschätzt sein 
dürfte, wenn man in (3) n=1,5 setzt. So wird 


Bei Messing 63 sind die Unterschiede noch geringer. — Selbst wenn #:%,., | würde, wäre 
erst 9, =1.2:0y. Nun ist nicht zu erwarten, daß die Spannungsverteilung beim Rundstab 


sehr viel ungleichmäßiger ist. Man kann daher die Reißkohäsion um etwa 10° höher als 
die gemittelte Spannung schätzen, wenn Reißbruch vorliegt, und die Scherkohäsion bei Scher- 
bruch zu etwa 0,55 - oy. 


| Hl-Gebref Zugunter 
| N | | | [ebene Ver Ag/mm? 
pr A form, ) | 
| 
45 
| 
Sen Ränder des Hlackstabes 30! 
S 30! N 
S 20) 
S 
| 
155% Betrag des Verformungsdenators D 
R51013.7) Dehnung & SHlat y2 u, 2 
Abh. 3. Messing. Abb. 4. Messing. 


?), Über die Entwieklung der günstigsten Stabform, die Versuchseinriehtung und die Versuchsergebnisse vgl. die 
Dissertation von G. Baranski, Zeitschr. f. Metallkunde 26, 1034. S. 173. 
s) Vgl. die im gleiehen Institut durchgeführten Versuche von E. Sehapitz: Zeitsehr. f. Physik 70, 1031, S. 641. 
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Um die absolute Größe der Spannungen zu bestimmen, wurde für mehrere Stäbe von 
verschiedener Breite, aber genau gleicher Randversteifung, die Zugkraft in Abhängigkeit von 
der Dehnung gemessen. Die Kraftdifferenz für gleiche Dehnung verschiedener Stäbe ist die 
auf «lie Breitendifferenz entfallende Kraft zur ebenen Verformung; Division durch Breiten- 
differenz und Dicke gibt die gemittelte Spannung. Ihre Auftragung über der Dehnung liefert 
eine Darstellung der Verfestigungskurve (Abb. 3). Andere Darstellungen ergeben sich durch 
Auftragung der Intensität”) des Spannungsdeviators über der des Verzerrungsdeviators (Abb. 4) 


- 


oder über der plastischen Arbeit (Abb. 5). In diesen letzten Auftragungen fand Schmidt’) 


10 | 
20% 
last. Arbeit an der Volumeneinheit A ? dt min Dehnung & 
Abb. 5. Messing. Abb. 6. Stahl. 


eine Übereinstimmung verschiedenartiger Versuche, bei denen er jedoch nur bis etwa 15" 
Verformung ging. Abb. 4 und 5 zeigen, daß Unterschiede zwischen ebener Verformung und 
einachsigem Zug hier noch in der Fehlergrenze liegen, bei größeren Verformungen aber 
deutlich werden. Dies zeigt auch die a-Kurve, die bei Deckung der Verfestigungskurven in 
Abb. 4 und 5 konstant a=15,5°/, (wie bei Gültigkeit der v. Misesschen Fließbedingung) 
sein müßte; sie nähert sich aber dem Wert a=-0 (wie bei Gültigkeit der Mohrschen Be- 
dingung). Dies könnte dazu verführen zu sagen, daß die Fließbedingung von Messing mit 
zunehmender Verfestigung von der v. Misesschen in die Mohrsche Form hinüberwechselt, 
ähnlich wie Lode") für Flußeisen ein Wechseln von der Mohrschen in die v. Misessche 


HL, 
N 25 
20 
| 
15 | 
70 %% N 
0 Fer. 75 2 % 23 
Betrag des Verformungsdeviarors D N 
nr, DW oz 95 
Abb. 7. Stahl. Abb. 8. 


Form ausgesagt hat, also in umgekehrter Richtung, wie sie auch beim Böhlerstahl gefunden 
wurde (Abb. 7). Eine solehe Aussage täuscht aber insofern etwas über den wahren Sach- 
verhalt, als die Darstellung verschiedenartiger Versuche in verschiedener Weise beim Wechsel 
der als Abszisse gewählten Größe verzerrt wird; daher nehmen auch die a-Kurven völlig 
verschiedenartige Formen an. Man muß vielmehr anstreben, neue Größen für Abszisse und 
Ordinate zu finden, in deren Darstellung die Kurven voll zur Deckung kommen. Dann stellt 
die Abszisse ein geeignetes Maß der Vorgeschichte dar, und die Ordinate gibt eine unveränder- 
liche Form der Fließbedingung an. 
», Vel. U. Fromm: Ing. Arch. 4, 1933, S. #32 und ZAMM 13, 1933, S. 427. 


10), R. Schmidt: Ing. Arch. 8, 193°, S. 215. 
1, Lode:Der Einfluß der mittleren Hauptspannung auf das Fließen der Metalle, Forsch. Ing. Wes., H. 303, 1928. 
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Diese noch bestehende Unsicherheit der Darstellung macht sich auch beim Vergleich 
des Bruches nach ebener Verformung und am Rundstab bemerkbar. Immerhin ergibt sich 
deutlich, daß die Bruchverformung beim Rundstab größer ist; so ergab z. B. die Messung 
der Brucheinschnürung beim Böhlerstahl: 


rund eben !?) 
Quersehnittsverminderung 67.6 56.3 0), 
mittlere Höchstdehnung . 208 129 
Deviatorintensität . . . Da 1.37 1.17 
Mittelspannung . ....0°0y 76.8 78.0 ke/mm? (Reißbruch) 
Wendetangente . . . . 0,28 (= 16°) = 0,24 (= 14°) 


Die Mittelspannung beim Bruch kann durch Extrapolation zur Quersehnittsverminderung 45 
ermittelt werden, auch wenn nur wenige Punkte hinter der Höchstlast gemessen sind, da der 
weitere Verlauf bei Darstellung der Spannung über der Quersehnittsverminderung der Tangente 
im Höchstlastpunkt nahe folgt und diese Tangente konstruiert werden kann, wenn der Höchst- 
lastpunkt (HL) festliegt'*) (Abb. 8). 

Nach (4) ist nun die Reißkobäsion bei ebener Verformung auf 1.1:78,0 kg/mm? 
zu schätzen. Beim Rundstab ist sie noch höher zu erwarten, da der größeren Verformung 
eine höhere Reißverfestigung entsprechen dürfte. Diesem Umstande ist es übrigens zu ver- 
danken, daß die Auswertung der Versuche überhaupt möglich war, da sonst der Bruch nicht 
im ebenen Mittelstück, sondern in der Randversteifung begonnen hätte, die sich wesentlich 
wie ein Rundstab verformt (vgl. auch Abb. 3). 

Die vorstehenden Mitteilungen sollen noch dureh weitere Versuche geprüft und ergänzt 
werden. Auch sollen nieht, wie bisher, nur einartige Versuche durchgeführt werden. Versuche 
mit Artwechsel sind in Vorbereitung. 475 


Das Zustandekommen der turbulenten Austauschgröße. 
Von Hans Gebelein in Stuttgart. 


eit Boussinesq setzt man die mittlere turbulente Strömung in Analogie zur laminaren 
S Strömung. Der Unterschied besteht darin, daß bei der turbulenten Strömung an Stelle 
der kinematischen Zähigkeit eine andersartige Änstauschgröße, die turbulente Austauschgröße, 
wirksam ist. Lautet z. B. für die einfache stationäre Rohrströmung die Gleichung der laminaren 
Bewegung 


du 

Yy 
so lautet die entsprechende Gleichung für die mittlere turbulente Strömung 

dü 


wobei @ die mittlere Geschwindigkeit bedeutet. Das Problem ist die Ermittlung der turbu- 
lenten Austauschgröße e, die ebenso wie die kinematische Zähigkeit die Dimension em?/see 
hat, die aber keine Konstante ist, sondern vom Ort in der Strömung und von der mittleren 
(Greschwindigkeitsverteilung abhängt. 

Bekanntlich machte Prandtl für die Größe & aus Gl. (2) den erfolgreichen Ansatz 


„au 
d y . . . . . . . . . . . . . (3), 


wo die Länge ! auf Grund anschaulicher Überlegungen in Analogie zur kinetischen Gastheorie 
als „Mischungsweg* bezeichnet wird. Tatsächlich ist dieser Ansatz demjenigen von Boussinesq 
wesentlich überlegen, da die Länge ! im wesentlichen nur eine Funktion des Ortes ist. Zwar 
hängt ! noch schwach von der Reynoldsschen Zahl ab, aber mit wachsendem Re nähert 
!(y) sich von oben her einer Grenzkurve, so daß für unendliche Reynoldssche Zahl der 
Mischungsweg exakt nur noch vom Orte abhängt. Die Prandtlsche Theorie macht keine 
Aussage über die Funktion, die diese Abhängigkeit vom Orte beschreibt. Aus der Erfahrung 
wissen wir nur, daß !(y) in Wandnähe proportional mit dem Wandabstand y ansteigt. 


12) Baranski gibt hierfür kleinere Werte an, die vor dem Bruch liegen. 
13, Vgl. z. B. G. Sachs: Mechan. Technologie der Metalle, Leipzig 1925. 8. 26. 
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In der neuen Theorie der Hydrodynamik auf der Grundlage der physikalischen Statistik ') 
tritt die turbulente Austauschgröße als gleichartiges Gegenstück zur kinematischen Zähigkeit 
neben dieser auf. Sie ist die mittlere Komponente ,, , eines Tensors, des „Turbulenztensors“, 
dessen Berechnung auf Grund einfacher hydrodynamischer Vorstellungen mit den Schlußweisen 
der Wahrscheinlichkeitsrechnung geschieht. Die Gleichung der stationären Rohrströmung 
lautet in dieser Theorie 


du 
+ „) d . (4). 
Der Aufbau ist für » und z,, der gleiche, nämlich 
| 
sen. 


wo c eine Geschwindigkeit und #* eine für die Statistik grundlegende Zeit, die Verweilzeit, 
darstellt. Der Faktor 1/3 rührt von den drei Dimensionen des Raumes her. Bekanntlich ist 


ce/, wo ce die Schwirrgeschwindigkeit der 


ce? 


nach der kinetischen Gastheorie — 


Moleküle bedeutet und f*: Z die mittlere Zeit ist, die die Moleküle ohne Zusammenstoß 
fliegen, d. h. die Zeit der mittleren freien Weglänge. 

Die physikalische Grundlage für die Entstehung der turbulenten Austauschgröße und 
damit der Turbulenzstatistik sind nun nicht die Vorstellungen der kinetischen Gastheorie, 
sondern einfache hydrodynamische Erfahrungstatsachen. Es ist dies die Tatsache, daß jede 
Strömung in einer zähen Flüssigkeit mit der Zeit zur Ruhe kommt, wenn nicht an den Wänden 
Schubspannung übertragen wird. Dies äußert sich in der dauernden Entstehung neuen Wirbel- 
materials. Die Achse der neuentstehenden Wirbel liegt im Mittel parallel zur Wand und 
senkrecht zur Strömungsriehtung. Bei der laminaren Strömung weist das neuentstehende 
Wirbelmaterial konstante Diehte auf, und die Riehtung der Wirbel ist exakt die angegebene. 
Bei der turbulenten Strömung dagegen ist diese Gleichförmigkeit durch elementare Vorgänge 
in der Grenzschicht, Wirbelablösungen an Zacken und dergleichen gestört, und eine Augen- 
blieksaufnahme der neuentstehenden Wirbelfäden ähnelt anschaulich der Maserung des Holzes, 
aus dem die Wand gefertigt sein mag, wenn die Fasern senkrecht zur Strömungsrichtung 
liegen. Diese zufallartigen Abweichungen der neuentstehenden Rotation vom Mittelwert sind 
die Grundlage der Turbulenztheorie und geben den Anlaß zur Entstehung der turbulenten 
Austauschgröße. Natürlich können im Rahmen eines Vortrags nur die Hauptgedanken mit- 
geteilt werden. 

Die nächste Aufgabe der Theorie ist die Ermittlung der Größen e und #*. Zunächst 
soll das. Zustandekommen der Geschwindigkeit e erklärt werden. e ist der Mittelwert der zur 
mittleren Geschwindigkeit # an der betrachteten Stelle infolge der turbulenten Schwankungen 
hinzutretenden Streuungsgesehwindigkeiten. In der laminaren Strömung ist e gleich Null. 
Diese Zusatzgeschwindigkeiten rühren davon her, daß die Rotation an jeder Stelle der Flüssig- 
keit dauernd um ihren Mittelwert schwankt. Diese fehlerhaften Größen der Rotation sind 
es, die im Aufpunkt die fehlerhafte, zur mittleren Geschwindigkeit hinzutretende Streuungs- 
voeschwindigkeit erzeugen. Dieser Gedanke liefert bei der exakten Durchreehnung mit den 
Hilfsmitteln der Wahrscheinlichkeitsreehnung das Ergebnis, daß die mittlere Zusatz- 
vseschwindigkeite proportional der vierten Wurzelaus dem Absolutwert der 
mittleren Rotation an der betreffenden Stelle ist: 


Das Zustandekommen der Verweilzeit f* ist nun in folgender Weise zu übersehen. Die 
Verweilzeit ist ein Maß dafür, ob der Strömungsvorgang an der betreffenden Stelle mehr oder 
weniger lange so erfolgt, wie es die deterministische Mechanik auf Grund der Anfangsbe- 
(dlingungen voraussagt. An den Wänden, an denen Schubspannung übertragen wird, entsteht 
zufallartig Wirbelmaterial, das in der Flüssigkeit nicht deterministisch erfaßbare, infinitesimale 
Zusatzgeschwindigkeiten hervorruft. Diese überlagern sich dauernd den deterministisch be- 
stimmten Geschwindigkeiten und führen zur allmählichen Auflösung des Determinismus und 
zu einer endlichen Verweilzeit. Sind keine Wände vorhanden, an denen Schubspannungen 
übertragen werden, so ist die Verweilzeit unbegrenzt groß. Wie die ausführliche Rechnung lehrt, 
ist die Verweilzeit # proportional dem Quotienten aus der kinetischen Energie, die in den 
deterministisch bestimmten Streuungsgeschwindigkeiten mit dem Mittelwert e steckt und der 


!, HM. Gebelein: Turbulenz. Im Verlag Julius Springer, Berlin, im Erscheinen begriffen. 
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kinetischen Energie in den zufallartigen, durch unregelmäßiges, neuentstehendes Wandwirbel- 
material erzeugten Zusatzgeschwindigkeiten. Diese Geschwindigkeiten aber sind um so kleiner, 


je größer die Entfernung der betrachteten Stelle von der Wand ist. Daher ergibt sich 


Aus den Beziehungen (6) und (7) folgt dann zusammen für die turbulente Austauschgröße 


Genauer liefert die Theorie im Falle der unendlich ausgedehnten ebenen Wand für den 
Tensor der turbulenten Austauschgrößen, den Turbulenztensor, den Ausdruck 


0 


Dabei sind € und a zwei Konstante, die die Theorie unbestimmt läßt, über deren Größe aber 
die Experimente erschöpfend Auskunft geben. Insbesondere ist a eine kleine Länge von der 
Größenordnung der Zackenhöhe der Wandrauhigkeiten, deren Größe für die Wirkung der 
Wandrauhigkeit entscheidend ist. Wir setzen sie im folgenden gleich Null. 

Es ist noch auf den Zusammenhang mit dem Prandtlschen Mischungsweg hinzuweisen. 
Die Austauschgröße & aus (3) ist nach (9) 


y\rotd'+r. 
du 
Beim ebenen stationären Problem ist |\rotb': “Fr und daher 
IW=Cyl/1+ . (li). 


Daraus ist zu ersehen, daß der Mischungsweg bei zunehmender Reynoldsscher Zahl von 
oben her gegen eine allein vom Ort abhängige Grenzfunktion strebt, die im einfachsten Fall 
der unendlichen ebenen Platte im Einklang mit der Erfahrung eine lineare Beziehung mit 


dem Wandabstand ist. 

Die neue Theorie liefert also ein Verfahren zur Berechnung des Turbulenztensors und 
damit auch der Abhängigkeit des Mischungsweges vom Ort. Der Turbulenztensor ist außer- 
dem nach dieser Theorie ganz allgemein ein Funktional der die Strömung begrenzenden Wände, 
an denen Schubspannung übertragen wird, und daher sind die Turbulenzgleichungen zum 
Unterschied von den Navier-Stokesschen Gleichungen andere und andere für die ver- 
schiedenen hydrodynamischen Probleme. 

Wir bringen nun zum Schluß eine kurze Tabelle der Turbulenztensoren für einige ein- 
fache Probleme. 


a) Strömung an der unendlich ausgedehnten ebenen Wand -—- Grundlage für die Grenz- 
schiehttheorie. (y= Wandabstand): 


00 
Hy) rotd 10 ) 


b) Strömung in der Entfernung x vor der Vorderkante einer unendlich ausgedehnten 


Halbebene — Grundlage für die Theorie des Freistrahls: 


VD) 0 0 
IH [VO 10 


2 
4 


ce) Strömung zwischen zwei parallelen Platten im Abstand 2b. (y gemessen von der Mittel- 
ebene aus): 


000 


b? — ı?)? 
II(ı er | [0 10 
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d) Strömung im Kreisrohr. (4° Abstand von Rohrachse mit der Einheit des Rohrradius r): 


00 
Il (ıy) tv [oO 1 0 
00 

f: 


Sr? \ -2yeosg)": Sr" \ —2ycosgp)’:' 


> 


speziell ıst 


a) in der Nähe der Wand theoretischer Verlauf des 
Mischungswegs im Kreisrohr 


- m) rotd (O1 (vgl. Beispiel a), 


5) in der Rohrachse 


Rohrachse, 


000 
= rot5 [0 10). 


Der theoretische Verlauf des Mischungswegs ist für diesen 
Fall aus der Abbildung zu ersehen. 476 


Zur Zähigkeit wäßriger Suspensionen'). 
Von R. Hermann in Aachen. 


ie von seiten der Keramiker über die Zähigkeit von Suspensionen aus Ton, Kaolin (Por- 
D zellanerde) oder Zementrohmehl vorliegenden zahlreichen Untersuchungen sind hin- 
siehtlich der benutzten Meßmethoden vom hydrodynamischen Standpunkt recht wenig be- 
[riedigend. Es wurde deshalb eine Meßmethode entwickelt (Ausfluß ‘aus einem langen, 
dünnen Glasrohr), die hydrodynamisch einwandfrei eine absolute Bestimmung der Zähigkeit 
von Tonsuspensionen erlaubt und den besonderen dabei vorliegenden Verhältnissen Rechnung 
trägt. An 3 wäßrigen Tonsuspensionen von 10 gewichtsprozentiger Konzentration in zwei 
Rohren verschiedenen Durehmessers bei starker Änderung der Druckhöhe bzw. der Durch- 
Hußmenge sowie der Reynoldsschen Zahl (1,4 bis 980) wurden die folgenden Ergebnisse 
gewonnen. 

Wäßrige Tonsuspensionen besitzen keine Zähigkeit im üblichen Sinne, sie befolgen 
den Newtonschen Schubspannungsansatz nicht. Der nach dem bekannten Hagen- 
Poiseuilleschen Gesetz berechnete Zähigkeitswert (Mittelwert über den Querschnitt) bzw. 
die Fluidität ist keine Konstante, sondern eine Funktion der Schubspannung an der Rohr- 
wand. Die Fluidität nimmt mit abnehmenden Schubspannungen ebenfalls ab und fällt bei 
den kleinsten Schubspannungen auf einen Bruchteil ('/s bis '/s) Ihres Wertes bei größeren 
Schubspannungen. Bei großen Schubspannungen nähert sich ihr Verhalten den homogenen 
Flüssigkeiten (konstante Fluidität).. Auch bei gegebener Schubspannung ist die Fluidität 
keine Konstante, sondern nimmt mit zunehmendem Alter dauernd ab („Ansteifung*). Nach 
'/; jähriger Alterung scheint die Ansteifung mit etwa der 2- bis 3fachen Zähigkeit des Aus- 
eangsstadiums ihr Ende erreicht zu haben. Die Suspensionen haften an der Rohrwand ge- 
nau wie die üblichen homogenen Flüssigkeiten, es tritt keine Gleitung auf. Extrapolation 
der gemessenen Kurven bis zur Schubspannung Null machen folgendes Verhalten wahr- 
scheinlich. ‚Junge Tonsuspensionen bleiben bis zu den kleinsten Schubspannungen, bis zum 
Ruhezustande, „flüssig“, wenn auch mit sehr hoher Zähigkeit, nämlich bis zur 12fachen Zähig- 
keit der bei großen Schubspannungen. Altere Tonsuspensionen dagegen sind unterhalb 
einer Schubbeanspruchung von I bis 2 dyn/em? nicht mehr „flüssig“, ihre Zähigkeit ist dort 
unendlich groß. Sie beginnen erst nach Überschreitung dieser Grenze („Fließgrenze*) zu 
ließen. 


I, Die Arbeit wurde ausgeführt im Institut Für Gesteinshüttenkunde der Techn. Hochschule Aachen mit freund- 
lieher Unterstützung der Notzemeinsehaft der Deutschen Wissenschaft, der ieh auch an dieser Stelle 
meinen Dank ausspreche. Die ausführliehe Veröffentliehung erfolgt in „Forschung auf d. Geb. d. Ing.-Wesens‘, 
voraussiehtlieh im ‚Januar 1935, 
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Schwingungen elastischer Seile. 
Von Rudolf Höger VDI in Prag. 


Ä 


ie Untersuchung befaßt sich mit den Schwingungen eines elastischen Seiles, welches 
D zwischen zwei gleichhoch gelegenen Punkten (Abstand 2 W) aufgehängt ist. Während 
‚lie bisher vorliegenden Arbeiten nur Verschiebungen der Seilpunkte in der vertikalen Rich- 
tung berücksichtigen oder in Anlehnung an die Theorie der schwingenden stark gespannten 
Saite den infolge der Schwingung auftretenden Spannungszuwachs vernachlässigen, sollen 
hier ohne diese einschränkenden Voraussetzungen die Verschiebungen r und f der Seilpunkte 
in radialer und tangentialer Richtung untersucht werden. 

Da der Durchhang h gegenüber der Spannweite 2 W als klein angenommen wird, möge 
die Kettenlinie durch einen Kreisbogen ersetzt werden. Zur Aufstellung der das Eigen- 
wertproblem charakterisierenden Differentialgleichung wird vom Hamiltonschen Prinzip 
auszegangen. Man gelangt so zunächst zu zwei simultanen Differentialgleichungen 2. Ord- 
nung, die auf eine Differentialgleichung 4. Ordnung zurückgeführt werden. Diese hat die 
Form 

or +zor" +20 ko? 


Hierin ist r der nur vom Orte abhängige Faktor der Radialverschiebung r; die Ableitungen 
beziehen sich auf den Neigungswinkel der Seiltangente. *, bezeichnet die bereits im Ruhe- 
stand vorhandene Dehnung, x hingegen eine von den mechanischen und geometrischen Eigen- 
schaften des Seiles (BElastizitätsmaß, Querschnitt usw.) abhängige Konstante. «© bedeutet die 
Kigenfrequenz. Eine analoge Gleichung besteht für r. 

Die Lösungen für r und + haben daher die Form einer Kombination von Kreis- und 
Hyperbelfunktionen. Die Erfüllung der Randbedingungen führt endlich zu zwei transzenden- 
ten Frequenzgleichungen, deren Lösungen zwei Reihen von Eigenwerten liefern. Man gelangt 
dergestalt zu einem System von Eigenschwingungen mit Schwingungsformen, welche bezüglich 
der Mittellotrechten symmetrisch sind und zu einem zweiten Schwingungssystem mit un- 
symmetrischen Schwingungsformen. 

Anschließend wurden noch die infolge der Schwingung auftretenden dynamischen Zusatz- 
spannungen berechnet, welche für die technischen Anwendungen (elektr. Fernleitungen, Seil- 
bahnen usw.) bedeutungsvoll sind. 475 


Finwirkung von Bodenschwingungen auf Pfeiler. 
Von Konrad Ludwig in Hannover. 


m Jahre 1933 wurde in dem Forschungsinstitut für Straßenbau von Herrn Prof. Dr.-Ing. Risch 
T:; der Technischen Hochschule in Hannover die Dissertation von Herrn Dipl.-Ing. Otto Brötz., 
Erschütterungen in Gebäudeteilen aus Backsteinmauerwerk (Mauerwerkspfeiler), abgeschlossen. 
Die Bodensehwingung, die als einfach harmonisch und ungedämpft vorausgesetzt wird, wird 
in eine Horizontal- und eine Vertikalkomponente zerlegt. Die Horizontalkomponente erregt 
im elastisch vorausgesetzten Pfeiler Biegungsschwingungen und die Vertikalkomponente 
Dehnungsschwingungen. In der Dissertation von Herrn Brötz wird die Theorie auf die er- 
zwungenen Schwingungen beschränkt. Ich erweitere diese Theorie durch die Beantwortung 
folgender aus der Praxis sich ergebenden Fragen: 

Wie sind die Eigenschwingungen beschaffen 

Wie sind die Resonanzschwingungen beschaffen ? 

Wie sind die erzwungenen, freien und Resonanzschwingungen beschaffen, wenn das 
obere Pfeilerende nicht frei ist, sondern eine Einzelmasse trägt? 

Die Differentialgleichung der Dehnungsschwingungen ist 


> 
0x of 


Ist das obere Pfeilerende frei, so sind die Randbedingungen 


= 


die Eigenfrequenzen 


n=1,2,3,... 


Ztschr. f. angew. 
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und die gesamte Schwingung 


n . 
sinwot—2—- vo °’ n (-1)" 
COS 


Wird die Bodenfrequenz einer Eigenfrequenz (der j) gleich, so wachsen die absoluten 
Amplituden der erzwungenen und der betreffenden freien Schwingung über alle Grenzen. Die 


dureh Grenzübergang oder den Ansatz 
D;(x) tcos w;t+ E; sin 


zu berechnende Resonanzschwingung ist 


u (— 1) cos -teosw;t— (- —(xsin cos |sin 
. 
cos, sin 
h — dr 


Trägt das obere Pfeilerende die Masse m, so wird die Randbedingung für das obere Ende 


F (F Inhalt des Pfeilerquerschnitts) 


und die Frequenzengleichung 
me 


eds 


Für die Amplituden .1, gilt die belastete Orthogonalität 


mn" 
Die Anfangsbedingung 
Q 
erfordert die Entwicklung 
oh 
F Esın + ın m COS 
x > . 
o,„B,.sın sın COS 
oh . oh 


FE cos „meosm 


Die zu entwickelnde Funktion erfüllt die Randbedingung für das obere Ende 


(Ay)x FE (A n)x 


A, eliminiert 


wenn hierin der Parameter »,° mit Hilfe der Differentialgleichung »’ A," 
Daher darf die Entwieklungsgleichung für das obere Ende zur Bestimmung der Ent- 
Demnach ist 


wird. 
wieklungskoeffizienten neben der belasteten Orthogonalität verwandt werden. 


die gesamte Schwingung 


FEsın + ım m cos 
sın COS sın 
oh " " 
FE cos 
») 
+ sin ""sino,„f 
+5, sın 


® | 
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nd die Resonanzschwingung bis auf die übrigen Eigenschwingungen 


" 
. 50; h . 509; h 
h+ 5. sin 3. sin 
0; | 2m; 
\ 
luten + cos, |sinmjt. 
. 50; ) 
Die h+ sin 
\ 
Ist das obere Pfeilerende frei, so ist die gesamte Biegungsschwingung 
sin Sin — cos sin, ei, +cos > 
/ h er / 
= | a Coj .) sın A, sın .) Coj A, sın un n 
h n a)" . PN n h A, n 
sin’ 5 
Ende 
cos? 
wobei die Frequenzengleichungen 
sind. Im Falle » = 1; ist bis auf die übrigen Eigenschwingungen 
„Ash 
2 „Ash x Bindıza 
Co) —; sin 5 
| sin | 
COS ja . 
+ sınoıjf. 
Arzt 
tcos 
Trägt das obere Pfeilerende eine Masse mit dem Trägheitsmoment ©, so ist die belastete 
Orthogonalität 
hi 
n/v A,dıe= Au 4 o Pfeilerdichte) . 
\ n Fir (@ 479 
Ent- 


Theorie der räumlichen Kerbwirkung. 
Von H. Neuber in München. 


n früheren Arbeiten wurden vom Verfasser grundlegende Probleme der ebenen und achs- 
I symmetrischen Kerbwirkung') behandelt. Durch Aufsuchen von Spannungsfunktionen in 
einem jeweils geeigneten krummlinigen Koordinatensystem konnte der elastische Zustand in 
verschiedenartig geformten Kerben bei verschiedenartiger Beanspruchung und für beliebige 
Kerbkrümmung rechnerisch erfaßt werden. 


ı) H. Neuber: Elastisch-strenge Lösungen zur Kerbwirkung bei Sch iben und Umdrehungskörpern, ZAMM, 
Bd. 13, 1933, S. 439. — H. Neuber: Zur Theorie der Kerbwirkung bei Biegung und Schub, Ing.-Arch., Bd. V, 1934, S. 238, 
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Ein vor kurzem vom Verfasser veröffentlichter neuer Ansatz zur Lösung räumlicher 
Klastizitätsprobleme’) ermöglicht es, diese Gedankengänge auch auf das für den Ingenieur 
so überaus wichtige Gebiet der räumlichen Spannungsverteilung in Kerben anzuwenden, 
welches bisher noch so gut wie unerforseht geblieben ist. 


Durch den neuen Ansatz lassen sich nämlich die Raumkomponenten des elastischen 
Verschiebungsvektors durch einfache Differentiation aus drei harmonischen Funktionen ableiten, 


IN 
“r 
0, 
IR 
| 
= Kons J 
0 
Abb, 1. Spannungsverteilung Abb. ?2. Die Höchstspannungen 
in einer Umdrehungsaußenkerbe bei Biegzung. in Umdrehungsaußenkerben bei Biegung. 


Schwbrichtung 


Schubrichtung 


Zımen gleicher Schubspannung 


Schnit u=0 Schnitt u=0 


Abb. 3 und 4. Spannungsverteilung in einer Umdrehungsaußenkerbe bei Schub, 


NG, wobei eine einfache Kombination dieser 
Funktionen die Rolle einer „räumlichen 

/W |  Spannungsfunktion“ spielt; die Theorie 
| so auf ein Spannungsfunktions- 
_ | verfahren zurückgeführt. Dabei ist 
2 von ausschlaggebender Bedeutung, daß für 

den Rechnungsgang auch in krummlinigen 


Koordinaten keine prinzipiellen Schwierig- 
FI] keiten bestehen. 
N | Bl Die Ableitung der einzelnen Lösungen 


wird an anderer Stelle ausführlich ver- 
öffentlicht werden. 
Die Abbildungen I bis 5 veranschau- 
ae lichen den Fall der Umdrehungsaußenkerbe 
Abb. 5. Die Höchstspannungen 
in Umdrehungsaußenkerben bei Schub. bei Biegung und Schub. 480 


®, H. Neuber: Ein nener Ansatz zur Lösung räumlicher Probleme der Elastizitätstheorie. Der Hoh'kegel unter 
Kinzellast als Beispiel, ZAMM, Bd. 14, 1954, S. 20, 
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Die Schwingungsvorgänge in Systemen mit zwei Freiheitsgraden. 
Von W. Quade in Karlsruhe. 


ie einfachsten Schwingungsprobleme, die bei Systemen mit einem Freiheitsgrad auftreten, 
führen, wenn man sich auf die freien Schwingungen beschränkt, auf eine lineare, 
homogene Differentialgeleichung 2-ter Ordnung mit konstanten Koeffizienten: 


dr 


oder (aD +bD+ . . 2.2.22. 


Dabei sind «a, b, e positive Konstanten, D ist der symbolische Differentialoperator, und x kann 
beispielsweise gedeutet werden als die von der Zeit f abhängige Ablenkung eines schwingenden 
\lassenpunktes aus der Ruhelage. 
Entscheidend für die Lösungen, welche (1) überhaupt haben kann, ist die Natur der 
Wurzeln der charakteristischen Gleichung von (1): 


Danach unterscheidet man folgende Fälle: 
Klasse [1 1]. Die Wurzeln von (2) sind voneinander verschieden. 
a) Die Wurzeln sind reell, aperiodischer Fall. 
b) Die Wurzeln sind konjugiert komplex, periodischer Fall. 
Klasse [2]. Die Wurzeln von (2) sind gleich, aperiodischer Grenzfall. 

Die hier gegebene Einteilung weicht ein wenig von der herkömmlichen ab. Die Zu- 
gehörigkeit zu einer Klasse wird dureh die Vielfachheit der Wurzeln von (2) entschieden. 
Innerhalb der einzelnen Klassen wird dann noch nach reellen und komplexen Wurzeln 
unterteilt. 

Damit sind alle in einem frei schwingenden System mit einem Freiheitsgrad möglichen 
Vorgänge aufgezählt. Die Unterscheidung der Vorgänge erfolgt nach formal-algebraischen 
Gesichtspunkten, nach den Eigenschaften der Wurzeln der algebraischen Gleichung (2). 

Das Problem der Aufzählung aller ın einem System mit mehreren Freiheitsgraden 
möglichen Vorgänge ist ungleich schwieriger. Zur Konstruktion von Fallunterscheidungen 
hat man sich in Ermangelung eines Besseren vielfach an das Nächstliegende gehalten, an die 
Art der Kopplung. Diese liefert aber kein Kriterium für den Ablauf des Vorganges. 

Der unten zur Klassifikation der Schwingungsvorgänge in Systemen mit 2 Freiheits- 
geraden eingeschlagene Weg beruht auf einer Verallgemeinerung der bei Systemen mit einem 
Freiheitsgrad schon von jeher gebräuchlichen formal-algebraischen Methode. Die verallge- 
meinerte Methode löst das Problem der Aufzählung der Schwingungsvorgänge für Systeme 
mit beliebig vielen Freiheitsgraden. 

Die freien Schwingungen eines Systems mit 2 Freiheitsgraden, dessen kleine Ablenkung 
aus der Ruhelage durch die Koordinaten .,, ., eharakterisiert wird, werden regiert durch 
ein System von linearen, homogenen Differentialgleichungen 2-ter Ordnung mit konstanten 
Koeffizienten: 

(a, 
(a, D+b,, D+e,)®, +(a,, D?+b,, D+ e,,).r, 
Die symmetrische Matrix des Systems (3) Ist 
D’+b,D+c,, | | 
D’+b, D+c„, 
und seine charakteristische Gleichung erhält man durch Nullsetzen der Determinante von M(D): 

Bei den Systemen mit einem Freiheitsgrad ist entscheidend für die 
Art des Vorgangs der Ausdruck aD®?+bD+e, bei den Systemen mit 
2 Freiheitsgraden ist es die Matrix M(D). 

Die Matrizes M(D) teilen wir zunächst in 3 Gruppen ein. 

I. Gruppe: Die Elemente von M(D) sind teilerfremd. 
Il. Gruppe: Der größte gemeinsame Teiler aller Elemente von M(D) ist ein Poly- 
nom vom 1-ten Grade in D: 


D s,)(a,D  0,), (D -s,)(a,,D 
= 
— s,)(a,, D—0,), (D—-s,)(u,, D 


MiD)— AD’+BD+C, 


Mi 


(4) hat hier immer eine reelle Doppelwurzel. 


x 
10 
18 
15 
105 
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III. Gruppe: Der größte gemeinsame Teiler aller Elemente von M(D) ist ein Polv- 
nom vom 2-ten Grade in D: 
s,)(D—s,), a,(D—s,)(D 
la, (D—-s,)(D- s,), a,(D—s)(D s,)| 
(4) hat hier immer 2 Doppelwurzeln. 

Innerhalb der einzelnen Gruppen unterteilen wir nach Klassen, wobei wir wie bei 
den Systemen mit einem Freiheitsgrad die Vielfachheit der Wurzeln von A(D) =0 ent- 
scheiden lassen. 

I. Gruppe: [1111], [211], [22], [31], [4]. 

Die in eckige Klammern gesetzten Ziffern sind die Exponenten der verschiedenen Linear- 
faktoren, in welche sieh .I(D) zerlegen läßt. 

II. Gruppe: [A [A [51]. 
Die Summe der in runde Klammern gesetzten Ziffern ergibt jeweils den Exponenten des 
Linearfaktors D-- s,. 

III. Gruppe: bb], [22]. 
Ks gibt hier nur 2 Klassen: s,$s, und s, 

Damit sind die Matrizes M(D) in 11 Klassen eingeteilt. 

Bei dieser Einteilung ist ein Umstand unberücksiehtigt geblieben, der für die bei 
Schwingungsproblemen auftretenden Matrizes M(D)=4AD’+BD-+C charakteristisch ist. 
A, B,C müssen symmetrisch sein und Matrizes von positiv definiten quadratischen Formen, 
nämlich von 27T, 2F, 27, wo T, F und U die kinetische Energie, die Rayleighsche 
Zerstreuungsfunktion und die potentielle Energie des Systems sind. Matrizes dieser Art 
eibt es aber für die Klassen | (11)2] und [(3 1)] nieht, eine eigentümliche Erscheinung. 

Bleibt noch zu untersuchen, ob und in welcher Weise innerhalb der einzelnen Klassen 
komplexe Wurzeln auftreten können. In der folgenden Tafel ist für jede Klasse (mit Aus- 
nahme der oben ausgeschlossenen [(1 1)2] und [®1)]) die Zahl der komplexen Wurzelpaare, 
die auftreten können, angegeben. 

I. Gruppe: [I 111], [211]. [22], [31], 


0 0) 0 0 0 
| | 
2 2 
I. Gruppe: [21]. 
0 0 


II. Gruppe: [22]. 
0 0) 


Bei einem System mit 2 Freiheitsgraden ist die Gesamtzahl der 
wesentlich voneinander verschiedenen Schwingungsvorgänge 14. 

Eine reinliche Scheidung zwischen periodischen und aperiodischen Vorgängen ist nicht 
möglich, denn in der Lösung des Systems (3) können gleichzeitig Funktionen mit periodischem 
und aperiodischem Verhalten auftreten. Man kann höchstens einteilen in aperiodische, periodische 
und gemischte Vorgänge und kommt dabei auf 9 aperiodische, 3 periodische und 2 gemischte. 

Die mathematischen Grundlagen dieser Untersuchungen sind einem speziellen Kapitel der 
linearen Algebra entnommen, der von Weierstraß begründeten Elementarteiler- 
theorie. Gewisse bei der Untersuchung der Matrix M(D) vorkommende Potenzen von 
Linearfaktoren nennt man KElementarteiler, und diese Teiler sind Invarianten, Invarianten 
nieht nur gegenüber linearen Transformationen mit konstanten Koeffizienten, sondern sogar 
gegenüber gewissen Transformationen, die Differentialquotienten der Koordinaten enthalten. 
Diese Invarianten charakterisieren diejenigen Eigenschaften eines schwingenden Systems, 
welehe von der zufälligen Wahl der Koordinaten unabhängig sind, insbesondere den Ablauf 
des Vorgangs, denn dieser kann ja nicht von der Wahl des zugrunde gelegten Koordinaten- 
systems abhängen. Die Art der Kopplung ist dagegen nicht invariant. Man kann weder 
von der Art der Kopplung auf den Ablauf des Vorganges, noch von der Klasse des Vor- 
ganges auf die Art der Kopplung schließen '). 481 


!) Vgl. des Verfassers Arbeit: „Klassifikation der Schwingungsvorgänge in gekoppelten Stromkreisen“ 
(Habilitationssehrift, Karlsruhe 1933). 
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Angenäherte Berechnung des kleinsten Eigenwertes zusammengesetzter 
Systeme. 


Von H. Schaefer in Hannover. 


ei der Behandlung von Eigenwertproblemen wird man häufig auf die Untersuchung der 

Abhängigkeit der Eigenwerte von Parametern geführt, die in das betreffende Problem 
eingehen. Die Methoden der Störungsrechnung versagen, wenn es sich darum handelt, 
solehe Abhängigkeiten im „Großen“ zu betrachten. Sind z. B. in der Differentialgleichung 


die Parameter A, und A, so festzulegen, daß die Lösungen gewissen Randbedingungen genügen, 
so liefern die jeweils zusammengehörigen, in einer (A,, -1,)- Ebene aufgetragenen Parameter- 
werte die „Eigenkurven“, deren Gestalt (Auftreten von Extremwerten und Wendepunkten) 
aus der Differentialgleichung heraus zu begründen noch nicht gelungen ist. Außer einigen 
vein formalen Ansätzen von van den Dungen') scheint über solche oder ähnliche Frage- 
stellungen keine Literatur zu existieren. Man kann jedoch auf Grund der Minimaleigenschaft 
des kleinsten Eigenwertes wenigstens über die Eigenkurven niedrigster Ordnung einiger, bei 
Eigenwertproblemen häufig auftretender Differentialgleichungen Aussagen von gewisser Trag- 
weite machen. Mit Hilfe der Greenschen Formel zeigt man leicht, daß die Gleichung der 
Tangente im Punkte A,=«a, JAJ,=P einer Eigenkurve von (1) (Abb. 1) bei solchen Rand- 
bedingungen, für die eine symmetrische Greensche Funktion existiert, die Gestalt 


1, 

A, A, | (2) 


0 0 


besitzt, wenn  () die zu diesen Parameterwerten gehörige Eigenfunktion ist. Aus der bekannten 
Minimaleigenschaft der rechten Seite von (2) folgt sofort, daß die Eigenkurve niedrigster 
Ordnung, deren Ordinaten die jeweils kleinsten zu 4, gehörigen Eigenwerte 4, sind, immer 
unterhalb ihrer Tangente verläuft (konvex gekrümmt nach positiver -1, - Riehtung), monoton 
fallend bei nichtnegativen y, und g.. Da der zweite Summand 'in (2) im allgemeinen nur 
schwach auf eine Änderung von w(xr) reagiert, kann man die Eigenkurve in Annäherung 


Aa 
4 
oo 
/ oo p 
-/----- A 
— 
A p 
Abb. 1. Abb. 2. 


(dureh eine Gerade ersetzen. Kennt man z. B. die kleinsten Eigenwerte 4, und /, zweier 
schwingender Systeme gleicher Steifigkeit p, mit den voneinander verschiedenen Massen- 
belegungen q, und q,. so folgt aus der für positive A, und A, gültigen Abschätzung 
A 
mit 1, = und A,= u,/ für den kleinsten Eigenwert eines dritten Systems 
+ 4, dessen Massenbelegung sich linear aus denen der beiden ersten 
zusammensetzt, die bekannte Dunkerlysche Formel 
2 
möglichkeit sich aber noch auf eine ganz andere Art von Zusammensetzungen erstreckt. 
Dureh Einführen der neuen abhängigen Veränderlichen py’ =z, die wir wieder als Aus- 


deren Anwendungs- 


F.H. vanden Dungen: Les equations integrales A plusieurs parametres et la teehnique des vibrations. 
Verhandlungen des 2. intern. Kongr. f. teehn. Mechanik, S. 113. 


= 
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lenkung deuten wollen, geht das System (py’’+7q49y=V über in das ihm „adjungierte*“ 


| :) +4 z=0, das dieselben Eigenwerte besitzt. (Man kann einen unsymmetrischen 

Kern angeben, dessen adjungierte Eigenfunktionen im Sinne der E. Schmidtschen Theorie 
die Eigenfunktionen dieser beiden Systeme sind.) Die für den kleinsten Eigenwert des 


l % 1 
zusammengesetzten Systems F 2’) + 7 —0 gültige Abschätzung . =. +. bleibt 


beim Übergang zum adjungierten System 
/ 
Pı TP> 


bestehen. Das einfache Beispiel der Torsionsschwingung einer homogenen Welle mit 
Kreisquerschnitt, am einen Ende frei, am anderen eingespannt, soll diese Zusammensetzung 
von Steifigkeiten beleuchten (Abb. 2). Nach Vorschrift von (3) denken wir uns dieses System 
aus den Teilsystemen A, und B, zusammengesetzt, die stückweise unendlichgroße Steifigkeit 
besitzen und den Systemen 4, und B, gleichwertig sind. Eine zweimalige Anwendung der 
Dunkerlysechen Formel gibt für den Eigenwert 4, der proportional dem Quadrat der kleinsten 


- 


Frequenz ist, den Wert AD =, [wahrer Wert: 1a). In Anlehnung an L. Föppl*) wollen 


wir solche Zusammensetzungen Hintereinanderschalten von Steifigkeiten nennen. Ganz ent- 
sprechend kann man, von der Differentialgleichung [(4,p, + 4A, p.)y’ ausgehen, 
Nebeneinanderschalten von Steifigkeiten und Hintereinanderschalten von Massen definieren. 
Hierbei wird man auf die Southwellsche Formel 424, +4, geführt. 

Selbstverständlich lassen sich solche Überlegungen sofort auf die Gleichung 4. Ordnung 
des schwingenden Balkens übertragen. 

Auf kompliziertere Zusammensetzungen wird man bei der Betrachtung der Differential- 
eleichung 


veführt. Erwähnt sei nur, daß man bei positiven p,r,g wieder durch Untersuchung der 
Eigenkurve niedrigster Ordnung im Gebiete positiver 1, und A, die einfache Abschätzung 
+ ; -1 gewinnt, wobei 4, und 4, wieder kleinste Eigenwerte der Teilsysteme sind, 
die aus (4) entstehen, wenn dort „1, bzw. 1, Null gesetzt wird. 

Im Bilde wurde an zwei Beispielen gezeigt, wie man einen sehr schmalen, langgestreckten 
Bereich abgrenzen kann, in dem die Eigenkurve verlaufen muß. 

Die hier angestellten Überlegungen führen auf eine neuartige Berechnung der oberen 
und unteren Grenze kleinster Eigenwerte, worüber an anderer Stelle ausführlich berichtet 
werden soll. 482 


Laminare Kanaleinlaufströmung. 
Von MH. Schlichting VDI in Göttingen. 
(Aus dem Kaiser Wilhelm-Institut für Strömungsforschung:.) 


m nach schon mehrfach mit Erfolg angewandten Methoden ') die Berechnung der 

Turbulenzentstehung für «die Strömung in einem Kanal mit parallelen Wänden durch- 
führen zu können, ist es erforderlich, die Geschwindigkeitsprofile der laminaren Anlauf- 
strömung im Kanal (ebenes Problem) mit einiger Genauigkeit zu kennen. Die früher von 
Prof. Schiller?) für die Berechnung der Anlaufströmung verwendete Methode geht anstatt 
von den Diflerentialgleichungen der Bewegung von dem Impulssatz aus. Dabei werden die 
Anlaufprofile durch zwei Parabelhälften an den Wänden und einen dazwischen liegenden 
Kern konstanter Geschwindigkeit angenähert. Den Widerstand im Anlauf, der ja ein Integral- 
wert über den ganzen Anlauf ist, gibt diese "Theorie in befriedigender Übereinstimmung mit 
den Versuchen; die Geschwindigkeitsprofile kann sie jedoch mit der für eine Stabilitäts- 
untersuchung erforderlichen Genauigkeit nicht liefern. Es wird deshalb im folgenden 
gezeigt, wie man diese direkt durch Integration der Bewegungsgleichungen erhalten kann. 


*, L,. Föppl: Über das Ausknieken von Gittermasten usw., ZAMM, Bd. 13, 19533, 8. 1. 

', W. Tollmien: Nachr. Ges. Wiss. Göttingen, Math.-Phys. Kl. 1929 und H. Schliehting: Die Naturwissen- 
sehaften 1934, S. 370. 

?, |. Sehiller: Forsehungsarbeiten auf dem Gebiet des Ingenieurwesens, Heft 248, 1922. 
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Eine Anlaufrechnung von G. Vogelpohl?), bei der die Differentialgleichung der Bewegung 
in eine Integralgleichung umgeformt wird, hat mit einer älteren Rechnung von Boussinesq') 
‚las gemeinsam, daß beide Entwicklungen von hinten her sind, d. h. von der Poiseuille- 
schen Parabel als erster Näherung ausgehen und die Abweichungen von dieser berechnen. 
Ganz vorn, wo die Profile sich erst wenig von dem rechteckigen Eintrittsprofil unterscheiden, 
treten dann Konvergenzschwierigkeiten auf. 

Für unsere Reehnung legen wir die «-Achse in die Kanalachse und die y-Achse in den 
Kintrittsquerse hnitt. Wir machen dann, um den ganzen Anlauf mit guter Genauigkeit zu be- 
kommen, zwei Entwicklungen: I. eine erste vom Einlauf «==-0 und von den Wänden aus, 
und II. eine zweite asymptotische Entwicklung für große x von der Kanalmitte aus. In einem 
(uerschnitt, wo beide Entwicklungen gelten, werden sie aneinander angeschlossen. 

I. Ganz vorn, unmittelbar nach Eintritt der rechteckigen Geschwindigkeitsverteilung U, 
in den Kanal, bildet sich an den Wänden die Reibungsschicht nahezu so aus, wie es von 
Blasius für die tangentiell angeströmte Platte in der unendlichen Flüssigkeit berechnet 
wurde. Weiter stromabwärts jedoch wird die Entwicklung der Reibungsschieht durch die 
Beschleunigung der Kernströmung 1’(“) immer stärker beeinflußt. 

Das Rechenverfahren ist das in der Prandtlschen Grenzschichttheorie°) übliche. Von 
den beiden Bewegungsgleichungen wird nur die der Hauptströmungsrichtung berücksichtigt: 


— kin. Zähigkeit), wobei die Kernströmung U(#) durch die Kontinuität 
Ula—ö*)=T,-a (U u)dy= Verdrängungsdieke; a = halbe Kanalbreite‘ 
gegeben ist, oder 
Nach Blasius®) ist für die nr angeströmte Platte in der unendlichen Flüssigkeit: 
- 
a 


Man erhält also eine Entwicklung der ENG nach Potenzen der für den Anlauf charak- 


teristischen dimensionslosen Größe & -V 


a?l 
(Ra). 
Führt man noch für die y-Koordinate, wie bei Blasius, einen neuen Maßstab ein durch 
„= (für negative y), so lassen sich Stromfunktion y und Geschwindigkeit « in den 


neuen Variablen e und ;, folgendermaßen entwickeln: 


Durch Einsetzen in die Differentialgleichung (1) erhält man dann, wenn man nach Potenzen 

Die Differentialgleichung für die erste Näherung £, ist die schon von Blasius gelöste 
mit den Randbedingungen = für „= und für große 7 (7 = 
Die Differentialgleichungen für die folgenden Näherungen sind sämtlich linear, inbomogen und 
von 3. Ordnung. Für die zweite Näherung erhält man 


mit den Randbedingungen: 


IK für 


>) G. Vogelpohl: ZAMM, Bd. 13, S. 466, 1933. 

#4, ). Boussinesq: Comptes Rendus, Bd. 113, S. 9 u. 49, 1801. 

5) 1. Prandtl: Verhdlg. d. III. Int. Math. Kongr., Heidelberg 1904. 
6, H. Blasius: Z. Math. Phys. Bd. 56, 1908. 


| 
24 


| IST I stco 1911 VITO 
\ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
R 0, 0 0 0 0, 0 0, 0, 0 0 0, 0 0, 0, 
= 
= 
[PT 
0 0 0 0 0 0 
N N n n 


L 


» 


1.44 


1.440 


1.904 


1.309 


1.090.) 


1.010 


— 


1,20 


1.10.09 
1.181 

1.193 
1.203 


.206 
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0.474 
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Eine ähnlich gebaute Differentialgleichung ergibt sich für die dritte Näherung £,, usf. Die 
Lösungen £,, £, ... lassen sich in der Umgebung von 7 =0 in eine Potenzreihe entwickeln, 
und für große n lassen sich asymptotische Lösungen angeben, ähnlich wie bei Blasius. 
Die vollständige Lösung für alle 7 ergibt sich daraus durch Zusammenfügung. Wir haben 
über das Blasiussche Glied hinaus zwei weitere Glieder £, und £, berechnet und das nächste [, 
abschätzend berücksichtigt ; damit kommt man etwa bis — — 0,016, wo mitte 1,239, 
0 0 
wo also der Anlauf etwa halb vollendet ist. 


Il. Für große Entfernungen vom Einlauf, wo das Profil sich schon allmählich der 
Poiseuilleschen Parabel nähert, machen wir eine nach #°, d.h. nach x asymptotische Ent- 
wicklung, indem wir setzen: 

or U af 
die Poiseuillesche Parabel ist und u“ eine Zusatzgeschwindigkeit, für die ein Ansatz 


% 


( a 
gemacht wird. 

Unter Berücksichtigung nur der größten Glieder erhält man für /, die lineare Differential- 
eleichung 4. Ordnung 


a? 


mit 4 als Eigenwert. Die Lösung dieser Differentialgleichung läßt sich nach Potenzen von yla 
entwickeln und ist durch die folgenden Randbedingungen eindeutig bestimmt: 


da €, noch frei verfügbar ist. 


Die Konstante €, wird dann noch dadurch bestimmt, daß für & = 0,016, wo sowohl die Ent- 
wieklung vom Einlauf her wie auch die von hinten her gilt, beide aneinander angeschlossen 
werden. Man findet ©, = — 0,3485 und A, = 18,75, und es zeigt sich, daß man mit nur einem 
Glied der Entwicklung (6) praktisch auskommt. 


N + + 4 + + | + 4 + + + 4 4 i $ + 
\ 
0oor 0002 0004 0006 0008 0000072 0.016 0020 002» 0,028 002 0» 00 006 008 05 
+ T + T - + + + + - ] + + 


I 1 
-02 + + + + 4 + + + 
| 


} 
% 

Abh. 1. Die Gesehwindigkeitsprofile der laminaren Anlaufströmung im Kanal 
(gestrichelte Kurve = Poiseuillesche Parabel). 


Die Resultate der Rechnung sind durch die nebenstehenden Abbildungen und Zahlen- 
tafeln wiedergegeben. In Abb. 1 und Tabelle I sind die Geschwindigkeitsprofile von 


=; 7, 0,0005 bis 0,016 durch die Entwicklung vom Einlauf her erhalten, die übrigen durch 


0 
die asymptotische Entwicklung für große x. Die ersten Profile der asymptotischen Rechnung 


zeigen noch deutlich die von der Reibung praktisch noch nicht erfaßte Kernströmung mit 
konstanter Geschwindigkeit in der Mitte. 


24° 


angew. 
"Mech. 
| 
| 
YA + r + + + | 
| | 
4 
| | | 
-710, | - | 
| 
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Abb. 2 zeigt den Verlauf der Geschwindigkeit in Ebenen parallel zur Wand „la = const. 


Die Geschwindigkeit in der Mitte steigt von [ I mit vertikaler Tangente zuerst sehr 
0 
vu - .. . 
schnell und dann langsamer an. Bei gr: 0,15 beträgt die Abweichung von der 


Asymptote nur noch rund 2°/o, hier ist der Anlauf also nahezu beendet. Nahe an der Mitte, 
ya=V1l,—=0,2, steigt die Geschwindigkeit monoton an bis zu dem asymptotischen Wert der 
Parabel: in größeren Abständen, bei ya =05 z. B., hat man zunächst Anstieg und später 
wieder Abfall. Das rührt daher, daß man sich bei kleinen Abständen vom Einlauf noch in 
der beschleunigten Kernströmung befindet und daher Geschwindigkeitsanstieg hat, während 
ınan ın größeren Abständen vom Einlauf in die durch die Reibung verzögerte Schicht gelangt, 
so daß die Geschwindigkeit wieder abnimmt. 

— Tabelle Il 


10 0 0,040 0,688 2,86 
I I | | | | 0001 | 01100 | 1833 | 0,050 | 0,772 | 2,57 
08 | | | | 0.002 | 0,1600 13,33 | 0,060 | 0,850 | 2,36 
u 7.0004 0.2322 968 0.080 1.002 | 2,085 
8.07 0,100 1.145 1,910 
0008 0338 | 705 | 0,185 18315 1,755 
77.0010 | 0,378 629 0.150 1.478 1,641 
0.012 0,41 571 0250 2101 1.400 
== 0016 | 0,466 485 0,500 3,601 1.200 
_ 0.082 0,620 3.23 
002 0.04 0.06 008 0,0 016 
on rz Abb. 2 (nebenstehend). Schnitte parallel zu den Wänden 
dureh dieGesehwindigkeitsprofile derKanalanlaufströmung. 


Ferner kann man noch den Druckverlauf ausrechnen. In Tabelle II ist der Druck- 
unterschied zwischen irgendeiner Stelle und dem Einlauf wo Ist, angegeben, 


und zwar dimensionslos Pe —i als Funktion von WU, Der Druck fällt zunächst sehr schnell 
/ 
Jo 
0 


ab, wie [ «=. Weiter hinten, wo der Anlauf beendet ist, hat man dann den linearen Poi- 
seuilleschen Abfall. Wäre von vorn an die ausgebildete Strömung vorhanden, so hätte 
Po. vr 


man einfach den linearen Abfall ! 6 : Der zusätzliche Druckverlust infolge der 
oO 
5 
Beschleunigung im Anlauf beträgt Ip: 0,601 5, 0°,? bezogen auf den Druck p, im Einlauf, 


wo u=ÜU,’). 
Aus dem Druekverlauf ist noch der Widerstandskoeffizient / ermittelt worden. Seı 4 
definiert dureh die Gleichung 


. 6 „a 
so ist für die ausgebildete Poiseuillesche Strömung |. 


Infolge des Anlaufs ist der Widerstandskoeffizient größer. Schreiben wir 
— . . . . . . . . . . . . . . (9), 


so Ist x ein von der dimensionslosen Anlauflänge #° = FE abhängiger dimensionsloser Faktor, 
der für große &° gegen I gehen muß. In Tabelle II ist x als Funktion von e’ angegeben für 


‘) Nach Mitteilung von Prof. Schiller ergibt seine Näherungsrechnung für den ebenen Fall Ip = — 0,614 ) Uo:. 


» 0 
also nur sehr geringe Abweichung von unserer genauen Reehnung (beim Kreisrohr Näherungswert Ip 1,16 5 Uo?): 


\ 
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0 
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den Fall, daß die eine Druckmeßstelle (,) im Eintrittsquersehnitt «= 0 liegt. Für &° = 1,0 weicht x 
noch etwa 10°, von seinem asymptotischen Wert ab. Liegt die eine Druckmeßstelle nicht 
zerade im Eintrittsquersehnitt, sondern an irgendeiner anderen Stelle im Anlanfgebiet, so läßt 
sich die entsprechende Widerstandszahl 4 mit den gemachten Angaben leicht errechnen. 

Schließlich kann man noch den Widerstand einer tangentiell laminar angeströmten 
ebenen Platte ausrechnen, die parallel zu den Kanalwänden in der Mitte eines Kanals der 
Breite 2a steht. Man findet für den Widerstandskoeflizienten: 


13284, , 696123 ur \ 


wo A, die mit der Plattenlänge ! gebildete Reynoldssche Zahl : bedeutet [e,.= 
UF 
= Widerstand, F= benetzte Fläche|. Das erste Glied ist die bekannte Blasiussche 


Formel, die folgenden geben die Widerstandsvermehrung durch die endliche Kanalbreite. Sie 
können unter Umständen von beträchtlicher Größe sein. 

Die hier für das ebene’Problem durchgeführte Rechnung läßt sich ohne weiteres auf 
das rotationssymmetrische Problem der Rohreinlaufströmung übertragen. 483 


Umkehrung von Integraltransformationen. 
Von Günther Schulz in Berlin. 
ei der mathematischen Behandlung physikalischer und technischer Fragen wird man 
B öfters auf Integraltransformationen der Gestalt 


+% 


geführt. Meist handelt es sich um die Umkehrung dieser Transformationen, bei gegebenem 
Kern K (x) und gegebener rechter Seite f(#) die Funktion y(#) zu bestimmen. Solche Integral- 
gleichungen erster Art sind bereits mehrfach behandelt worden, ohne daß die Ergebnisse für 
die praktische Auflösung nutzbar gemacht worden sind. 

Vorausgesetzt wird, daß K(#) im Unendiichen genügend stark verschwindet. Fast 
immer ist A (x) nichtnegativ und in den Veränderlichen x und E symmetrisch; es ist also 
K(x)= K(— Ferner sei 


Die Transformation ist also eine Mittelbildung mit der Gewichtsfunktion K(*). Die 
Funktionen y(#) und f(x) sollen periodisch mit der Periode 2 7 sein. 

Den im folgenden angegebenen Beispielen ist gemeinsam, daß das Kurvenbild (x) einer 
physikalischen Größe durch örtliche oder zeitliche Nachbareinflüsse „verundeutlicht“ wird, 
so daß bei der Kurve f(“#) die Maxima von (x) erniedrigt und die Minima erhöht sind, und 
zwar um so weniger, je ausgeprägter das Maximum von K(x#) an der Stelle «0 ist. 


Beispiele: 


Schwärzung photographischer Schichten K(x)=ce-*i* 
Il. Optische Abbildung eines selbstleuchtenden Objekts mit sinz.’ 
.. . . N (.r) P 
Berücksichtigung der Beugzung') 


III. Wirkung des Dopplereffekts auf die Intensitätsverteilung 


im Spektrum?) K (x) = ce-»?* 


IV. Verundeutlichung des Spektrums infolge der endlichen K(x)=0, für -o...—u 
Breite des Spektrometerspalts, 
Verfälschung der aufgenommenen Intensitätskurve im + für —a...+a 
Spektrum infolge der endlichen Breite des Bolometer- ns 
streifens’?) —-0, für +a...+®. 


ı) Vgl. L. Mandelstam: Festschrift Heinrich Weber gewidmet, Leipzig und Berlin 1912, S. 228 bis 241. 

?) Dieses Problem ist bereits von G. Doetseh ersehöpfend behandelt worden: Ztsehr, f. Phys. 49, 1928, S. 705 bis 730. 

3) Das Problem IV ist von €. Runge behandelt worden: Ztschr. f. Math. u. Phys. 42, 1897, S. 205 bis 213. In 
der hier verwendeten Ausdrucksweise wird dort die Integralgleichung (}) mit dem zum Kern (IV) zehörigen zweiten 
iterierten Keru betrachtet und ihre Lösung in einer nach geradzahligen Ableitungen von f fortschreitenden Reihe an- 
zereben. Dieses Verfahren hat den Nachteil, daß es sieh schwer auf andere Kerne übertragen läßt und auch über die 
Kindeutigkeit der Lösung keine Auskunft gibt. 
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Auch viele Probleme des Wärmeausgleichs sowie Nachwirkungserscheinungen an 
Registrierinstrumenten gehören hierher. 
Die Behandlung der Gl. (1) wird übersichtlicher, wenn an Stelle von K(x) ein 
periodisierter Kern eingeführt wird, der durch 
v=— 
definiert ist. Nach einer weiteren leichten Umformung ergibt sich 


ra 


Alle drei Funktionen K(z), y(#), f(*) sind „zyklisch“, hängen also nur von der Differenz 
zweier Veränderlichen ab und sind in dieser Differenz periodisch mit der Periode 2x. 
Man sieht, daß es im wesentlichen dasselbe Problem ist, K aus f und y wie y aus K und f 
zu bestimmen. Wirken zwei Einflüsse, gekennzeichnet durch die Kerne K, und K,, dann 


ist der durch Komposition entstehende Kern - 
....Ö) 


einzuführen, der ebenfalls zyklisch ist. Gl. (5) hat übrigens denselben Bau wie Gl. (4). 

Bei der praktischen Durchführung der Rechnung benutzt man entweder die Werte der 
vorkommenden Funktionen an n äquidistanten Stellen im Intervall —=...-+7 oder deren 
Fourierentwieklungen. Im ersten Fall wird man auf die Gleichung KY)= 5 geführt, in 
der sämtliche Matrizen zyklisch sind. (Eine »-reihige Matrix heißt zyklisch, wenn die Ele- 
mente nur von der Differenz der Indizes abhängen und wenn die Elemente, deren Indizes 
kongruent mod n sind, gleich sind.) Systeme linearer Gleichungen mit zyklischer Matrix 
sind bequem auflösbar und deren Eigenwerte leicht angebbar. Die durch XX,=x,X%, 
definierten Eigenmatrizen einer symmetrischen, zyklischen Matrix St sind von St unabhängig. 
Dem entspricht, daß ein symmetrischer, zyklischer Kern stets cosyx= und sinv«= zu Eigen- 
funktionen hat. 

Bei Benutzung der Fourierentwieklungen gilt für die Koeffizienten bei passender 
Bestimmung der Konstanten 


Dabei sind die Kosinuskoeffizienten von K y(#), f(x) und die Sinus- 
koeffizienten. Das trigonometrische Polynom »-ten Grades mit den so bestimmten Koeffizienten 
gibt die beste Approximation im Mittel. Mit „(#) ist auch y(a@)—+h(x) eine Lösung von 
Gl. (1), wenn h(x) eine Lösung der zu Gl. (1) gehörigen homogenen Gleichung mit ver- 
schwindender rechter Seite ist. Da dieser die Beziehungen 


zwischen den Fourierkoeffizienten von K(x) und h(.) entsprechen, so ist die Lösung von 
Gl. (1) immer dann vieldeutig, wenn in der Fourierentwicklung von A (x) verschwindende 
Koeffizienten vorkommen. 

Der Fall aperiodischer Funktionen kann meistens auf den hier behandelten Fall der 
Periodizität zurückgeführt werden, da die in den Anwendungen vorkommenden Funktionen 
außerhalb eines genügend großen Intervalles gleich Null angenommen werden dürfen. 

Eine ausführliche Darstellung mit durchgerechneten Beispielen soll erscheinen. 484 


Kavitationsversuche am Massachusetts Institute of Technology und ihre 
Deutung. 


Von E. W. Spannhake in Karlsruhe. 


s wird besonders die in Cambridge gründlich studierte Periodizität der Kavitation hervor- 
F gehoben. Ein Film wird vorgeführt, in dem gezeigt wird, daß das Totalvolumen der 
Kavitation periodisch als Ganzes gebildet wird und als Ganzes zusammenstürzt. Entgegen 
(der bisher üblichen Auffassung der Bildung, des Stromabtransportes und darauffolgenden 
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/usammensturzes von Einzelblasen wird die Vermutung ausgesprochen, daß der hier beob- 
achtete Vorgang ein für die Kavitation genereller sei. Es wird gezeigt, daß schlagartige 
Erschütterungen von der Periode der zusammenstürzenden Totalvolumina elektrisch gemessen 
werden können, und daß der Materialschaden tatsächlich nur an den Stellen auftritt, wo 
dieser Zusammensturz im Film beobachtet werden konnte. Versuche haben ergeben, daß 
(lie folgende Abhängigkeit der charakteristischen Größen der Kavitation von der Geschwindigkeit 
besteht. (Natürlich vorausgesetzt, daß der Kavitationszustand, d.h. der Ort des Zusammen- 
sturzes derselbe bleibt.) 


1. Der Energieverlust, verursacht durch die Kavitation, ist proportional dem Quadrat 
der Geschwindigkeit. 

2. Die Frequenz der zusammenstürzenden Totalvolumina ist der Geschwindigkeit direkt 
proportional. 

3. Verkleinert oder vergrößert man den Strömungsraum im Gebiete der Kavitation 
ähnlich, so ist die Frequenz umgekehrt proportional einer Länge. 


Ferner wurde festgestellt: 


Die Gestaltung des vor und hinter der Versuchsstrecke liegenden Rohrsystems, seine 
elastischen und trägen Konstanten haben keinen Einfluß auf die Frequenz. 


Aus allen diesen Gründen wird die Ansicht ausgesprochen, daß die Ursache des 
periodischen Charakters eine hydrodynamische und keine thermodynamische sein müsse, 
da ja die Gesetze der mechanischen Ähnlichkeit als gültig festgestellt wurden, und daß 
diese Ursache allein im Kavitationsgebiet selber, nicht aber in dem davor- oder dahinterliegenden 
Rohrsystem liegt. Als tatsächliche Ursache wird der Jlabile Charakter der bei Dampf- 
bildung vom Profil abgelösten Strahles vermutet, der bei geringer Druckerhöhung an der 
Ursprungsstelle sofort für eine gewisse kleine Zeit wieder zum Anliegen kommt und so eine 
Art „plunger* stromabwärts jagt, der dann das alte Kavitationsvolumen zerquetscht. Als 
Grund dieser periodischen Druckerhöhung wird im beobachteten Fall der periodische 
Zusammensturz des alten Volumens angegeben, im allgemeinen Fall aber vermutet, daß jede 
starke Verzögerung des „plungers*“ rückwirkend auf die Ursprungsstelle der Kavitation als 
momentane Druckerhöhung wirken kann. Es wird gezeigt, daß, während der „plunger“ strom- 
abwärts läuft, neues Wasser aus dem Hauptstrahl in ihn hineinströmt und so seine Masse 
ständig vergrößert. Der Materialschaden wird in Verbindung gebracht mit der verfügbaren 
kinetischen Energie des „plungers“ und mit der Art und Weise, wie er auf die kompakte 
Wassermasse stromabwärts unter Zerquetschung des alten Kavitationsvolumens auftrifft. Die 
Möglichkeit wird erörtert, den Einfluß der Profilformen auf die Ausbildung des „plungers“ 
systematisch zu studieren und auf diese Weise vielleicht Konstruktionsvorschriften zur Ver- 
meidung besonders gefährlicher Formen zu gewinnen. 486 


Ein Kriterium der Instabilität von laminaren 
Geschwindigkeitsverteilungen. 


Von W. Tollmien in Göttingen. 


cehon in den Anfangszeiten der Turbulenzforscehung ist die Entstehung der Turbulenz auf 
S eine Instabilität der Laminarströmung zurückgeführt worden. Doch blieb dieser Gesichts- 
punkt lange Zeit nur Programm. Erst in den letzten 5 Jahren gelang es, die Stabilitätsgrenze 
von einigen Laminarströmungen zu berechnen. So erfreulich das Ergebnis dieser Unter- 
suchungen ist, so muß man doch sagen, daß sie sich durchweg auf sehr spezielle Geschwindigkeits- 
profile beziehen. Es besteht in diesem Gebiet ein fühlbarer Mangel an allgemeinen 'Theoremen, 
die eine rasche Klassifikation von Geschwindigkeitsprofilen hinsichtlich ihrer Stabilität ge- 
statten. Im folgenden soll nun ein solches sehr allgemeines Instabilitätskriterium aufgestellt 
werden. 


Wir betrachten zweidimensionale Strömungen, in denen überdies die Geschwindigkeit 
im wesentlichen nur von der Koordinate y quer zur Strömung abhängen soll. Als ersten 
Typus ziehen wir Profile in Betracht, die etwa in schwach divergenten Kanälen auftreten 
können und symmetrisch zur Mitte sind (Abb. 1, U,,., — maximale Geschwindigkeit, b = halbe 
Kanalbreite). Wir schließen aber Profile mit Rückströmung oder Ablösung aus. Im übrigen 
sind die Profilformen wenig beschränkt. Vor allem kann in jeder Profilhälfte gerade ein 
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Wendepunkt auftreten. Als zweiten praktisch besonders wichtigen Typus betrachten wir 
Grenzsehichtprofile (Abb. 2, 9 = Grenzschichtdicke). 


Diesen Strömungen werden Störungen überlagert, die gemäß der Methode der kleinen 
Schwingungen als Wellen, die sieh in der Strömungsriehtung « fortpflanzen, angenommen 
werden. Auch die Störungen sollen zweidimensional sein, was nach einer neueren Unter- 
suchung von H. B. Squire keine Einschränkung bedeutet. Die Störungsbewegung läßt sich 
dann durch eine Stromfunktion beschreiben. Wegen der vorausgesetzten Kleinheit kann die 
Störung aus Partialstörungen aufgebaut werden. Die Stromfunktion einer Partialstörung 

ist in komplexer Schreibweise a(reell) = Wellenlänge; (komplex) 


— Pr zeitliche Kreisfrequenz, — lo- 
garithmisches Dekrement, 5; >0 bei Anfachung; 


Cr+ie;,c, Wellengeschwindigkeit. 

In erster asymptotischer Näherung für sehr 
eroße Reynoldssche Zahlen kann man, wie im 
einzelnen gezeigt wird, das explizite Auftreten 
der inneren Reibung unterdrücken. Freilich kann 
man dann nicht mehr eine kritische Reynolds- 
sche Zahl berechnen. Dagegen erhält man gute 
Auskunft über die Existenz und die Eigen- 
schaften der gefährlichen angefachten Störungen. 


Mathematisch formuliert lautet das Problem 
so. U(y) stelle die Grundströmung dar. Dann 
ist die Differentialgleiehung 


u” 

- a = = 

P- 7% 

Umax 
b b 

Abb. 1. Abb. 2. 


mit den Randbedingungen 9=0 an beiden Rändern (Wände oder das Unendliche) zu 
lösen. Dies Randwertproblem ist ungewöhnlich. Die reelle Größe a ist als gegeben anzu- 
sehen. e wird aus einem nichtlinearen Randwertproblem bestimmt. Nicht zu jedem a exi- 
stiert eine Lösung, obwohl e komplex sein kann. 


Wenn man von der vertieften Begründung des Problems durch asymptotische Betrach- 
tungen für sehr große Reynoldssche Zahlen absieht, so ist das eben formulierte Problem 
keineswegs neu. Bereits im Jahre 1550 veröffentlichte Lord Rayleigh eine Untersuchung, 
die ein wichtiges Resultat lieferte. Lord Rayleigh zeigte, daß Profle ohne Wendepunkt 
im oben formulierten Sinne stabil sind. Obwohl Lord Ravleigh selbst immer wieder (bis 
ins Jahr 1913) dies Problem angriff und nach ihm andere Forscher in verschiedenen Ländern 
diese Bemühungen fortsetzten, gelang dennoch keine Aussage über Profile mit Wendepunkt. 
Für diese stellen wir nun den einfach auszusprechenden Satz auf, daß sie sämtlich instabil 
sind. Dieser Effekt enthält die innere Reibung nicht explizite und ist daher bedeutend stärker 
als die Anfachung, die bei manchen Profilen ohne Wendepunkt in zweiter asymptotischer 
Näherung durch eine kleine Reibungskorrektur erhalten wurde. Zweifellos spielt dieser 
Effekt bei der Entstehung der Turbulenz eine fundamentale Rolle. 


Über dies erstaunlich umfassende Theorem hinaus werden im Verlauf des Existenzbe- 
weises in gewissem Umfange auch allgemeine Formeln für die Größe der Anfachung auf- 
gestellt. 485 
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KLEINE MITTEILUNGEN 


Mathematische Theorie des Schwingen- 
fluges 11. Den Ausführungen im ersten Teil 
meiner Arbeit über den Schwingentlug (vgl. diese 
Zeitschrift Bd. 14 8. 163 bis 172) sei unter Bei- 
behaltung der Bezeichnungen dieser Arbeit noch 
folgende kurze Ergänzung hinzugefügt: 

Die wichtigste Frage ist die nach dem Vortrieb, 
ler dureh den schwingenden Flügel erzeugt wird, 
weil von ihm in erster Linie der Erfolg der ganzen 
Vorrichtung abhängt. Es soll daher im folgenden 
eine Näherungsformel für den reinen mittleren 
Vortrieb V, also den Vortrieb ohne Widerstands- 
bestandteile, hergeleitet werden, die den Vorzug 
außerordentliceher Kürze und Einfachheit besitzt: 


Unter gewissen, unten näher angegebenen Voraus- 
setzungen gilt für den reinen mittleren Vortrieb /’ 
näherungsweise: 


( 


Zum Beweise formen wir die a.a. 0. angegebene 


Ausgangsformel für den mittleren Vortrieb T (der 
den Widerstandsanteil noch mit umfaßt): 

9m 


wie folgt um: Wegen 


(2) 
—= svto®*- 2 
Ist 
v 
2 

00 


a* 
Machen wir nun die Annahme, daß Er klein 
c 


zegen © ist, was bei nicht zu schnellen Schwingungen 
einigermaßen zutreffen dürfte. so kann im Integral 
näherungsweise I’ durch zeta* ersetzt werden, 
so daß man erhält: 


Ir 
a)sinpgdydr 
00 
In 
a)sinpydydı. 
An J. 


In analoger Weise kann der Ausdruck für den 
mittleren Auftrieb vereinfacht werden und liefert: 


= 
| [F'sin ...60). 
47T 


Wir stellen nun die analoge Frage, die in der 
gewöhnlichen Tragflächentheorie die elliptische 
Auftriebsverteilung liefert. nämlich die Frage 
nach dem maximalen mittleren Vortrieb 
bei fest vorgegebenem mittlerem Auf- 
trieb. 


Wegen (3) ist, wenn !=2bv- I A,(r) sinng: 
1 


IT In 
v 
0 


und somit wegen 


SL 
hveosvr cosq sinng 
= \ 4.(t) 
sin 


—— \ A,„sin 


00 
In 
sinng\ 
\ n An — sing 
n— | 
v 
177 
| hvcosvr \ '—1) * A, 
3 \ n Ay; A) 
Damit dies ein Maximum werden kann, muß 
er 2hvcosvtr, 
(4) 
2 hvcosvr-(- 1) 
Au= 
—4) 


sein. Danach schwingt also die Zirkulation um 
die dureh den Mittelwert von 4, gegebene elliptische 
Verteilung herum, was außerordentlich plausibel 
erscheint. Der Vortrieb selbst wird, wenn man je- 
weils den konstanten Bestandteil wegläßt, der den 
üblichen induzierten Widerstand liefert. 


657 657 \ 
0,797 
= | \ n 4? \ n Ay?) dt 
u 1 
Ir 


v 
ovv?’b? 
| \ n A, dr 
4 {nr 
6 


v 
b? [ 
12 92 2 2 


1 


4)? 


1 

Die Berechnung der unendlichen Summe liefert ge- 

g 50 daß schließlich (1) 


entsteht, wie behauptet wurde. 

Es sei ausdrücklich bemerkt. daß dieser Näherungs- 
wert, der i. a. größere Werte liefert als die a.a. 0. 
angegebenen Formeln, nur dann erreicht werden 
kann, wenn 


klein gegen ist, also die Schwingungen 


nügend genau den Wert 


nicht zu rasch sind, 

2. Flügeltiefe und Anstellwinkel zemäß 
Formel (2) in ihrer räumlich zeitlichen Ver- 
änderlichkeit dem durch (4) gegebenen Gesetz 
für die Zirkulation angepaßt werden, was aber 
möglich erscheint. 

Bemerkenswert an der Formel (1) ist vor allem 

die Unabhängigkeit von v® und «a. 
Breslau. Werner Schmeidler. 464 
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Berechnung der Eigenschwingungs- 
zahlen von Wellenleitungen durch Tei- 
iung der Feder in den Knotenpunkten. 
Um für die Berechnung der Eigenfrequenzen von 
Wellenleitungen eine anschauliche Rechnungs- 
methode zu geben, wurde von H. Wydler') das 
zegebene n-Massensystem durch Teilung der Massen 
auf n-1 Zweimassensysteme zurückgeführt, welche 
alle mit der gleichen Eizenschwingungszahl schwin- 
een. Eine andere anschauliche Methode gibt 
0. Föpp!?) an, welcher außer den Massen noch 
die Federlängen teilt und dadurch das »-System 
durch 2n-2 Einmassensysteme ersetzt. 

Vernachlässiet man den Einfluß der Wellenmasse, 
so genügt eine Teilung der Federlängen. um das 
zegebene System auf n-Einmassensysteme zurück - 
zuführen. Das Schwingungsbild sei gegeben; dann 
ist auch die Lage der Knotenpunkte bekannt, welche 
man als Einspannstellen auffassen kann. Das 
System zerfällt in „-Einmassensysteme nach Abb. 1, 


Ci 


N 

Mi; 

\A457.1] 

Abh. 1. 


die alle mit der gleichen Eigenschwingungszahl © 
schwingen. Für die Anordnung nach Abb. 1 gilt 
die Differentialgleichung 


[23 


Darin sind m; die öte Masse, bzw. c'i, i+1 
die Federkonstanten bezogen auf die Einspann- 
stellen (Knotenpunkte). Dann gilt für die KEigen- 
schwingungszahl 


ci — 
m; 


/,u bemerken ist nur noch, daß die Federkonstanten 
sowohl positiv als auch negativ sein können (vgl. 
Rausceh?®)). je nachdem ob der Schwingungsknoten 
in bezur zur Masse auf seiten der Feder oder 
auf der entzerengesetzten Seite zu liegen kommt. 
Sind Drehschwinzungen zu berechnen und bezieht 
man sich auf ein konstantes /» (Trägheitsmoment des 
Querschnittes), dann ist, wenn Pi,i+1 die 


zu c'i,i+1 zugehörigen Federlängen sind. 
Y 
 IMi-ni’ 


1 1 
und die Knotenpunkte sind nun so zu bestimmen, 
daß 
1 | 1 
5 . . (1) 
m; 
wird. Es ist somit ein Probierverfahren gegeben, 
um die Eigenschwingungszahl zu erhalten. Man 
braucht nur die Schwingungsform und damit die 
Knotenpunkte so lange zu ändern, bis die Konstante 
in Gl.(1) für alle Massen die gleiche geworden ist. 
Wien. Fritz Söchting. 457 


MH. Wydler: Drehschwingungen in Kolbenmaschinen 
und die Gesetze ihres Ausgleiches. Berlin 1922. 

", 0. Föppl: Dreh- und geradlinige Schwingungen von 
Wellen und Massen, ZAMM 1921, S. 367. 

>) E. Rauseh: Graphisehes Verfahren zur Bestimmung 
der Eigenfremenzen bei mehrgliedriger Schwingungs- 
anordnung, Ingenieurarechiv 1930, S. 20%. 


Über die Form eines Pfeilers gleicher 
Festigkeit bei gleichzeitiger Druck- und 
Biegebeanspruchung. Die folgende Unter- 
suchung beschäftigt sich mit der Frage, wie bei 
einem Pfeiler, der außer durch sein Eigengewicht 
noch durch eine an seinem oberen Ende angreifende 
Kraft belastet ist, die Breite sich mit der Höhe 
ändern muß, wenn die maximale Druckspannung 
in allen Querschnitten die gleiche sein soll („Pfeiler 
gleicher Festigkeit”). Der Pfeiler soll dabei recht- 
eckigen Querschnitt haben: die Ebene durch die 
Pfeilerachse und die Einzelkraft am oberen Ende 
gehe durch die eine Symmetrieachse des Rechtecks, 
senkrecht zu dieser Ebene habe das Rechteck eine 
konstante Tiefe, die wir ohne Einschränkung der All- 
eemeinheit eleich eins annehmen wollen. 

1. Am oberen Ende wirke eine horizontale 
Kraft P. Ist 4 die gesuchte Breite des Quersehnitts 
und y das spezifische Gewicht des Pfeilermaterials, 
so ist 

x 


0 


das Eigengewicht, das auf einem Schnitt in der 
Entfernung x unter dem oberen Ende lastet (siehe 


Abb. 1). Dasselbe erzeugt eine Druckspannung ; 


p 


24321 
Abb. 1. 


zu dieser kommt die Druckspannung infolge der 
8oll die 
maximale Druckspannung R für alle x dieselbe 
sein, so mub 


Biegungswirkung der Einzelkraft 


werden. Multiplizieren wir (2) mit 4 und differen- 
zieren, so erhalten wir 


6Prdy 6Pdre 
yydıe— Rdy . (9 
bzw. wenn zur Abkürzung 
a 


gesetzt wird. 


Für diese Differentialgleiehung findet man einen 
nur von 4 abhängigen integrierenden Faktor 


1 


(Y 1 a? y? 


mit dessen Hilfe man das Integral von (5) zu 


y? b 1 | 
a? =c..®W 


Y a 


| 
| 
aEZZZA 
| 
| 


.angew. 
d Mech. 


icher 
und 
Unter- 
vie bei 
'ewicht 
:ifende 
Höhe 
ınnung 
Pfeiler 
recht- 
ch die 
ı Ende 
ıtecks, 
k eine 
er Äll- 


ontale 


"hnitts 
erials, 


(1) 


n der 
(siehe 
() 


der 
die 


selbe 


(2) 


rell- 


nen 


(6). 


(7) 


‚and 14, Heft # 
Dezember 1934 


Kleine Mitteilungen 379 


rmittelt. Die Integrationskonstante ec ermittelt 
san aus der Bedingung, daß für #0 nach Gl. (2) 


sit Q=0 sowohl y als verschwinden müssen. 


Das ergibt 


wird 


+a?y 


IntaytVi-tayy. . (9. 


Abb. 2. 


Man trägt der in der Biegungstheorie vernach- 
lässigten Schubspannung Rechnung, indem man die 
Breite am oberen Ende nicht auf Null abnehmen 
läßt, sondern dem obersten Stück eine durch die 
zulässige Schubspannung ZA, bestimmte konstante 
Breite £ gibt. Abb. 2 stellt das Ergebnis für ein 
numerisches Beispiel dar, welchem die Zahlenwerte 


„2 R=1W 
m m“ 
. (10) 
to 25 
pP = 3 to 
p 


8 


_ 
ar 


Abb. 3. 


zugrunde liegen. Die Gewichtsersparnis gegenüber 
einem prismatischen Körper gleicher Breite ist be- 
deutend. Vergleicht man in Abb. 2 den oberhalb 
des Schnittes m — n liegenden Teil mit einem 
Körper konstanter Breite, so ist bei gleicher maxi- 
maler Druckspannung das Gewicht des Körpers 
rleicher Festigkeit um 46% geringer als das Ge- 
wicht des prismatischen Körpers. 


2. Am oberen Ende wirke eine vertikale Kraft P 
in der Entfernung ce von der Mitte (Abb. 3). Die 
größte Druckspannung wird 


Y 
Nach Multiplikation mit y und Differentiation er- 
gibt sich 


6Pedy 


was sich ohne weiteres integrieren läßt: die auf- 
tretende Integrationskonstante bestimmt sich aus 
der Bedingung, daß am oberen Ende 2=0, wo 


Q=0 ist. nach Gl. (11) die Breite £ durch die 
Gleichung 


r ...x.:... (3 
bestimmt ist. Mit dem aus (13) zu ermittelnden 
Werte von £ erhält man 

Kiew. W.Swida. 49 


(14). 


Graphische, näherungsweise Bestim- 
mung von Kippziffern mittels der Seil- 
kurve. Für viele praktisch vorkommende Kipp- 
fälle ist es noch unmöglich, die genaue Kipp- 
last anzugeben. So z.B. für die Belastung mit 
mehreren Einzelkräften, Einzelmomenten und für 
sprungweise veränderliche Steifigkeiten. In den 
meisten Fällen ist die Differentialgleichung des 
Problems gegeben. Ist diese Differentialgleichung 
für den Kippwinkel $ von zweiter Ordnung und 
linear, dann läßt sieh — wenn die Randbedingungen 
durch ® allein ausdrückbar sind — aus einer 
Näherung P,, für den wahren Verlauf von ® mittels 
der Seilkurve eine verschärfte Näherung P,„;, g2e- 
winnen. Aus 2, und ®,,, wird hier dureh .Ver- 
eleich der Ordinatenflächen“ ein Näherungswert 
für die Kippziffer gefunden |1]. — Dies ist auch 
möglich, wenn eine die Verdrehung einzelner Quer- 
schnitte hindernde Stützungz nach Abb. 1 zu be- 
rücksiehtigen ist (Flugzeugholm und Rippe). 


Abh. 1. 


Das an sich bekannte Verfahren soll für den 
Fall des Kragträgers mit zentrisch angreifender 
Endlast kurz erläutert werden. Die Differential- 
gleichung lautet |2]: 


d2® pen 


(z ist die dimensionslose Längskoordinate, ge- 
zählt vom Angriffsquerschnitte der Endlast P, ! be- 
deutet die Trägerlänge. A die kleinste Biegungs- 
steifigkeit und € die Torsionssteifigkeit des kon- 
stanten Trägerquerschnittes). — Die Randbedin- 
gungen für das Kippen des Kragträgers sind 


dx 


P(r=1)=0 (2a) und 
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S| nahme 2, entsprechenden Ordinaten 9, zu erhalten 
7 S; U und damit weiters in gehöriger Größe gezeichnet 
A werden kann, wurde zweckentsprechend mit 10 facher 
2 bei 1:1.7facher Poldistanz gearbeitet. Der 
1? erste Näherungswert für die Kippziffer ergibt sieh 
S 
| | 
| d.r (4) 
| | | S In+1 « 
I IE 1 NS RS v 
=) | | yo 
| Sal | Vielfaches der Lastvergrößerung 
| — u Beim zweiten Schritt wird der Verlauf von », zur 
Berechnung der Belastungsintensität verwendet usf. 
Abb. 2. In der weiter unten angegebenen Tabelle ist 
” für jeden der betrachteten 6 Fälle das Belastungs- 
Die gewählte eı BR schema. die Differentialgleichung und die nach der 
(3) besprochenen Methode ermittelten Näherungswerte 


erfüllt schon die Randbedingungen (2a) und (2b). 
Als Einheit der Zeichnung (siehe Abb. 2) wurde 
100) mm im Original angenommen. Die Ordinaten 
unter der zgestrichelt gezeichneten Kurve ent- 
sprechen einer „Belastungsintensität” 7° 2,, welche 
mit einer Poldistanz / =1 eine im Verhältnis 1: X? 
kleinere Verdrehung ®,' ergeben würde. Um die 


für die Kippziffer angegeben. Als erste Annahme 
für ?P wurde immer die quadratische Parabel nach 
61. (3) angenommen. Es gelten jeweils die Rand- 
bedingungen (2a) und (2b). nur im Falle e) (siehe 
die Tabelle) ist der von der Einzelstützung gegen 
Verdrehen an der Stelle #0 ausgeübte Drall in 
der Form 


Ordinaten des Seilzuges 7, der Zeiehnung in der (7 1 B(r—0). (?e) 
eleiehen Größenordnung wie die der ersten An- dr 
Tabelle. 
= 
BER I SER | 
Belastungsschema Differentialgleichung | | | 3° 
E - 
- 1% 4.15 3.4 
| 
2 7A 
| — 2 3.54 
— pen 3,43 
L 
1 No d .r? ( | 5.20 [3 1.5 
dN= \ 
m 1,35 
2—udu AU 3.54 3,36 | 5,35 
u\. 3.41 [5] 1.49 
f wenn v=1—r.B=d, 4’ ist 
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u berücksichtigen. In den drei letzten Fällen ist 
uch die dort notwendige erste Annahme für X 
anzereben. weiters ist in den Fällen a), d) und 
'\ der wahre Wert von X neben dem Fehler der 
\nnäherung eingetragen. — 

Sieht man vom Zeichenfehler ab. dann muß dieses 
Verfahren ebenso konvergieren wie das bei der 
Herleitung von Näherungswerten für X mit Hilfe 
der dureh die entsprechende Integralgleichung ver- 
schärften Näherung für die Eigenfunktion |3| durch 
len Vergleich der Ordinatenflächen. Für den Fall a) 
wurde rein rechnerisch X, = 4.183 gefunden. Der 
Unterschied zwischen diesem und dem auf gra- 
phischen Wege gefundenen analogen Wert ist 
1,79%, er ist auf die Ungenauigkeit der Zeichnung 
zurückzuführen. Beim zweiten Schritt ist dieser 
Unterschied 0,34%. Wie man sieht, ist die Aus- 
wirkung des Zeichenfehlers auf die Größe der 
Kippziffer gegenüber der Abweichung vom wahren 
Wert für die ersten Schritte nicht bedeutend. Bei 
einer oftmaligen Wiederholung des Verfahrens 
könnte es sich wegen des mit jedem Sehritt mit- 
reschleppten Zeichenfehlers von einem gewissen 
Sehritt an wohl herausstellen. daß die Näherungs- 
werte nicht mehr gegen die wahre Kippziffer kon- 
verzieren. Diese oftmalige Wiederholung ist aber 
für technische Zwecke nieht notwendig. Gewöhn- 
lieh dürfte bei zweckmäßiger Wahl von ®, (und 
evtl. X,) schon nach dem 2. oder 3. Schritt ein 
brauchbarer Näherungswert für X erreicht sein. 

Bei rein reehnerischer fortschreitender Verschär- 
funz von Näherungen für die Eigenfunktion be- 


stimmen sieh die Kippziffern am genauesten mit 
Hilfe der Energiebedingung |4|. So wird z.B. für 
den Kragträger mit zentrisch angreifender Endlast 
die wahre Kippziffer X = 4,0126 durch die Folge 
4.153, 4,0157, 4,012607 .... angenähert |3]. Die Be- 
nützung der graphisch gefundenen Näherung P,„ in 
der Energiegleichung erfordert aber mehr mit Ab- 
lese- und Zeichenfehlern behaftete Sehritte als die 
Gewinnung der Kippziffer aus und dureh 
Vergleich der Ordinatentlächen. Es bleibt noch zu 
untersuchen, unter welchen Bedingungen beim 
graphischen Verfahren die Benützung der Energie- 
gleiehung ebenfalls von Vorteil ist. 
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Dr. M. LAGALLY, o. Prof. an der Techn. Hoch- 
schule Dresden, Mechanik und Thermo- 
iynamik des stationären Glet- 
schers. (Sonderdruck aus Ergebnise der kosmi- 
schen Physik, Bd. IL) 80 S. m. 13 Fig. Leipzig 
1934. Akademische Verlagsees. m. b. H. Preis 
broseh. 7 RM. 

Der Verfasser ist seit über 20 Jahren in der 
Gletseherforschung tätig: er kennt die Gletscher- 
erscheinungen der Alpen von Grund aus und hat 
sich an der Vermessung der Ostalpengletscher er- 
folgreieh beteiligt. Wiederholt hat er im einschlägi- 
sen Schrifttum das Wort ergriffen: zuerst in be- 
schreibender und geschiechtlicher Art, später mit der 
Erörterung mechanischer und wärmekundlicher Auf- 
raben, aus deren Zusammenfassung das vorliegende 
Werk hervorgegangen ist. Vor mehr als 40 Jahren 
hat der Schreiber dieser Zeilen geometrische Ideen 
in die damals nenaufblühende Gletscherkunde der 
Ostalpen eingeführt. Es waren Fragen nach dem 
Zusammenhang (der Gletschersehwankungen mit der 
Witterung. die «damals die Aufmerksamkeit der 
Naturforscher auf die Gletscher lenkten und zur 
Beeriffsbildlung des „stationären“ Gletschers führ- 
ten. (dessen ortsfeste Strömung durch die Schwere 
bei gleichförmiger Schneezufuhr im Firngebiete und 
Abschmelzung im Zungengebiete unterhalten wird. 
Nimmt man zu diesem Begriff die Forderung des 
lückenlosen Zusammenhangs der bewegten Massen 
und der Stetigkeit der Bewegung hinzu. so ist der 
Ablauf der Bewegung bei bekannter Gestalt des 
(letschers und seines Bettes bereits in recht enge 
(renzen gebannt. Aus diesem Ablauf läßt sich eine 
Fülle teils wohlbekannter, aber auch bis dahin 
weniger  beachteter  Gletschererscheinungen her- 
leiten. Diese kinematische Gletscherlehre wird im 
vorliegenden Werk kurz auseinandergesetzt und 
dient als Grundlage für weitergehende dynamische 
und thermodynamische Betrachtungen. die tiefer in 
lie Ursachen des zeometrisch erkannten Be- 
weerungsablaufes hineinführen. Von den Navier- 
Stokesschen Differentialgleichungen für die Be- 


werune einer zähen Flüssigkeit ausgehend, nützt 
hagally ein von Somigliana angegebenes 
Integral, das für die stationäre Bewegung in einem 
zeneieten zylindrischen Bett mit gleichgeneigter 
ebener freier Oberfläche gilt, nach Möglichkeit aus. 
In (dieses Integral außer der Zähigkeitszahl 
eine Funktion einer komplexen Veränderlichen ein, 
die entweder dureh die Verteilung der Geschwindig- 
keit über die Oberfläche eines Querschnittes oder 
dureh die Form des Querschnittes selbst bestimmt 
werden kann und eine Abhängirkeit der oberfläch- 
lichen Geschwindigkeitsverteilune von der Quer- 
schnittsform vermittelt. So gehört z.B. zu einem halb- 
elliptischen (trogförmigen) Querschnitt eine parabel- 
förmige  Geschwindigkeitsverteilung mit einem 
Hochwert in der Mitte und annähernd linearem 
Abfall gegen die Ränder, während bei einem drei- 
eckigen Querschnitt die Geschwindigkeit vom Ranıl 
weg, wo sie Null ist, mit dem Quadrat des Ab- 
standes steigt. Kennt man an einem wirklichen 
stationären Gletscher sowohl die Verteilung der 
Oberflächengeschwindigkeit als auch die Quer- 
schnittslorm, so läbt sich die Zähigkeitszahl («des 
Gletschereises bestimmen. Durch die von A. Blümke 
und H. Heß vor 25 Jahren ausgeführten Tiefboh- 
rungen am Hintereisgletscher im Ötztal und die 
eleichzeitigen Geschwindiekeitsmessungen an der 
Gletscheroberfläehe sind «die Grundlagen für eine 
solche Bestimmung der Zähigkeitszahl gegeben. Der 
Wert. den Lagally hierfür auf «doppelte Weise 
findet, ist größer, als er nach Laboratoriumsver- 
suchen zu erwarten war, kann aber mit einigem 
Vertrauen auf die Ermittlung der Eismächtigkeit 
und Querschnittsform von Gletschern angewandt 
werden. von denen nur die Oberflächengeschwindig- 
keit bekannt ist. 


Viel schwieriger ist die Untersuchung der Spal- 
tenbildung der Gletscher. Dazu bedarf es der 
Kenntnis des Spannungszustandes des be- 
wegten Eises, wenn man von der naheliegenden 
Voraussetzung ausgeht, dab das Aufreißen der 
Spalten erfolgt, wenn die Zugspannungen die Festig- 
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keitsgrenze überschreiten. Lagally vergleicht beginnen mit der Berechnung des Energieumsatzes 


nun den Spannungszustand in einer zähen Flüssig- 
keit mit dem in einem elastischen Körper und ins- 
besondere die ähnlichgearteten Beziehungen, die in 
beiden Fällen zwischen den 6 Spannungskomponen- 
ten einerseits und den 6 Komponenten der Defor- 
mationszeschwindirekeiten bzw. der Verzerrung 
andererseits bestehen. Er betrachtet nun das Glet- 
schereis als einen elastischen Körper, in dem die 
Schubfestirekeit dauernd überschritten ist und in dem 
die Winkelverzerrunz mit der Zeit weiter anwächst. 
wobei ein schon von Maxwell gerebener Ansatz 
für plastische Körper zugrundegelegt wird. Der von 
der Zeit unabhängize elastische Anteil an der 
Winkelverzerrung wird durch den elastischen Schub- 
modul des Eises, die Geschwindigkeit des sich mit 
der Zeit ändernden Anteils durch die Zähigkeits- 
zahl des fließenden Eises geregelt. Werden Schub- 
modul und Zähierkeitszahl einander gleichgesetzt. so 
lassen sich bei geeigneter Wahl des elastischen Deh- 
nunesmoduls (die vorhin genannten Beziehungen 
eleichzeitie erfüllen. und zwar die elastischen für 
rasch ablaufende Vorgänge (z. B. Wellenbewegun- 
een), die plastischen für langdauernde, bei denen 
die kleine elastische Formänderung gegenüber der 
mit der Zeit eroßen Formänderung des Fließens 
vernachlässigt werden kann. Es wird so der Schub- 
modul des Gletschereises aus der Zähiekeitszahl be- 
stimmt. Der Dehnungsmodul kann auf diesem 
Wege nicht ermittelt werden, da beim Fließen des 
Gletschers das Eis als unzusammendrückbar voraus- 
eesetzt wird. während das elastische Verhalten eine, 
wenn auch sehr geringe Zusammendrückbarkeit vor- 
aussetzt. Man braucht da neue Versuchsmöglich- 
keiten, die sich aus der Beobachtunz der Ausbrei- 
tunz von künstlich dureh Zerknall von Spreng- 
körpern erzeugten Wellen ergeben. Es treten dabei 
L,äneswellen mit 3580 m und Querwellen mit 1670 m 
Sekundengesehwindirkeit auf. Aus diesen Fort- 
pflanzungszeschwindigekeiten rechnet sich der Deh- 
nungesmodul des Gletschereises im technischen Maß 
zu 673 ke/mm?, der Schubmodul zu 261 ke/mm?. 
Diese Zahlen verdienen als technische Mittelwerte 
für natürliches Gletschereis im  Schmelzzustande 
mehr Vertrauen als solche, die aus schwierig anzu- 
stellenden Laboratoriumsversuchen hervorgegangen 
sind. 

Nunmehr wird der Spannungszustand an der ebe- 
nen freien Oberfläche des Somiglianaschen 
Gletschers betrachtet. Eine Hauptspannung senk- 
recht zu der freien Oberfläche verschwindet. die 
beiden anderen lieren in derselben und schließen 
mit der Strömungsriehtunge den Winkel von 45° 
ein. Das gleiche muß dann auch von den aufreißen- 
den Spalten gelten, die entstehen. wenn eine jener 
Hauptspannungen eine Zugspannung ist. die die 
Festirkeit des Eises überschreitet. In der Tat zeigen 
(lie Gletscherspalten am Rande durchweg die von 
der Theorie zeforderte Eigentümlichkeit. von diesen 
wer unter 45° nach oben zu zu verlaufen. Unter 
Voraussetzung eines troeförmiren Querschnittes sind 
beim Somiglianaschen Gletscher die Randge- 
biete Zugspannungsgeobiete mit Überschreitung der 
Zugfestiekeit: es foleen zeren die Mitte zu Gebiete 
mit geringerer Zugspannung, während in (der Mitte 
selbst in allen 3 Hauptriehtungen Druckspannung 
und demnach Spaltenfreiheit herrscht. Die Ent- 
stehunz von Spalten «dureh Überschreitung der 
Schubfestigkeit des Eises wird von Lagally prak- 
tisch »brelehnt: für die an Gletschern beobachteten 
(leitflächen. auf welehe Philipp eine neue Be- 
werungstheorie der Gletscher aufzubauen gesucht 
hät, werden von ihm auf Gleitbewegungen längs 
Reißspalten zurückgeführt. Die wohldurchdachten 
Ausführungen über das Gletscherkorn und den inne- 
ren Aufbau des Gletschers seien hier übergangen. 
dagegen müssen die thermodynamischen 


Betrachtungen besonders hervorgehoben werden. Sie 


in einer strömenden zähen Flüssigkeit, deren Be- 
wegung dauernd so langsam ist, daß die lebendig» 
Kraft vernachlässigt werden kann, was für den 
(Gletscher in der Regel zutrifft. Der Energiesatz 
läßt sich dann so ausdrücken, daß die sekundliche 
Wärmeentwicklung in einem Raumteilchen des Eises 
gleich ist der sekundlich geleisteten Fallarbeit, plus 
der Arbeit des isotropen Druckes, plus der Rei- 
bungsarbeit, welch letztere sich aus der Arbeit der 
resultierenden Reibung und der eigentlichen Form- 
änderungsarbeit zusammensetzt. Für den ganzen 
(letscher ist die entwickelte Wärmemenge gleich 
der Fallarbeit, da wegen der Annahme, daß (die 
Bodengeschwindigkeit Null ist, die beiden anderen 
Arbeitsanteile im stationären Gletscher sich auf- 
heben: für einzelne Teile des Gletschers dageren 
nicht. Das wird sehr anschaulich für den Sonderfall 
des Somiglianaschen Gletschers gezeigt, dessen 
Querschnitt ein Rechteck mit langer wagrechter unil 
kurzer senkrechter Seite ist, und in dem die Strö- 
munesgeschwindigekeit parabelförmig von der Ober- 
fläche geeren die Tiefe zu auf Null abnimmt. Hier 
ist die Fallarbeit an der Oberfläche wegen der star- 
ken Geschwindiekeit eroß und nimmt nach der Tiefe 
zu ab, während bei der Wärmeentwicklung das Um- 
zekehrte zutrifft, da die Formänderung nach unten 
zunimmt. 

Hier setzt nun ein von Lagally zum ersten 
Male besonders betonter Begriff des „kalten“ Glet- 
schers ein. dessen Temperatur überall unter der je- 
weilieen. mit dem Druck abnehmenden Schmelz- 
temperatur des Eises liegt. In einem solchen bildet 
sieh infolge der Wärmeentwicklung im Innern ein 
Temperaturgefälle gegen die Oberfläche zu aus, wo- 
durch Wärme nach außen dureh Leitung abgeführt 
wird. Außerdem ist noch die Einfuhr von Wärme 
aus dem Erdinnern am Gletschergrunde zu be- 
achten, sowie jene aus der Außenluft. deren Wärme- 
gerad zwar nicht. wie es für einen stationären Glet- 
scher nötie wäre, gleichbleibend. vielmehr jährlichen 
und täglichen Schwankuneen unterworfen ist. Das 
thermische Verhalten eines solchen „kalten“ Glet- 
schers gleieht dann im ganzen dem eines Fels- 
bodens mit einer verminderten geothermischen Tie- 
fenstufe (20 m auf 1° ©). In den Polargegenden ist 
Vorkommen solcher kalter Gletscher von mäßiger 
Gesamtfallhöhe anzunehmen. Da jedoch die Fall- 
arbeit des Eises im Durchschnitt auf je 100 m Fali- 
höhe die Temperatur des Eises um 0.5° C erhöht. 
so nähert sich die Innentemperatur des Gletschers 
beim Abwärtsfließen vom Grunde her, wo Reibungs- 
wärme und Erdwärme zusammenwirken, der 
Schmelztemperatur des Eises unter dem jeweiligen 
Druck. Es wird dann bei großen Polargletschern 
wenigstens in der Tiefe der Zustand herrschen. der 
für die Gletscher der gemäßigten Zonen kennzeich- 
nend ist. Da die Zähierkeit des Eises unterhalb der 
Schmelztemperatur größer ist, muß bei großen 
Polargletschern außer dem Sprung im Temperatur- 
eefälle am Grunde der kalten Oberschicht des Glet- 
schers auch ein soleher in der Formänderung des 
Eises beim Fließen angenommen werden, der aber 
jenem entgegengesetzt ist. Es ist also in der kalten 
Oberschicht das Temperaturgefälle groß, gegenüber 
jenem der schmelzenden Unterschicht, in der es 
nur in geringem Maße vom Druckanstieg bedingt 
wird, während in der Oberschicht die Formänderung 
weren der größeren Zähirkeit des kalten Eises 
gering ist gegenüber jener in der Unterschicht mit 
veringer Zähirkeit. 

Das sind nur einige Stichproben aus dem Buche 
Lagallys, das sehr knapp, aber mit viel Um- 
sicht und genauer Kenntnis des einschlägigen 
Schrifttums abgefaßt ist, wovon schon das Schrif- 
re mit über hundert Nummern Zeugnis 
ablegt. 
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D. HILBERT und P. BERNAYS, Grundlagen 
\er Mathematik. Bd. 1. (Die Grundlagen d. 
math. Wiss. in Einzeldarstell. mit besonderer Be- 
‚icksiehtigung d. Anwendungsgeb. Hrsg. v. R. 
'ourant. Gemeinsam mit W. Blaschke, M. Born u. 
B. L. van der Waerden. Bd. 40.) Berlin 1934, 
‚Julius Springer Verlag. XII + 471 8. Preis 36 M. 


Insofern das Buch aus dem Rahmen konkreter und 
ınwendbarer mathematischer Forschung herausfällt, 
ist es für den Praktiker der Mathematik ohne Be- 
ang, insofern es sich darin um den Wahırheitswert 
ler Mathematik handelt, hat es für Mathematiker, 
Physiker und Philosophen, ja für jeden denkenden 
Menschen gleich hohes Interesse. — Die moderne an- 


vebliche Grundlagenkrise ist freilich in weiten Kreisen 


ihrer vermeintlichen Schwierigkeit und Nutzlosigkeit 
wegen übel beleumundet. Wer einseitig mitten in 
mathematischer Forschungsarbeit drinnen steht, 
neigt leicht zu der beschränkten Meinung, dab hier 
durch schlecht angebrachten Fleiß ein klarer Tat- 
bestand dunkel, verwirrt und unbrauchbar gemacht 
wird. In Wahrheit ist das Gegenteil der Fall: die 
Probleme sind nicht das Erzeugnis mathematischer 
Überkultur, sondern gehören seit Eudoxos zum 
ewigen Bestande menschlicher Erkenntnissorgen und 
treten bei der heutigen Erkenntnislage der Mathe- 
matik und Physik zwangsläufig in den Vordergrund. 
Zwei Schulen stehen heute einander gegenüber, die 
man durch die Schlagworte Formalismus und In- 
tuitionismus zu bezeichnen pflegt. Während der In- 
tuitionismus — kurz gesagt — nur solche Urteile in 
ler Mathematik und Physik gelten läßt, deren Sinn 
in der Anschauung vollziehbar ist. und die einen 
evidenten Tatbestand ausdrücken (Beispiel: „finite“ 
Mathematik, Zahlentheorie), erträgt der Formallsmus 
so enge Grenzen, die schon in der klassischen 
Mathematik, z. B. bei Einführung idealer, ohne 
weiteres nicht in der Anschauung aufweisbarer Ele- 
mente überschritten werden, nicht und dehnt die 
Erkenntnis auf transzendente, unabhängig von 
unserem Wissen um irgendwelche Phänomene er- 
schaffene Gegenstände aus (Beispiel: die transfiniten 
Ordnungszahlen der Gantorschen Mengenlehre). 
Während Tatbestände einander nicht widersprechen 
und die darauf bezüglichen Theorien zwar mensch- 
lichem Irrtum unterworfen, aber grundsätzlich 
widerspruchsfrei sind, ist für die transzendenten 
Theorien der Beweis ihrer Widerspruchslosigkeit zu 
ihrer sachlichen Rechtfertigung notwendig, haben 
sich doch bei der uferlosen Mengenbildung an den 
Grenzen der Mengenlehre Widersprüche ergeben. 
Der Beweis der Widerspruchslosigkeit wird dadurch 
seführt, daß das Verfahren eines mathematischen 
(finiten oder transfiniten) Beweises formalisiert wird 
und auf Grund der evidenten finiten Mathematik (in 
dieser Anwendung Metamathematik genannt) ge- 
zeigt wird, daß in dem formalisierten Beweisgang 
aus den Axiomenformeln durch den mathematischen 
Prozeß niemals die Formel des Widerspruchs ent- 
stehen kann, so wie beim Schachspiel aus der Aus- 
gangsstellung nach den Zugregeln niemals eine mehr 
als neun Damen enthaltende Stellung gewonnen 
werden kann. Es wird also gleichsam durch eine 
Art logisches Relais, die Formalisierung, die finite 
Mathematik zur Begründung der transfiniten be- 
nutzt. — Das von Bernays abgefaßte Buch hat als 
das derzeitige Standardwerk des Formalismus zu 
gelten. Auf die -angedeuteten philosophischen 
Hintergründe geht es fast nicht ein. Es ist im 
Gegensatz zu manchen anderen Veröffentlichungen 
über den Gegenstand lesbar und wohltuend ausführ- 
lich geschrieben. Doch wird nur der mit Erfolg es 
in die Hand nehmen, der eine gewisse Kenntnis der 
axiomatischen Methode der exakten Wissenschaften 
hat. Denn der Formalismus ist nichts anderes als 
die konsequente Vollendung der Axiomatik. 
Hierauf, sowie auf die Fortschritte näher einzu- 


gehen, die das Buch innerhalb der ausgedehnten 
Disziplin des Formalismus, insbesondere in der Be- 
weistheorie, erzielt, ist an dieser Stelie nicht der 
Ort. Es seien nur einige Paragraphen-Übersehriften 
angeführt: Das Problem der Widerspruchsfreiheit in 
der Axiomatik als logisches Entscheidungsproblem; 
die elementare Zahlentheorie; das finite Schließen 
und seine Grenzen; die Formalisierung des logischen 
Schließens in Aussagen- und Prädikatenkaikül: 
Widerspruchsfreiheit unendlicher Individuenbereiche:; 
die rekursiven Definitionen. — Noch sei hervor- 
gehoben, daß sowenig die Relativitätstheorie etwas 
mit philosophischem Relativismus zu tun hat, so- 
wenig der Formalismus mit dem sich berührt, was 
in neuester Zeit von einer journalistischen Ein- 
stellung her als toter Formelkram abgewertet wird. 
Bisweilen sagt man statt Formalismus treffender 
Symbolismus, da durch Anschauung nieht ohne 
weiteres adäquat erkennbare Gegenstände durch das 
Symbol der Formel dargestellt werden. 
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A. HEYTING, Mathematische Grund- 
lagenforschung, Intuitionismus, Be- 
weistheorie. (Ergebnisse der Mathematik und 
ihrer Grenzgebiete. Hrsg. v. d. Schriftleitung des 
„wentralblatt für Mathematik“, Ill. Bd., 4. Heft.) 
Berlin 1934, Julius Springer Verlag. X + 73 >. 

Verfasser, der zu den Intuitionisten (vgl. vor- 
stehendes Referat) zählt, gibt einen umfassenden, 
übersichtlichen enzyklopädischen Bericht über den 
derzeitigen Stand der Grundlagenfrage. Aus der 
Einleitung: „In diesem Referat behandle ich die 
mathematische Grundlagenforschung nach der intui- 
tionistischen und formalistischen Richtung. Es ist 
offenbar unmöglich, auf diesem Gebiet die philoso- 
phischen und die mathematischen Betrachtungen 
voneinander zu trennen. Im Gegenteil bin ich bei 
der Begrenzung des Stoffes von der engen Verbin- 
dung der Philosophie mit der Mathematik ausge- 
gangen. Nach der 'philosophischen Seite habe ich 
mich auf das für das Verständnis der matlemati- 
schen Probleme Notwendigste beschränkt; von den 
mathematischen Problemen behandle ich nur die- 
jenigen, die noch deutlich mit der Antwort auf die 
Frage ist Mathematik?" zusammenhängen.“ 

Dresden. W. Threlfall. 471 


Dr.-Ing. KURT BEYER, ord. Professor an der 
Technischen Hochschule Dresden. Die Statik 
im Eisenbetonbau. Ein Lehr- und Handbuch 
der Baustatik. Verfaßt im Auftrage des Deutschen 
Beton-Vereines. Zweite, vollständig neubearbeitete 
Auflage. Erster Band. Mit 572 Abbildungen im 
Text, zahlreichen Tabellen und Rechenvorschriften. 
Berlin 1933, Verlag von Julius Springer. VIII + 
389 S. Preis geb. 32.50 M. 

Der erste Band des Handbuchs der Baustatik um- 
faßt die statisch bestimmten Tragwerke, die Be- 
rechnung der Formänderung gerader und gekrümm- 
ter Stäbe und die statisch unbestimmten Tragwerke. 
Der Inhalt ist knapp formuliert und den heutigen 
Kenntnissen der Statik weitgehend angepaßt. Der 
reichhaltige Stoff, die vielen Tabellen und Literatur- 
angaben werden das Buch dem  rechnenden 
Ingenieur wertvoll machen !). 

Wien. J. Ratzersdorfer. 


Ferner sind bei der Sehriftleitung folgende 
Bücher eingegangen (ausführliche Besprechung 
bleibt vorbehalten): 

KARL H. GROSSMANN, Zürich, Elemente 
derelementaren Mechanik. 1. Teil, 99 S. 
Zürich 1934, Verlag Karl H. Großmann. Preis 
brosch. 3 Fr. 


1) Vgl. auch die Bespreehung der ersten Auflage, Bd. 8 
(1928), S. 153. 
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Dr. phil. HERMANN FRICKE, Öberiegierungsrat 
und Mitgl. des Reichspatentamts, Dieim Innern 
erdähnliche Sonne, eine neue Anschauung 
von Äther, Schwerkraft und Sonne, 738.—+1WAbb. 
Weimar 1934, Verlag R. Borkmann. Preis brosch. 
1,25 M. 


Dr. KARL GEY, Oberstudiendirektor in Leipzig, 
und Dr. HORST TEICHMANN, Privatdozeut a. d. 
Techn. Hochseh. Dresden, Einführung in die 
hehre vom Schuß (Ballistik, 
phys. Bibl. Reihe U, Bd. 11.) 113 S. + 50 Fig.. 
>2 Tafeln. Leipzig und Berlin 1934, Verlag B. 6. 
Teubner. Preis geb. 3.20 M. 


Dr. PH. LÖTZBEYER, Uberstudiendirektor in 
Berlin. Theorie und Praxis der Tafeln 


und des Tafelrecehnens. >. Dresden 
1934, Verlag L. Ehlermann. Preis geb. 3.20 M. 


KARL MARBE, Grundfragen der ange- 
wandten Wahrscheinlichkeitsrech- 
nung und theoretischen Statistik. 
177 S. München und Berlin 1934, Verlag U. H. 
secksche Verlagsbuchhandlung. Preis brosch. S M. 


Dr.-Ing. HUBERT SCHARDIN, Berlin, Das 
Toeplersche Schlierenverfahren. 
Grundlagen für seine Anwendung und 
quantitative Auswertung (Forschungsheit 
367, Beil. zu „Forschung auf dem Gebiete des Inge- 
nieurwesens"“, Ausg. B, Bd. 5, Juli/Aug. 1934). 
32 S. + 79 Abb., 5 Zahlentaf. Berlin 1934, VDI- 
Verlag G.m.b.H. Preis 5 M. für VDI-Mitgi. 4,50 M. 


NACHRICHTEN 


Tafelwerk höherer mathematischer Funktionen. 

Zur Zeit wird im Institut für angewandte Mathe- 
matik an der Universität Berlin im Rahmen der 
Wissenschaftlichen Akademikerhilfe der Notgemein- 
schaft der Deutschen Wissenschaft an der Berech- 
nung von Tafeln höherer mathematischer Funk- 
tionen gearbeitet. Es handelt sich, was den Kreis 
der Beschäftigten anbelangt, um Maßnahmen zur 
rhaltung der Arbeitsfähigkeit erwerbsloser Aka- 
demiker, denen Gelegenheit geboten werden soll, 
durch Arbeit auf ihren Fachgebieten ihre Kennt- 
nisse zu erhalten, zu erweitern und dadurch sich 
das Fortkommen in ihrem eigenen Berufe zu er- 
leichtern. Die Förderungssätze liegen über denen 
der Arbeitslosen-Unterstützung. Es sind noch 
Plätze frei. Interessenten können sieh im Institut 
für angewandte Mathematik an der Universität 
Berlin im Zimmer 256 oder Zimmer 88 melden. 

Es ist zunächst mit der Polynome der Hermite- 
schen und Laguerreschen Polynome und Ortho- 
onalfunktionen begonnen worden. Darüber hin- 
aus sollen Funktionen tabuliert werden, die nume- 
risch näher zu kennen das Bedürfnis besteht, für die 
aber bisher keine ausreichenden Tafeln vorliegen. 
Die Leitung «dieser Notarbeit, die in Händen (des 
Direktors des genannten Instituts liegt, ist gern be- 
reit, Anregungen und Wünsche in dieser Richtung 
zu berücksichtigen. A. Klose. 501 


Gesellschaft für angewandte Mathematik und 
Mechanik. 
Prager Mitglieder. 

Am 25. Oktober 1034 hielt Hr. Ing. Dr. techn. 
Alfred Eekert einen Vortrag über „Berechnung 
der Wärmestrahlung in der Heiztechnik“. 

Am 8. November 1934 sprach Hr. Ing. Dr. techn. 
Heinrich Treml über den „Anfangszustand bei 
Gebäudeschwingungen infolge von Erdbeben“. 


Persönliches 
Prof. Dr. Fr. A. Willers. o. Prof. der Mathe- 
matik an der Bergakademie Freiburg (Sachsen), 
wurde mit dem 1. Oktober 19534 von seinen amt- 
lichen Ptliehten entbunden (emeritiert). 


Prof./Dr. Wilhelm Müller, bisher Professor an 
der Deutschen Techn. Hochschule in Prag, hat einen 
Ruf als 0. Prof der Mechanik an die Techn. Hoclh- 
schule in Aachen erhalten und angenommen. 

Hr. Dr. Fritz Noether, bisher o. Prof. der 
Mechanik an der Techn. Hochschule Breslau, und 
Dr. Stefan Bergmann, früher Privatdozent in 
Berlin. sind ın das Forschungsinstitut für Mathe- 
matik und Mechanik an der Kujbyscheff-Universität 
Tomsk (USSR) berufen worden. Sie vertreten dort 
die Fächer der Funktionentheorie, mathematischen 
Physik und angewandten Mathematik. 

Der n. b. a. o. Prof. Dr. Alfred Klose in Berlin 
ist beauftragt worden, in der philosophischen Fa- 
kultät der Universität Berlin im Wintersemester 
1934/35 die bisher von Ministerialdirektor Prof. Dr. 
Vahlen innegehabte Professur für angewandte Ma- 
thematik sowie die Leitung des Instituts vertre- 
tungsweise zu übernehmen. 

Frau Dr. Hilda Pollaczek-Geiringer, bis- 
her an der Universität Brüssel, früher in Berlin, ist 
als außerordentlicher Professor (Charge de cours) 
an die Universität Istanbul berufen worden. 

Hr. Dr. Hans Fromm, früher Privatdozent uni 
seit September 1934 a. 0. Prof. an der Techn. 
Hochschule Berlin. hat einen Ruf auf den Lehrstuhl 
für Baustoffkunde und Festiekeitslehre an der 
Techn. Hochschule Danzig unter Ernennung zum 
o. Professor erhalten und angenommen. 


Der von Domherrn. Dr. Dr. Ing. Alfred Acker- 
mann-Teubner. im Jahre 1922 der Universität 
Leipzig gestiftete „Alfred-Ackermann-Teubner-Ge- 
dächtnispreis zur Förderung der mathematischen 
Wissenschaften” ist für das Jahr 1934 durch das 
Preisgericht dem Professor Dr. Dr. ing. E. h. Erich 
Trefftz an der Technischen Hochschule in Dresden 
für seine Arbeiten auf dem Gebiete der angewandten 
Mechanik zuerkannt worden. 

Die Universität Brüssel hat bei Gelegenheit ihrer 
100-Jahr-Feier Hrn. Prof. Dr. P.Debve. Direktor 
des Physikalischen Instituts der Universität Leipzig. 
und Hrn. Prof. Dr. R.v. Mises., früher in Berlin. 
jetzt in Istanbul, den Titel eines Dr. h.e. verlieben. 


Einbanddecken für den Jahrgang 1934. 


Erst dureh das Einbinden wird der beendete Jahrgang zu einem handlichen und übersichtlichen 
Nachschlagewerk. so daß man seinen Inhalt jederzeit bequem auswerten kann. Wir haben deshalb auch 
für den Jahrgang 1934 der Zeitschrift für angewandte Mathematik und Mechanik wieder Einbanddecken 
herstellen lassen, die zum Preise von 2.25 RM (für VDI-Mitglieder 2.—) durch jede Buchhandlung 


bezogen werden können. 
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